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ANALYSIS  OF  SEASAT  SAR  IMAGERY  COLLECTED 
DURING  THE  JASIN  EXPERIMENT 


1 

INTRODUCTION 

This  report  presents  results  from  a  continuing  study  for  the 
Office  of  Naval  Research  of  synthetic  aperture  radar  imagery  col¬ 
lected  by  the  Seasat  satellite  (see  Kasischke,  1980;  Kasischke  and 
Shuchman,  1981;  Vesecky,  et  al.,  1982,  1982a).  This  study  investi¬ 
gated  in  detail  Seasat  SAR  imagery  collected  over  the  North  Atlantic 
Ocean  during  the  Joint  Air-Sea  Interaction  (JASIN)  experiment  during 
the  summer/fall  of  1978.  The  results  of  an  associated  JASIN/Seasat 
study,  an  investigation  of  the  use  of  Seasat  SAR  data  to  estimate 
surface  wind  speed,  are  presented  in  a  separate  document  (see  Lyden, 
et  al.,  1983). 

This  study  considered  the  Seasat  SAR's  ability  to  detect  and 
measure  ocean  gravity  waves  and  its  potential  for  imaging  large, 
mesoscale  surface  features  such  as  internal  waves,  frontal  boundaries 
and  weather  related  patterns.  This  report,  which  discusses  the  vari¬ 
ous  analyses  performed  during  the  study,  is  organized  into  five  chap¬ 
ters,  including  this  introduction,  plus  four  appendices. 

Chapter  2  discusses  the  background  of  the  JASIN  experiment  and 
briefly  reviews  the  Seasat  SAR  system.  Chapter  3  discusses  the  addi¬ 
tional  analyses  performed  on  the  Seasat  SAR  wave  imagery  beyond  those 
conducted  during  the  initial  study  reported  by  Kasischke  (1980). 
This  chapter  includes  sections  on  correlations  between  sea  truth  and 
occasions  when  the  Seasat  SAR  imaged  gravity  waves,  the  testing  of 
SAR  image  enhancement  routines  on  Seasat  data,  and  assessment  of  a 
technique  to  extract  wave  height  information  from  SAR  imagery.  Chap¬ 
ter  4  discusses  the  Seasat  SAR's  ability  to  detect  three  mesoscale 
surface  features  present  in  the  JASIN  test  area:  internal  waves. 


frontal  boundaries,  and  surface  features  associated  with  meteorologi¬ 
cal  phenomena.  Finally,  Chapter  5  presents  a  summary  of  this 
investigation. 

Appendix  A  to  this  report  presents  diagrams  outlining  the  loca¬ 
tion  of  each  Seasat  SAR  pass  used  during  this  study.  Each  diagram 
also  summarizes  the  wind  and  wave  conditions  present  on  that  date. 
Appendix  B  presents  the  fast  Fourier  transforms  which  were  generated 
from  the  JASIN  Seasat  SAR  imagery.  Appendix  C  presents  a  paper  by 
Thomas  (1982)  which  proposes  a  technique  to  extract  wave  height  in¬ 
formation  from  Seasat  SAR  data.  This  technique  is  evaluated  in  Chap¬ 
ter  3.  Appendix  D  to  this  report  contains  copies  of  articles  which 
resulted  from  the  ONR-sponsored  research  of  the  Seasat  SAR  data  col¬ 
lected  during  JASIN. 
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2 

BACKGROUND 


This  chapter  briefly  reviews  the  Joint  Air-Sea  Interaction  exper¬ 
iment  and  the  synthetic  aperture  radar  onboard  the  Seasat  satellite. 
For  a  more  in-depth  discussion  of  the  JASIN  experiment,  see  Allan 
and  Guymer  (1980).  For  a  review  of  the  Seasat  SAR's  capabilities  to 
image  oceanographic  phenomena,  see  Beal,  et  al .  (1981),  Shuchman,  et 
al.  (1981)  or  Fu  and  Holt  (1982). 

2.1  THE  JOINT  AIR-SEA  INTERACTION  EXPERIMENT 

Data  used  in  this  study  were  collected  during  the  Joint  Air-Sea 
Interaction  (JASIN)  experiment.  The  JASIN  experiment  was  proposed 
in  1966  as  part  of  the  Global  Atmospheric  Research  Project  and  was 
sponsored  and  organized  by  the  Royal  Society  of  the  United  Kingdom. 
The  primary  aims  of  the  JASIN  experiment  were  to  "observe  and  distin¬ 
guish  between  the  physical  processes  causing  mixing  in  the  atmo¬ 
spheric  and  oceanic  boundary  layers  and  relate  them  to  mean  proper¬ 
ties  of  the  layers;  and  to  examine  and  quantify  aspects  of  the  momen¬ 
tum  and  heat  budgets  in  the  atmospheric  and  oceanic  boundary  layers 
and  fluxes  across  and  between  them"  (Royal  Society,  1979).  To  this 
end,  fifty  principal  investigators  from  nine  countries  collaborated 
on  an  intensive  experiment  from  mid-July  to  mid-September,  1978. 
Fourteen  oceanographic  research  vessels  were  on  site  at  one  time  or 
another  during  this  period.  The  JASIN  experiment  took  place  in  the 
North  Atlantic,  in  an  area  just  west  of  Scotland,  as  is  shown  in 
Figure  1.  Most  of  the  meteorologic  and  oceanographic  observations 
were  made  within  the  Oceanographic  Intensive  Area  (OIA)  (the  large, 
"dashed-line"  circle  in  Figure  1),  but  ships  operating  in  the  large 
triangular  region  made  observations  as  well.  Wave  data  were  obtained 
from  visual  observations  from  these  ships,  from  wave  recorders 
mounted  on  the  ships,  from  wave  buoys  deployed  by  the  ships,  and  from 
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the  returned  echoes  are  processed  to  determine  the  distance  to  the 
scene  and  to  obtain  the  measure  of  relative  strength  of  the  return 
(aQ).  In  real-aperture  systems,  the  resolution  in  azimuth  is  lim¬ 
ited  by  the  beamwidth  (i.e.,  the  physical  length)  of  the  antenna. 
Since  resolution  decreases  with  range  for  such  systems,  operation  at 
spacecraft  altitudes  is  not  possible  with  real-aperture  radars. 

Synthetic  aperture  radar  is  a  coherent  radar  that  uses  the  motion 
of  a  moderately  broad  physical  antenna  beam  to  synthesize  a  very 
narrow  beam,  thus  providing  fine  azimuthal  (along-track )  resolution 
(Brown  and  Porcello,  1969;  Harger,  1970).  Fine  range  (cross-track) 
resolution  is  achieved  by  transmitting  either  very  short  pulses  or 
longer  coded  pulses  which  are  compressed  by  matched-filtering  tech¬ 
niques  into  equivalent  short  pulses.  Usually,  the  coded  pulse  is  a 
waveform  linearly  modulated  in  frequency. 

Historically,  in  SAR  systems,  the  phase  history  of  a  scattering 
point  in  the  scene  is  recorded  on  photographic  film  as  an  anamorphic 
(astigmatic)  Fresnel  zone  plate.  The  parameters  of  the  zone  plate 
are  set  in  the  azimuth  direction  by  the  Doppler  frequencies  produced 
by  the  relative  motion  between  the  sensor  and  the  point  scatterer, 
and  in  the  range  direction  by  the  structure  of  the  transmitted 
pulses.  The  film  image  is  a  collection  of  superimposed  zone  plates 
representing  the  collection  of  point  scatterers  in  the  scene.  This 
film  is  used  by  a  coherent  optical  processor  (a  set  of  spherical  and 
cylindrical  Fourier  transform  lenses)  which  focuses  the  anamorphic 
zone  plates  into  the  points  which  produced  the  microwave  scatter  of 
the  scene  (Kozma,  et  al . ,  1972). 

Recently,  SAR  systems  have  employed  digital  techniques  to  both 
record  and  process  the  data.  Digital  processing  typically  uses 
matched  filtering  techniques  to  "dechirp"  the  signal  in  range,  while 
Fast  Fourier  Transform  (FFT)  techniques  achieve  the  required  azi¬ 
muthal  compression  of  the  SAR  Doppler  history  (Ausherman,  19c.pi. 
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permanently  anchored  wave  rider  buoys.  The  wave  rider  buoys  were 
anchored  within  the  OIA. 

Seasat,  launched  just  prior  to  the  JASIN  experiment,  offered 
investigators  a  unique  opportunity  to  collect  satellite  data  over  a 
well  instrumented  test  site.  A  Seasat  SAR  receiving  station  had  been 
set  up  at  Oakhangar,  England  for  the  European  Space  Agency  and  radar 
data  were  collected  on  twenty  different  occasions  when  Seasat  passed 
over  or  near  the  JASIN  study  area.  The  geographical  relationship 
(to  the  JASIN  area)  of  the  eighteen  Seasat  passes  examined  for  the 
JASIN  study  are  presented  in  Appendix  A  of  this  report.  Data  from 
the  other  two  passes  were  of  too  poor  quality  to  process. 

2.2  THE  SEASAT  SYNTHETIC  APERTURE  RADAR 

The  data  presented  in  this  report  were  collected  by  the  Seasat 
satellite.  Among  the  instrumentation  carried  by  Seasat,  which  was 
launched  during  June  of  1978,  was  a  synthetic  aperture  radar  (SAR). 
This  satellite  collected  over  500  passes  of  SAR  data  before  suffering 
a  catastrophic  power  loss  in  October  of  1978.  The  SAR  on  board 
Seasat  was  an  L-band  (23.5  cm  wavelength)  radar.  It  collected 
25  x  25  meter  resolution  imagery  with  a  ground  swath-width  of  100  km 
and  a  swath-length  of  up  to  4000  km,  and  viewed  the  surface  of  the 
earth  with  an  average  incidence  angle  of  20°. 

An  imaging  radar  such  as  the  Seasat  SAR  is  an  active  device  that 
senses  the  environment  with  short  wavelength  electromagnetic  waves. 
As  an  active  sensor,  the  Seasat  SAR  provided  its  own  illumination  in 
the  microwave  region  of  the  electromagnetic  spectrum  and  thus  was 
not  affected  by  diurnal  changes  in  emitted  or  reflected  radiation 
from  the  earth's  surface. 

Side-looking  radars  such  as  the  Seasat  SAR  use  pulse  ranging, 
whereby  the  radar  antenna  attached  to  a  moving  platform  illuminates 
a  section  of  the  ground  and  the  amplitude,  phase  and  polarization  of 


TABLE  1 

RADAR  PARAMETERS  OF  THE  SEASAT  SYNTHETIC  APERTURE  RADAR  (SAR) 


Frequency 

Wavelength 

Transmitted  Bandwidth 
Pulse  Duration 

Pulse  Time-Bandwidth  Product 
Transmitter  RF  Power 
Transmitter  Type 
PRF 

Satellite  Altitude 

Nominal  Range  (20°) 

Antenna  Beamwidth,  Elevation 

Antenna  Beamwidth,  Azimuth 

Antenna  Pointing  Angle 

Surface  Resolution 

Slant  Range  Resolution 
Azimuth  Resolution  (one  look) 

Integration  Time: 

Image  Swath  Width 
Image  Length 
Sensor  Power 
Satellite  Velocity 
Satellite  Altitude 
Nominal  Range  (20°) 


1274.8  MHz 

23.5  cm  (L-band) 

19  MHz 

33.8  ysec 
642 

800  W  Peak  -  46  W  Average 
Solid-State  Bipolar  Transmitter 
1647 
~800  km 
~850  km 
6° 

1° 

19-25°  off  nadir,  right  side 

25  m  x  25  m  (4  look  data) 

8  m 

6.25  m 

0.5  to  2.5  sec  depending  on 
resolution 

100  km 

250  to  4000  km 

60  W,  nominal  operation 

-7800  m/sec 

“800  km 

~850  km 
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After  appropriate  system  corrections,  digital  data  have  better  geo¬ 
metric  and  radiometric  fidelity  than  optical  data  and  are  subject  to 
fewer  image  processing  artifacts  such  as  noise  streaks.  Addition¬ 
ally,  digital  data  recorded  in  square  root  mode  can  store  40  dB  worth 
of  information.  In  contrast,  an  optical  SAR  image  has  only  an  18  dB 
dynamic  range.  A  disadvantage  of  digital  SAR  data  is  that  processing 
algorithms  are  quite  slow,  particularly  when  one-look  (full  resolu¬ 
tion)  data  are  required.  Due  to  long  processing  times  (i.e.,  8  hours 
for  a  100  x  100  km  Seasat  scene),  a  limited  amount  of  digital  Seasat 
SAR  data  exists.  The  optical  processing  time  for  a  100  x  100  km 
scene  is  approximately  one  minute,  excluding  set-up  or  film 
development. 

Specific  parameters  of  the  Seasat  SAR  system  are  presented  in 
Table  1.  The  wavelength  at  which  the  Seasat  SAR  system  operates, 
23.5  cm,  does  not  penetrate  more  than  a  few  centimeters  into  the 
water  surface.  Therefore,  the  recorded  signal  is  primarily  reflected 
from  the  water's  surface  and  is  a  function  of  the  surface  roughness, 
slope,  and  motion  of  the  scatterers.  The  23.5  cm  wavelength  will 
also  determine  the  ocean  wavelength  (Bragg  wavel  to  which  the  SAR  is 
most  sensitive.  The  Bragg  ocean  wavelength  for  resonance  for  the 
Seasat  SAR  parameters  is  approximately  36  cm.  The  bandwidth  and  in¬ 
cidence  angle,  given  in  Table  1,  determine  the  ground  range  resolu¬ 
tion.  The  range  to  target,  velocity,  azimuth  resolution  and  antenna 
dimensions  are  information  that  can  be  utilized  to  determine  the  syn¬ 
thetic  aperture  length.  Tf  this  information  indicates  an  azimuth 
resolution  of  6.25  me  *  >"ound  requires  a  synthetic  aperture 
length  of  approximately  1 "  asat  traveling  at  800  km  altitude 

had  an  orbital  speed  of  7.9  ,  and  therefore  required  about  two 

seconds  to  synthesize  the  requited  aperture. 
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The  initial  analysis  of  the  Seasat  SAR  imagery  collected  during 
the  JASIN  experiment  focused  on  two  areas  (see  Kasischke,  1980; 
Kasischke  and  Shuchman,  1981;  Vesecky,  et  al.,  1982,  1982a): 

(1)  comparison  of  dominant  wavelength  and  direction  of  gravity  waves 
derived  from  SAR  data  to  those  obtained  from  surface  measurement  de¬ 
vices;  and  (2)  analysis  of  the  conditions  under  which  the  Seasat  SAR 
imaged  surface  gravity  waves,  using  a  wave  contrast  measurement 
(peak-to-back ground  ratio  or  PBR). 

Since  then,  additional  analyses  have  been  performed  on  the  JASIN 
Seasat  SAR  wave  data: 

1.  Generation  of  digital  fast  Fourier  transforms  from  all  18 
JASIN  passes, 

2.  Updating  of  surface  truth  information  and  further  analysis 
of  the  wave  contrast  measurements, 

3.  Testing  of  motion  compensation  algorithms  on  Seasat  SAR 
images  of  gravity  waves,  and 

4.  Evaluation  of  a  technique  published  by  Thomas  (1982)  to  ex¬ 
tract  wave  height  estimates  from  Seasat  SAR  data. 

This  chapter  discusses  these  activities. 

3.1  FAST  FOURIER  TRANSFORM  ANALYSIS 

Previous  comparisons  of  dominant  wavelengths  and  directions  of 
gravity  waves  derived  from  Seasat  SAR  data  to  surface  truth  estimates 
were  based  largely  on  the  use  of  large  area  optical  Fourier  trans¬ 
forms  (OFTs).  Although  it  was  demonstrated  during  these  earlier 
studies  that  OFT  estimates  of  dominant  wavelength  and  direction  com¬ 
pare  favorably  to  surface  measured  estimates,  there  are  several 
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i  shortcomings  in  this  technique.  First,  using  OFTs  requires  manual 

measurement  of  the  distance  to  the  center  of  the  OFT  spectrum  and 
the  spectrum  orientation.  Thus,  this  technique  introduces  the  varia¬ 
bility  of  the  human  interpreter.  Next,  to  generate  the  OFTs  from 
I  the  JASIN  Seasat  data,  very  large  apertures  were  used  (25  by  25  km). 

In  areas  where  waves  are  undergoing  transformations,  such  as  refrac¬ 
tion  (in  shoaling  water)  or  diffraction  (around  islands),  it  is  not 
always  practical  to  use  such  large  apertures.  Finally,  in  the  near 
future,  nearly  all  SAR  data  will  be  digitally  processed.  To  produce 
OFTs,  it  will  be  necessary  to  first  generate  an  image  from  the  digi¬ 
tal  data,  adding  an  additional  step  to  the  analysis. 

The  purpose  of  this  portion  of  the  study  was  to  generate  fast 
Fourier  transforms  of  the  Seasat  SAR  data  collected  during  the  JASIN 
experiment. 

3.1.1  METHODS 

All  digital  data  discussed  in  this  report  were  generated  by 
ERIM's  hybrid  optical-digital  processor  (Ausherman,  et  al . ,  1975). 
The  SAR  image  films  (optically  processed  at  ERIM)  of  ocean  gravity 
waves  from  the  JASIN  Seasat  passes  were  digitized  using  a  12.5  meter 
by  12.5  meter  pixel  aperture.  A  1024  by  1024  line-by-point  area  was 
generated  and  the  data  recorded  on  a  computer  compatible  tape  (CCT). 
The  area  selected  for  digitization  was  the  portion  of  the  pass  which 
was  closest  to  the  position  where  the  concurrent  surface  gravity  wave 
measurements  were  collected  (see  Appendix  A). 

The  digital  SAR  ocean  wave  imagery  was  then  analyzed  using  a 
standard  computer  package  developed  by  ERIM  under  sponsorship  of  the 
Office  of  Naval  Research  (see  Shuchman,  et  al . ,  1979).  Using  this 
routine,  a  512  by  512  pixel  subset  of  the  digital  data  is  extracted 
from  the  CCT.  Usually  the  SAR  data  are  geometrically  rectified,  but 
since  the  ERIM  data  had  been  previously  geometrically  corrected  when 


the  data  were  optically  correlated,  this  step  was  not  necessary. 
The  data  are  then  corrected  for  long-period  variations  (such  as  power 
loss  or  antenna  gain)  to  remove  D.C.  bias,  smoothed  using  a  (sin  x)/x 
filter  to  reduce  radar  speckle  and  resampled  to  256  by  256  pixels. 
This  resampling  results  in  a  pixel  size  of  25  by  25  meters.  These 
data  are  then  inputted  to  a  fast  Fourier  transform  program  to  gener¬ 
ate  spectral  estimates  from  the  SAR  data. 

When  the  digital  SAR  data  are  fast  Fourier  transformed,  the  re¬ 
sults  are  typically  displayed  as  two-dimensional  contour  plots.  Pre¬ 
sented  in  Figure  2  is  a  typical  two-dimensional  contour  plot  produced 
from  an  FFT  of  SAR  ocean  wave  data  and  one-dimensional  spectra 
derived  from  it.  Each  of  the  three  contour  levels  in  Figure  2a  rep¬ 
resents  3  dB  of  energy  from  the  SAR  image.  These  three  contour 
levels  range  from  -3  dB  to  -9  dB  (i.e.,  one-half  to  one-eighth  of 
the  peak  value).  Although  this  FFT  produces  no  wave  height  informa¬ 
tion,  the  two-dimensional  contour  plot  does  portray  wave  number  and 
wave  direction  information.  The  x-axis  of  the  plot  represents  azi¬ 
muth  direction  and  the  y-axis  the  range  direction,  with  the  units  on 
the  axes  being  wave  number  (k).  By  finding  the  center  of  the  highest 
contour  level  in  the  spectrum,  and  measuring  k  and 

kazimuth»  an  estimate  °f  the  dominant  wave  number  (and  hence  wave¬ 
length)  can  be  calculated,  as  can  an  estimate  of  the  wave's  orienta¬ 
tion.  If  the  SAR  platform's  direction  is  known,  an  estimate  of  the 
dominant  wave  direction  can  be  calculated.  Note,  there  is  still  a 
180’  ambiguity  in  the  wave  direction  data. 

The  information  in  the  two-dimensional  FFT  spectrum  can  be  more 
rigorously  analyzed.  To  do  so,  one-half  of  the  plot  in  Figure  2a  is 
summed  digitally  over  specified  wave  number  ranges  to  produce  a  one¬ 
dimensional  wave  number  spectrum,  as  is  shown  in  Figure  2b.  Finally, 
a  one-dimensional  wave  direction  spectrum  (usually  at  the  peak  wave 
number)  can  be  plotted,  as  is  presented  in  Figure  2c.  The  digital 
output  from  these  plots  can  be  consulted  to  extract  exact  estimates 
of  dominant  wavelength  and  direction. 


Fast  Fourier  Transforms 


FIGURE  2 


Wave  Direction  Spectra 


Wave  Number  Spectra 


EXAMPLE  OF  TWO-DIMENSIONAL  FFT  AND  RESULTANT 
ONE-DIMENSIONAL  WAVE  DIRECTION  AND  WAVE  NUMBER 
SPECTRA.  (Seasat  Rev.  1087  Data.) 


During  the  JASIN  experiment,  surface  truth  (wave  data)  were  col¬ 
lected  by  a  variety  of  different  instruments.  This  instrumentation 
included  pitch  and  roll  buoys  (Stewart,  1977)  deployed  from  the  R/V 
Atlantis  II  and  R/V  Discovery,  a  set  of  moored  wave  rider  buoys  which 
were  within  the  Oceanographic  Intensive  Area,  wave  recorders  onboard 
the  R/V  John  Murray  and  R/V  Discovery  and  visual  readings  made  from 
the  R/V  John  Murray.  The  exact  location  of  the  research  vessels 
varied  day  to  day  depending  on  what  experiments  they  were  participat¬ 
ing  in,  but  all  obtained  wave  data  at  the  time  of  the  Seasat  over¬ 
pass,  or  at  worst,  within  an  hour  or  two  of  it.  The  wave  rider  buoys 
were  all  located  within  the  Oceanographic  Intensive  Area  (their  exact 
locations  are  marked  by  Figure  1)  and  the  data  from  these  were  col¬ 
lected  at  the  times  of  che  Seasat  overpasses. 

Wave  data  and  sea-truth  estimates  were  processed  using  standard 
oceanographic  techniques  at  the  Scripps  Institute  of  Oceanography. 
Included  in  this  sea-truth  data  set  were  estimates  of  dominant  wave¬ 
length  and  direction,  wave  height  (H^3)  and  surface  wind  speed 
and  direction.  The  two-dimensional  ocean  wave  frequency  spectra, 
F(e,a>)  obtained  from  the  buoys  and  wave  recorders  were  converted  to 

O  O 

wave  number  spectra  using  the  Jacobian  approximation  {q  lu , 
where  g  is  the  gravitational  acceleration  and  id  is  the  wave  frequen¬ 
cy)  appropriate  to  deep  water  waves  (Vesecky,  et  al.,  1981). 

Because  several  of  the  test  dates  had  more  than  one  available 
set  of  sea-truth,  it  was  necessary  to  decide  which  value  to  use  to 
compare  to  the  FFT  value.  Available  sea-truth  was  therefore  ranked 
according  to  the  desirability  of  the  data  source,  as  follows: 

BEST  1.  Atlantis  II  Pitch  and  Roll  Buoy  or  Discovery  Pitch  and 
Roll  Buoy. 

2.  Moored  OIA  Wave  Rider  Buoy. 

3.  John  Murray  Wave  Recorder  or  Discovery  Wave  Recorder. 
WORST  4.  John  Murray  Visual  Readings. 
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If  no  sea-truth  data  were  available  using  the  most  desirable  cate¬ 
gory,  then  the  next  highest  category  which  contained  the  desired 
sea-truth  data  was  used. 

3.1.2  RESULTS 

The  two-dimensional  fast  Fourier  transforms  and  their  resultant 
one-dimensional  wave-number  and  wave  direction  spectra  (when  appro¬ 
priate)  are  presented  in  Appendix  B  to  this  report. 

Table  2  summarizes  the  coincident  surface  measured  wave  and  wind 
data  along  with  the  SAR-derived  estimates  of  dominant  wavelength  and 
direction  and  wave  contrast  for  the  eighteen  SAR  overpasses  made 
during  the  JASIN  experiment.  Included  in  Table  2  are  the  Seasat  pass 
number,  the  date  and  time  the  SAR  data  were  collected,  the  source  of 
the  sea  truth  used  to  compare  to  the  SAR  data,  the  distance  between 
the  SAR  and  surface  observations,  the  surface  measured  values  for 

dominant  wave  direction  (©j),  dominant  wavelength  (x),  significant 
wave  height  ( / 3) ,  wind  speed  (U)  and  wind  direction  (oj)  and 

SAR  derived  estimates  of  wave  direction  (e^._0FT  and  eT-FFT), 
wavelength  (x-OFT  and  x-FFT)  and  wave  contrast  (PBR).  (See 
Kasischke,  1980  for  a  discussion  on  wave  contrast.) 

From  Table  2,  we  can  see  that  estimates  of  dominant  wavelength 
and  direction  were  obtained  from  Seasat  SAR  data  from  thirteen  passes 
using  OFTs,  and  from  seven  passes  using  FFTs.  Figure  3  plots  the 

SAR  derived  (FFT)  estimates  of  wavelength  and  direction  versus  the 

surface  measured  values. 

For  wavelength,  it  can  be  seen  that  both  the  estimates  derived 
from  FFTs  and  OFTs  tend  to  be  longer  than  the  sea  truth  values,  with 
the  FFT  values  being  somewhat  longer.  There  are  two  cases  where  the 
OFT  derived  value  varies  significantly  from  the  surface  measured 
value  of  direction  of  wave  propagation.  Other  than  these  two  cases, 
there  is  very  little  variation.  From  Table  2  and  Figure  3,  it  c*.n 
be  seen  that  the  FFTs  and  OFTs  give  essentially  the  same  result. 
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_  Expected  Line 

Slope  3  1.0 
Intercept  *  0 

£ _ A  OFT  Regression  Line 

~  Slope  3  1 .04 

Intercept  3  21  m 
R  3  0.98 

r\ _ rA  FFT  Regression  Line 

°  ^  Slope  3  0.97 

Intercept  3  23  m 

R  3  0.98 
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DIRECTION  OF  WAVE  PROPAGATION  COMPARISON 


—  Expected  Line 
Slope  3  1.0 
Intercept  3  0 

.  A0pT  Regression  Line 
Slope  =  0.79 
Intercept  3  42  m 
R  3  0.91 

.  .0  FFT  Regression  Line 
Slope  3  0.63 
Intercept  3  93  m 
R  3  0.74 


00  360.00 


FIGURE  3.  SEA  TRUTH  VERSUS  SAR  DERIVED  ESTIMATES  OF  DOMINANT  WAVELENGTH 
AND  DIRECTION  OF  PROPAGATION  FOR  JASIN  DATA. 
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It  is  interesting  to  note  that  the  OFT  technique  provided  spec¬ 
tral  wave  estimates  on  data  for  which  the  FFT  technique  showed  no 
dominant  wave.  Recall  that  an  aperture  of  25  by  25  kilometers  was 
used  for  each  OFT,  while  an  aperture  of  6.4  by  6.4  kilometers  was 
used  for  each  FFT.  Thus,  the  aperture  for  the  OFTs  was  over  15  times 
as  large  as  the  aperture  for  FFTs.  Hence,  the  OFTs  has  a  better 
spectral  resolution  than  the  FFTs  and  therefore  could  in  some  cases 
resolve  wave  spectra  features  (dominant  wavelength)  and  direction 
where  the  FFT  could  not.  This  large  of  an  aperture  can  be  used  in 
areas  where  little  change  in  wavelength  and  direction  are  occurring, 
but  would  probably  result  in  very  diffuse  spectra  in  areas  of  wave 
refraction  or  diffraction.  Also,  because  of  the  computation  time 
involved,  most  FFT  routines  used  today  are  limited  to  areas  of  256 
by  256  pixels. 

3.2  SAR  DETECTION  OF  GRAVITY  WAVES 

In  determining  whether  a  spaceborne  SAR  system  will  be  an  effec¬ 
tive  tool  for  gathering  wave  data,  an  area  of  important  consideration 
is  understanding  the  conditions  under  which  a  SAR  will  and  will  not 
image  gravity  waves.  The  JASIN  data  set  provides  a  valuable  oppor¬ 
tunity  to  examine  this  question  because,  not  only  were  a  wide  variety 
of  surface  measurements  made  which  can  be  compared  to  the  SAR  data, 
but  also  because  a  wide  spectrum  of  environmental  conditions  were 
encountered  during  the  time  span  of  the  JASIN  study. 

3.2.1  METHODS 

Both  qualitative  and  quantitative  techniques  were  applied  to  look 
for  correlations  between  SAR  imaging  of  gravity  waves  and  environ¬ 
mental  parameters.  These  parameters  were  either  directly  measured 
as  part  of  the  surface  truth,  or  derived  from  the  surface  truth  using 
theoretical  relationships. 


As  mentioned  in  the  previous  section,  wave  measuring  instruments 
collected  wave  period  and  directional  data  at  or  near  the  areas  where 
SAR  data  were  collected.  From  these  wave  data,  significant  wave 
heights  (^3)  were  also  calculated.  Other  surface  measurements 
of  potential  interest  are  surface  wind  speed  (U)  and  direction  (a). 

When  examining  directional  data  (wind  or  wave),  it  is  important 
to  know  how  the  waves  (or  wind)  are  traveling  with  respect  to  the 
radar  look  direction.  The  following  directional  variables  were 
cal cul ated: 

1.  Normalized  wave  direction.  When  ©A  =  0°,  the  waves 
are  traveling  in  the  same  direction  as  the  SAR  plat¬ 
form;  when  =  180°,  the  opposite  direction 

2.  Wave  direction  normalized  to  a  0°  to  90°  range  (where 
0°  is  an  azimuth  traveling  wave  and  90°  is  a  range 
travel  ing  wave) 


3.  a$:  Wind  direction  relative  to  the  radar  look  direction 

(normalized  to  a  0°  to  90°  range,  where  90°  would  be 
a  wind  blowing  perpendicular  to  the  SAR  look 
direction) 

4.  a,,:  Wind  direction  relative  to  wave  direction  (normalized 

to  a  0°  to  90°  range,  where  90°  would  be  a  wind  blow¬ 
ing  perpendicular  to  the  wave  direction  of  travel). 

It  has  been  theorized  that  factors  which  might  affect  SAR  imaging 
of  gravity  waves  include  the  local  tilt  of  the  water  surface  (Elachi 
and  Brown,  1977)  and  the  orbital  velocity  of  the  gravity  waves  which 
causes  a  bunching  of  the  smaller  capillary  and  ultra-gravity  waves 
(Alpers  and  Rufenach,  1979;  Phillips,  1981). 

Wave  slope  (S)  can  be  defined  as: 

S  =  H1/3/x  (1 ) 

where  x  =  the  dominant  wavelength. 
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(v0)  as 


Kinsman  (1965)  defines  the  orbital  velocity  of  a  gravity  wave 

“  H 1  /  3 


vo 


(2) 


where  o>  =  wave  frequency  ( 2irg /x)^. 


A  normalized  orbital  velocity  (v  ),  defined  as  the  orbital 
velocity  in  the  direction  of  radar  1 ine-of-sight,  can  be  calculated 
as: 


v  =  v  sin 
no  AC 


(3) 


The  surface  measurements  and  derived  variables  were  compared  to 
the  Seasat  SAR  data  in  both  a  qualitative  and  quantitative  manner. 
For  the  qualitative  comparison,  the  SAR  data  were  subjectively 
divided  into  three  visibility  categories,  and  common  denominators 
within  these  categories  for  surface  measurements  were  determined. 
For  the  quantitative  analysis,  a  wave  crest-to-trough  contrast  mea¬ 
surement  was  used  as  a  numerical  indicator  of  visibility  of  waves  on 
SAR  imagery  and  statistically  correlated  to  the  various  surface  mea¬ 
surements  and  derived  variables. 


The  SAR  imagery  collected  over  the  JASIN  area  was  visually  exam¬ 
ined  and  classified  as  having  waves  in  one  of  three  categories: 
(1)  GOOD  imagery,  where  gravity  waves  are  clearly  visible  on  the  SAR 
image;  (2)  POOR  imagery,  where  gravity  waves  are  present,  but  de¬ 
tected  only  upon  close  inspection;  and  (3)  NONE  imagery,  where  no 
gravity  wavelike  patterns  are  visible  on  the  SAR  image.  Upon  cate¬ 
gorizing  the  imagery,  the  surface  measurements  and  derived  data  for 
each  category  were  examined  to  determine  if  a  trend  could  be 
detected. 


The  data  set  was  then  more  rigorously  analyzed  by  statistical!'/ 
comparing  the  surface  measurements  and  derived  variables  to  a  recent¬ 
ly  developed  wave  crest-to-trough  contrast  measurement  known  as  a 
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peak-to-background  ratio  (Kasischke,  1980;  Kasischke  and  Shuchman, 
1981).  The  peak-to-background  ratio  (PBR)  is  measured  as  follows: 
First,  the  peak  intensity  of  the  Fourier  transform  spectrum  is  ex¬ 
tracted  (either  measured  from  an  optical  Fourier  transform  or  digi¬ 
tally  derived  from  a  fast  Fourier  transform).  Next,  a  minimum  in¬ 
tensity  from  the  same  k-space  as  the  peak  intensity  is  extracted. 
The  ratio  of  these  two  intensities  (i.e.,  peak/minimum)  is  defined 
as  the  peak-to-background  ratio.  For  this  portion  of  the  study,  the 
peak  and  minimum  intensities  were  measured  as  film  densities  of  the 
photographic  negative  of  the  spectrum  obtained  from  an  optical 
Fourier  transform  of  the  wave  data. 

The  peak-to-background  ratios  for  the  eighteen  JASIN  SAR  over¬ 
passes  were  linearly  correlated  (Draper  and  Smith,  1966)  to  the  vari¬ 
ous  surface  measurements  and  derived  variables  (such  as  wave  slope, 
orbital  velocity  and  normalized  orbital  velocity)  to  determine  where 
a  relationship  might  exist.  It  should  be  noted  that  a  linear  corre¬ 
lation  does  not  imply  that  a  cause/effect  relationship  exists  between 
two  variables  in  an  a  posteriori  examination  of  the  data  such  as  is 
presented  here;  it  only  suggests  where  such  a  relationship  might 
exist. 

3.2.2  RESULTS 

Table  3  lists  wavelength  (x),  wave  height  (H^3)  and  radar 
look  direction  normalized  to  the  direction  of  wave  travel  (e^) 
for  the  three  classes  of  SAR  wave  imagery  (i.e.,  GOOD,  POOR  and 
NONE).  For  e^,  a  0*  reading  means  the  waves  are  traveling  paral¬ 
lel  to  the  radar  1 ine-of-si ght  and  a  90°  reading  indicates  the  waves 
are  traveling  perpendicular  to  the  radar  line-of-sight.  An  examina¬ 
tion  of  Table  3  reveals  that  one  of  two  conditions  existed  when  waves 
were  not  visible  on  the  SAR  imagery:  (1)  the  wave  height 
was  less  than  1.3  meters;  or  (2)  the  waves  were  traveling  nearly 
parallel  to  the  radar  line-of-sight  (&AC  <  35°).  For  the  POOR 


20 


RADAR  DIVISION 


TABLE 

3 

SUMMARY  OF 

WAVELENGTH,  HEIGHT, 

AND  DIRECTION 

DATA  FOR  JASIN  DATA 

★ 

Revolution 

\ 

Hl/3 

eAC 

None 

590 

89-222 

1.3 

84  “-6° 

599 

71 

1.1 

15° 

633 

89 

1.1 

77° 

642 

169 

2.9 

15° 

1006 

105 

3.6 

33° 

Poor 

556 

149 

1.5 

64° 

714 

182* 

— 

59°* 

719 

164* 

— 

75°* 

834 

149 

2.8 

29° 

958 

100 

1.3 

41° 

Good 

547 

169 

1.7 

58° 

757 

215 

4.9 

25" 

762 

206 

4.3 

90° 

791 

132 

3.3 

85° 

1001 

63 

2.5 

73° 

1044 

301* 

3.5 

44° 

1049 

301 

3.5 

59° 

1087 

299* 

32° 

*A11  data  use  sea  truth 

estimates 

except 

for  those  denoted  with 

an  (*),  where  SAR  estimates 

were  used 

• 
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wave  imagery,  we  can  also  see  that  one  of  the  above  limits  was  being 
approached  on  three  out  of  the  five  cases. 


Table  4  summarizes  the  derived  variables  which  were  calculated 
to  supplement  the  surface  measured  data.  Included  in  this  table  is 
wave  slope,  orbital  velocity  (v  ),  normalized  orbital  velocity 
(v  ),  two  estimates  of  wave  direction  relative  to  the  SAR  look 
direction  (e^-sea  truth  and  e^-OFT),  wave  direction  normalized 
to  radar  look  direction  (e^-sea  truth  and  e^-OFT)  and  two  nor¬ 
malized  wind  directions  (a<-^  and  y) .  The  surface  measured 
and  derived  variables  were  then  linearly  correlated  to  the  wave  con¬ 
trast  peak-to-background  ratio  (PBR). 


Table  5  summarizes  the  linear  correlations  between  the  peak-to- 
background  ratio  and  the  various  surface  measurements  and  derived 
variables.  Significant  (o  <  0.05)  linear  correlations  existed  be¬ 
tween  the  wave  contrast  measurement  and  wavelength  (both  x-sea  truth 
and  x-OFT),  and  wave  height  (H^).  Table  6  summarizes  multiple 
linear  regression  coefficients  calculated  between  combinations  of  x, 
H-|  12  and  sin  vq.  Again,  these  coefficients  are  all  signif¬ 
icant  (0  <  0.05). 


From  the  analysis  of  the  sea  and  wind  conditions  during  the  JASIN 
study,  there  appears  to  be  a  positive  correlation  between  wave  de¬ 
tectability  on  SAR  imagery  and  both  wave  height  and  wavelength  of 
gravity  waves.  In  other  words,  the  longer  the  wavelength  and  the 
higher  the  wave,  the  greater  the  chance  of  detecting  that  wave  on 
imagery  collected  by  a  Seasat  type  SAR. 


3.3  SAR  MOTION  CORRECTIONS 

Because  the  imaging  mechanism  employed  by  the  synthetic  aperture 
radar  is  dependent  upon  the  phase  history  of  the  target  being  imaaed, 
any  motion  in  that  target  may  result  in  a  degraded  SAR  imaae.  This 
phenomena  was  first  reported  on  by  Raney  (1971)  and  more  recent! 7 


TABLE  5 

LINEAR  CORRELATIONS  BETWEEN  WAVE  CONTRAST  (PBR)  AND 
ENVIRONMENTAL  PARAMETERS  FOR  JASIN  SEASAT  DATA 


Parameter 

Linear 

Correlation 

Wave  Height 

0.53** 

Wavelength  (Sea  Truth) 

0.67** 

Wavelength  (OFT) 

0.71* 

Wave  Slope 

-0.02 

Wave  Direction  Relative  to 
Radar  Look  Direction  (Sea 
Truth) 

0.28 

Wave  Direction  Relative  to 
Radar  Look  Direction  (OFT) 

-0.05 

Normalized  Orbital  Velocity 

0.43 

Wind  Direction  Relative  to 

SAR  Look  Direction 

-0.06 

Wind  Direction  Relative  to 
Wave  Direction 

-0.10 

Wind  Speed 

0.30 

*Significant  at  o  <  O.Ol 

**Signif icant  at  o  <  0.05 
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TABLE  6 

MULTIPLE  LINEAR  CORRELATION  COEFFICIENTS 


Parameter 

PBR 

H 1  /  3  *  x 

* 

0.71 

* 

x,  sin  ©AC  vO 

0.82 

x,  sin  eAC  Vq,  H1/3 

0.83 

*Significant  at  o  <0.01 
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Shuchman  and  Zelenka  (1978)  and  Shuchman  (1981).  Shuchman  (1981) 
produced  a  set  of  equations  which  can  be  used  to  calculate  the 
adjustments  necessary  to  compensate  for  these  motion  distortions  when 
processing  the  SAR  signal  histories. 

Studies  using  aircraft  SAR  data  have  shown  that  the  visibility 
or  detectability  of  gravity  waves  is  often  sensitive  to  motion  com¬ 
pensation  adjustments  made  during  the  processing  of  the  SAR  signal 
histories  (Kasischke,  et  al.,  1979;  Kasischke  and  Shuchman,  1981; 
Shuchman,  et  al.,  1982).  Because  these  motion  compensation  adjust¬ 
ments  are  inversely  proportional  to  the  velocity  of  the  SAR  platform, 
it  was  theoretically  determined  by  early  Seasat  investigators  that 
the  motion  adjustments  may  be  necessary  for  optimum  contrast  wave 
imagery  if  the  SAR  data  are  processed  to  full  resolution.  Further¬ 
more,  it  was  determined  by  Shuchman  and  Zelenka  (1978)  that  the 
adjustment,  if  necessary,  would  be  quite  small  and  the  effect  on  the 
SAR  wave  imagery  quite  subtle. 

3.3.1  METHODS 

The  purpose  of  the  present  investigation  was  to  determine  the 
degree  of  sensitivity  of  Seasat  SAR  data  to  motion  compensation 
adjustments.  This  was  accomplished  in  two  separate  experiments. 
First,  wave  detectability  or  contrast  was  measured  as  a  function  of 
both  the  range  telerotation  adjustment  and  azimuth  teleshift  adjust¬ 
ment.  Next,  SAR  imagery  with  various  telerotation  and  teleshift 
adjustments  were  digitized  and  fast  Fourier  transformed  to  determine 
if  the  spectral  estimates  varied. 

A  target  velocity  in  the  range  (1 ine-of-sight)  direction  affects 
the  SAR  imaging  process  in  several  ways.  One  effect  is  an  azimuthal 
displacement  of  the  moving  target's  image  relative  to  a  stationary 
target's  image.  When  a  target  is  accelerating  in  the  range  direc¬ 
tion,  this  azimuth  displacement  changes  during  the  imaging  time. 
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resulting  in  a  smearing  in  the  azimuthal  direction.  Neither  of  these 
effects  is  correctable  during  processing. 

A  range  velocity  component  will  also  cause  a  smearing  or  defocus- 
ing  in  the  range  direction  due  to  a  rotation  of  the  phase  history  of 
the  target.  This  effect  can  be  corrected  for  during  correlation  of 
the  data  by  rotating  the  lenses  in  the  optical  processor  as  described 
below. 

Motion  of  a  SAR  imaged  target  in  the  along  track  direction  re¬ 
sults  in  a  defocusing  of  the  image  in  the  azimuth  direction.  This 
defocusing  can  be  compensated  for  during  processing  by  adjusting  the 
focal  length  of  the  azimuthal  lens.  For  a  more  detailed  discussion 
of  the  problems  associated  with  imaging  moving  targets  with  a  SAR, 
the  reader  is  referred  to  articles  by  Shuchman  (1981)  or  Shuchman, 
et  al.  (1981). 

Shuchman  (1981)  presented  equations  to  calculate  adjustments  to 
use  during  processing  of  the  signal  histories  for  the  range  and  azi¬ 
muth  velocity  components  present  in  gravity  waves.  These  velocity 
adjustments  are  directly  proportional  to  the  velocity  of  the  SAR 
platform.  The  motion  compensation  adjustments  are  of  two  types: 
telerotation  adjustments  which  compensate  for  motion  in  the  range 
(1 ine-of-sight)  direction  and  teleshifts  which  compensate  for  motion 
in  the  azimuth  direction. 

The  telerotation  adjustment  (Shuchman,  1981)  is: 


VAC  Q 


where 


vr  is  the  1 ine-of-sight  velocity  of  the  target, 
vA£  is  the  SAR  platform  velocity,  and 
P  and  Q  are  SAR  processing  parameters. 


In  a  similar  fashion,  an  adjustment  can  also  be  made  to  compen¬ 
sate  for  azimuth  velocity  distortions  in  the  SAR  data.  As  derived 
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by  Shuchman  and  Zelenka  (1978),  the  focal  length  (F  )  for  SAR 
processing  is  defined  as: 


'\T~ ) 
\  AC/ 


where  R  is  the  slant  range  to  the  target, 
x  is  the  radar  wavelength, 

XQ  is  the  optical  processor  wavelength, 

M  is  azimuth  demagnification  of  the  optical  processor,  and 
is  the  film  velocity. 

Using  the  focal  length  (FQ)  calculated  in  Eq.  (5),  the  shift 
in  focal  length  (5F)  produced  by  a  moving  target  is  defined  as 

VA 

6F  =  2F  (6) 

0  VAC 

where  vft  is  the  target  velocity  parallel  to  the  SAR  platform 
direction. 

It  can  be  seen  that  both  adjustments  are  inversely  proportional 
to  the  velocity  of  the  SAR  platform.  For  the  velocity  of  the  Seasat 
satellite  (~7800  m/s),  the  adjustments  needed  to  correct  for  the 
phase  velocity  of  the  gravity  waves  are  very  much  smaller  than  those 
needed  for  an  aircraft  SAR  (v^^  =  150  m/s). 

Telerotation  adjustments  were  made  on  Seasat  SAR  data  from  Revo¬ 
lution  762.  Surface  measurements  indicated  the  dominant  wavelength 
at  this  time  was  210  meters  (phase  velocity  =  C  =  18.1  m/s)  and  the 
SAR  data  indicated  a  241  meter  (C  =  19.4  m/s)  dominant  wave  was 
present.  The  telerotation  adjustment  (0)  calculated  for  the  Revolu¬ 
tion  762  SAR  data  assumed  a  19.4  m/s  dominant  wave  was  present.  SAR 
imagery  of  waves  from  Revolution  762  were  processed  using  the  follow¬ 
ing  telerotation  adjustments:  -9 0,  -70,  -50,  -30,  -20,  -0,  0  +  0, 
+20,  +30,  +50,  +70,  and  +90.  A  negative  (-)  telerotation  assumes 
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the  waves  are  moving  away  from  the  SAR  and  a  positive  (+)  telerota¬ 
tion  assumes  the  waves  are  moving  towards  the  SAR.  During  Rev.  762, 
the  waves  were  moving  towards  the  SAR. 

Seasat  SAR  imagery  from  Revolution  1087  was  optically  processed 
using  a  variety  of  different  azimuth  focus  shifts.  No  sea-truth  was 
collected  for  this  pass,  but  it  was  clear  from  an  examination  of  the 
SAR  imagery  that  the  waves  are  nearly  azimuth  traveling.  The  optical 
Fourier  transform  of  the  SAR  data  indicates  a  299  meter  dominant  wave 
is  present  (C  =  21.2  m/sec),  with  an  orientation  of  34°  relative  to 
the  satellite  ground  track.  A  teleshift,  P,  was  calculated  assuming 
the  waves  had  a  phase  speed  of  21.2  m/sec.  Revolution  1087  data  were 
processed  using  the  following  teleshifts:  -4P,  -3P,  -2P,  -P,  0,  P, 
2P,  3P,  4P.  A  negative  (-)  focus  shift  assumes  the  waves  are  moving 
in  the  same  direction  as  the  SAR  platform  while  a  positive  (  +  )  cor¬ 
rection  assumes  the  waves  are  moving  in  the  opposite  direction. 
Since  no  sea  truth  was  available,  the  assumption  was  made  that  the 
waves  were  propagating  towards  the  Scottish  Coast,  therefore  opposite 
the  SAR  platform  direction. 

To  explore  the  question  of  whether  or  not  the  telerotation  and 
teleshift  adjustments  significantly  alter  the  SAR  derived  estimates 
of  wavelength  and  direction,  digital  analysis  techniques  were  em¬ 
ployed.  SAR  imagery  from  both  Revolutions  762  and  1087  were  digi¬ 
tized  using  ERIM's  Optical  Hybrid  Digitizing  Facility  (see  Ausherman, 
et  al.,  1975).  For  Revolution  762,  imagery  from  the  -2 <t>,  -<t>,  0,  +<t> 
and  +20  telerotation  settings  were  digitized.  For  Revolution  1087, 
imagery  from  the  -2P,  -P,  0,  +P  and  +2P  azimuth  teleshifts  were  dig¬ 
itized.  These  digitial  data  were  then  fast  Fourier  transformed  and 
the  resultant  one-dimensional  wave  number  and  wave  direction  plots 
examined. 
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3.3.2  RESULTS 


Figure  4  and  Table  7  summarize  the  results  of  the  wave  contrast 
versus  telerotation  measurements.  Note  on  Figure  4  that  a  positive 
telerotation  is  compensating  for  wave  motion  away  from  the  radar 
while  a  negative  telerotation  adjusts  for  motion  towards  the  radar. 
Each  point  on  the  curve  represents  an  average  of  five  separate 
measurements. 

The  trend  in  Figure  4  is  clear.  Telerotation  adjustments  do  not 
appear  to  appreciably  degrade  Seasat  SAR  imagery  within  ±30  of  the 
stationary  focus  setting.  Studies  of  aircraft  SAR  data  have  shown 
that  waves  become  undetectable  on  SAR  imagery  when  a  telerotation  of 
±2 0  from  the  optimum  setting  is  used  (see  Kasischke,  et  al.,  1979; 
Kasischke  and  Shuchman,  1981). 

Figure  5  illustrates  the  two-dimensional  spectra  generated  by  a 
fast  Fourier  transform  of  the  digitized  SAR  data  from  Revolution  762* 
for  three  telerotation  adjustments  (-2 0,  0,  +20).  Figure  5  also 
presents  the  one-dimensional  wave  direction  and  wave  number  spectra 
generated  from  the  fast  Fourier  transform.  Table  8  summarizes  the 
dominant  wavelength  and  direction  of  propagation  for  the  various 
adjustments.  We  can  see  that  the  spectral  estimates  do  not  change 
when  a  non-zero  range  telerotation  is  used. 

Table  7  and  Figure  6  summarize  the  results  of  the  wave  contrast 
versus  the  azimuth  teleshift  adjustments.  Again,  each  point  repre¬ 
sents  an  average  of  five  points.  It  can  be  seen  that  the  highest 
wave  contrast  was  obtained  when  a  land  (0)  teleshift  setting  was 
used.  With  the  azimuth  teleshifts,  it  appears  that  the  Seasat  SAR 
wave  imagery  was  not  significantly  degraded  until  a  focus  shift 
greater  than  ±2P  was  used. 

Figure  7  illustrates  the  two-dimensional  spectra  generated  from 
the  Revolution  1087  data.  For  the  three  azimuth  focus  settinos  (-2°, 
0,  +2P),  Figure  7  also  presents  the  one-dimensional  wave  direction 
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FIGURE  4.  PEAK-TO-BACKGROUND  RATIO  VERSUS  TELEROTATION  ADJUSTMENT.  (Seasat  Rev.  762.) 


TABLE  8 

FAST  FOURIER  TRANSFORM  ESTIMATES  OF  DOMINANT 
WAVELENGTH  AND  DIRECTION  AS  A  FUNCTION  OF 
RANGE  TELEROTATION  ADJUSTMENTS  FOR 
SEASAT  REV.  762 
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DIMENSIONAL  FAST  FOURIER  TRANSFORMS  AND  ONE-DIMENSIONAL  WAVE  DIRECTION  AND 
NUMBER  SPECTRA  VERSUS  TELEROTATION  ADJUSTMENTS  FOR  SEASAT  REV.  7R 7  DATA 


Standard  Deviation 


FIGURE  6.  PEAK-TO-BACKGROUND  RATIO  VERSUS  AZIMUTH 
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and  wave  number  spectra  from  these  data.  Table  9  summarizes  the 
dominant  wavelength  and  direction  of  propagation  for  the  various 
teleshift  adjustments.  From  Figure  7  and  Table  9,  we  can  see  that 
the  spectral  estimates  generated  using  different  azimuth  focus  set¬ 
tings  change  very  little  when  non-stati onary  settings  are  used. 

It  appears  that  for  the  wave  conditions  present  during  Revolu¬ 
tions  762  and  1087,  motion  compensation  adjustments  do  not  apprecia¬ 
bly  improve  the  detectability  of  waves  on  the  Seasat  SAR  imagery. 
This  insensitivity  to  motion  effects  is  true  in  terms  of  both  detect¬ 
ability  of  waves  on  SAR  imagery  and  the  estimates  of  dominant  wave¬ 
length  and  direction  obtained  from  the  SAR  data.  This  is  not  always 
the  case  when  processing  data  from  air craft -mo unted  SARs  (see 
Kasischke,  et  al.,  1979;  Kasischke  and  Shuchman,  1981;  or  Shuchman, 
et  al . ,  1982)  and  may  not  be  the  case  for  all  ocean  waves  imaged  by 
Seasat. 

3.4  WAVE  HEIGHT  DETERMINATION  USING  DIGITALLY  PROCESSED  SEASAT  SAR 
IMAGERY 

This  section  of  the  report  examines  a  technique  proposed  by 
Thomas  (1982)  for  the  extraction  of  wave  height  information  from 
digitally  processed  Seasat  SAR  imagery  of  ocean  waves.  A  copy  of 
Thomas'  paper  is  included  as  Appendix  C.  The  discussion  presented 
below  will  include  a  brief  review  of  the  theory,  measurements  and 
results  presented  by  Thomas,  as  well  as  some  measurements  made  at 
ERIM  which  tend  to  confirm  his  results.  Finally,  a  brief  discussion 
of  the  limitations  and  possible  errors  associated  with  this  approach 
is  presented. 

3.4.1  THEORETICAL  DEVELOPMENT 

The  three  basic  assumptions  in  Thomas'  method  include:  (1)  that 
any  wave  modulation  in  the  SAR  imagery  is  caused  by  tilt  modulation 
effects,  (2)  that  all  image  intensity  variations  in  a  wave  imaoe  are 


TWO-DIMENSIONAL  FAST  FOURIER  TRANSFORMS  ANO  ONE-DIMENSIONAL  WAVE  DIRECTION  AND 
WAVE  NUMBER  SPECTRA  VERSUS  TELEROTATION  ADJUSTMENTS  FOR  SEASAT  REV.  1087  DATA 
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the  result  of  the  imaged  gravity  waves  (due  to  tilt  modulation)  or 
speckle,  and  (3)  that  the  speckle  characteristics  of  a  digitally 
processed  SAR  image  are  known.  Each  of  these  assumptions  will  be 
addressed  briefly  in  Section  3.4.2.  The  method  of  Thomas  relates 
the  contrast  in  a  digitally  processed  SAR  wave  image  to  a  change  in 
local  incidence  angle  (slope)  using  a  theoretical  scattering  model. 
By  combining  this  slope  measurement  with  wavelength  (easily  obtained 
by  FFT  analysis),  an  estimate  for  wave  height  can  be  calculated. 
Presented  below  is  a  review  of  this  process  adopting  Thomas' 
nomenclature. 

Denoting  the  ocean  surface  displacement  from  the  mean  level  by 
f(r,  t),  the  significant  wave  height,  H  can  be  expressed  as 

Hs  =  4  ,  (7) 

where  the  brackets  indicate  a  mean  value.  For  a  narrow-band  swell- 
wave  spectrum  centered  on  wave  number  kQ,  the  r.m.s.  slope  of  the 
sea  surface  is  given  by 

k  H  irH 

tan  er  =  — ^ —  =  Y\~’  (8) 

o 

where  xQ  is  the  swell  wavelength.  Therefore,  H$  can  be  deter¬ 
mined  if  tan  and  x„  are  known.  Determination  of  x„  is 

r  o  o 

easily  accomplished  by  taking  a  Fast  Fourier  Transform  (FFT).  The 
r.m.s.  slope  (tan  ©  )  of  the  ocean  surface  is  obtained  by  compar¬ 
ing  the  variation  of  the  radar  cross-section  attributed  to  waves  with 
a  theoretical  scattering  model  which  gives  the  variation  of  radar 
cross-section  as  a  function  of  incidence  angle.  The  theoretical 
variation  of  radar  cross-section  as  a  function  of  incidence  angle 
was  calculated  by  Valenzuela  (1978),  and  is  shown  in  Figure  8.  Al¬ 
though  the  data  in  Figure  8  are  for  428  MHz  and  Seasat  operated  at 
1.275  GHz  (L-band),  theory  predicts  similar  results  for  these  two 


FIGURE  8.  COMPARISON  OF  MEASURED  AND  THEORETICAL  CROSS-SECTIONS  OF 
THE  OCEANS  FOR  A  RADAR  FREQUENCY  OF  428  MHz,  HORIZONTAL 
POLARIZATION.  (After  Valenzuela,  1978.) 
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frequencies.  The  variation  in  radar  cross-section  attributed  to 
waves  can  be  related  to  a  change  in  local  incidence  angle  (r.m.s. 
slope)  by  using  the  slope  of  the  curve  in  Figure  8  at  that  particular 
incidence  angle. 

A  SAR  image  contains  intensity  variations  due  to  both  speckle 
and  natural  reflectivity  changes.  Speckle  is  a  characteristic  of 
coherent  imaging  systems  and  is  described  by  Porcello,  et  al .  (1976). 

Thomas  presents  a  somewhat  confusing  argument  where  he  refers  to 
the  term  I  or  the  speckle  intensity.  An  alternative  way  of 
deriving  his  results  is  to  describe  the  effects  of  speckle  by  a  mul¬ 
tiplicative  model  (Wu,  1980;  Lee,  1981)  which  may  be  written  as 

I  =  IS  (9) 


where  I  is  the  intensity  at  a  given  point,  I  is  the  mean  intensity 
in  the  neighborhood  of  that  point,  and  S  is  a  random  variable  with 
mean  1.0  and  standard  deviation 


o 


S 


(10) 


where  N  is  the  number  of  independent  looks  in  the  processor 
(Porcello,  et  al.,  1976).  We  may  also  view  I  as  the  "true"  signal, 
un contaminated  by  speckle.  For  the  case  in  point,  the  true  signal 
is  due  to  waves  on  the  surface,  so  we  change  notation  and  replace  I 
with  I  (the  subscript  referring  to  waves).  We  also  replace  I 
with  It  to  avoid  confusion.  Equation  (9)  then  becomes 

It  =  I«s  .  <") 

from  which  we  also  have 

ll  -  I2S2  .  (12) 

We  now  invoke  a  theorem  from  probability  theory  which  states  that 
the  mean  value  of  the  product  of  two  uncorrelated  variables  is  equal 
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to  the  product  of  the  mean  values  of  those  variables.  Applying  this 
to  Eqs.  (11)  and  (12),  we  have 

<it>  =  <iw><s>  03) 

and 

<It>  =  <!w><s2>  *  O4) 

Using  the  definition  of  the  variance,  we  have 

<!t>2  =  at  +  <!t>2  (15) 

9  9 

and  similarly  for  cO  and  o^.  Substutiting  these  into  Eq. 
(14)  and  dividing  by  the  square  of  Eq.  (13)  results  in 


which  is  the  same  as  Thomas'  Eq.  (11).  From  Eq.  (12)  and  <S>  =  1, 
the  last  quantity  in  brackets  is  equal  to  (1  +  1/N),  where  N  is  the 
number  of  looks.  Since  the  JPL  and  DFVLR  digital  images  are  proc¬ 
essed  to  four  looks,  Eq.  (16)  becomes 


which  is  Eq.  (4)  in  Thomas  (1982). 

2 

The  terms  and  <I-t>2  can  be  obtained  easily  from  the 

digital  image,  and  a  /< I  >  can  then  be  calculated  from  Ea.  (17). 

Cl)  CO 

To  relate  this  to  Valenzuela's  curve  (Figure  8),  the  quantity  V  =  10 

1  og-j q  [(<IU>  +  )  /( < I  >  -  a^)]  is  calculated.  This 

allows  the  determination  of  tan  er  which  leads  to  the  significant 
wave  height. 
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3.4.2  RESULTS 

Presented  in  Table  10  are  the  results  of  Thomas'  two  test  cases 
and  three  additional  measurements  made  at  ERIM  following  the  methods 
of  Thomas.  In  general,  the  wave  height  estimates  agree  relatively 
well  with  those  reported  from  sea  truth  measurements.  Thomas  does 
not  elaborate  on  his  actual  image  measurement  technique.  For  the 
analysis  presented  here,  image  statistics  from  JPL  digitally  proc¬ 
essed  SAR  data  for  areas  of  100  x  100  pixels  (~1.6  x  1.6  km)  were 
generated.  The  size  of  the  area  selected  to  implement  Thomas'  algo¬ 
rithm  could  affect  the  results.  If  the  area  chosen  is  too  small  to 
accurately  characterize  the  contrast  in  the  SAR  wave  scene,  the  wave 
height  estimate  will  be  in  error.  An  example  would  be  if  only  a 
portion  of  a  wavelength  was  studied.  In  this  case,  the  contrast  and 
resultant  wave  height  estimate  would  both  be  low.  Also,  if  there 
were  not  a  sufficient  number  of  points  to  "average"  out  a  few  out¬ 
liers,  the  contrast  measurement  and  resultant  wave  height  estimate 
would  be  too  high.  On  the  other  hand,  if  the  area  for  statistics 
generation  is  too  large,  the  incidence  angle  range  it  would  subtend 
would  not  allow  Valenzuela's  curve  to  be  used  in  the  manner  presented 
in  Section  3.4.1 . 

Thomas  makes  the  assumption  that  all  wave  modulation  in  the  SAR 
wave  scene  is  due  to  tilt  modulation.  This  (as  he  agrees)  is  both  a 
powerful  and  limiting  assumption  for  his  technique.  Any  modulation 
in  the  wave  scene  other  than  tilt  modulation  would  decrease  the  es¬ 
timated  wave  height.  For  waves  not  purely  range-travel! ing,  the 
combined  effect  of  tilt  modulation,  straining  and  velocity  bunching 
is  poorly  understood  and  it  is  doubtful  whether  the  effect  of  tilt 
modulation  (and  thus,  wave  height)  could  be  sorted  out.  Kasischke 
(1980)  showed  that  wave  contrast  was  more  closely  correlated  to  wave 
height  than  to  wave  slope  for  a  variety  of  wave  conditions  and  radar 
look  angles.  Perhaps  one  does  not  have  to  assume  solely  tilt  modula¬ 
tion  and  use  Valenzuela's  curves  to  recover  wave  height  information 
from  image  contrast  measurements. 
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TABLE  10 

WAVE  HEIGHT  CALCULATIONS 

Data  -yi- 
Source  <it> 

o  . 

“  V  xo 

<  V  (dB)  tan  ®r  (m) 

Hs 

Sea  Truth  H, 

(m) 

(m) 

Orbit 
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Thomas  also  makes  the  assumption  that  any  variability  in  image 
intensity  over  the  areas  studied  was  due  to  either  ocean  waves  or 
speckle.  This  assumption  could  be  violated  due  to  local  wind  varia¬ 
tions,  extreme  wave  slopes  causing  specular  reflections  or  breakinq 
waves.  He  also  makes  the  assumption  that  the  speckle  characteristics 
of  the  scene  are  known;  this  is  probably  valid  but  should  be  shown 
experimental  ly. 

In  summary,  the  approach  of  Thomas  offers  a  possible  means  of 
extracting  ocean  wave  height  information  from  digitally  processed 
Seasat  SAR  imagery  for  the  case  of  range-traveling  waves.  Clearly, 
more  work  is  necessary  to  fully  evaluate  its  usefulness  and  limita¬ 
tions.  Further  tests  should  include:  (1)  varying  the  size  of  the 
area  where  statistics  are  generated  to  determine  the  sensitivity  of 
the  wave  contrast  (height)  measures  on  this  parameter;  (2)  applying 
the  algorithm  to  a  wider  range  of  wave  conditions  and  radar  look 
directions  to  determine  whether  wave  contrast  can  be  related  to  wave 
height  for  waves  not  purely  range-traveling;  and  (3)  studying  the 
speckle  characteristics  of  digitally  processed  Seasat  SAR  data  to 
determine  its  characteristics  experimentally. 


4 

DETECTION  OF  MESOSCALE  FEATURES 


RADAR  DIVISION 


Terim 


During  the  early  examination  of  Seasat  SAR  imagery,  it  became 
evident  that  a  variety  of  patterns  were  present  on  the  oceans'  sur¬ 
face  which  were  being  detected  by  the  Seasat  SAR.  Because  of  the 
variety  of  surface  and  subsurface  measurements  which  were  made  at 
the  times  of  the  Seasat  overpasses  during  the  JASIN  experiment,  this 
data  set  offers  an  opportunity  to  gain  a  better  understanding  of  why 
many  of  the  mesoscale  features  occurred  and  were  detected  by  the 
Seasat  SAR.  Because  many  of  the  surface  and  subsurface  measurements 
are  just  now  being  prepared  for  publication,  the  present  study  con¬ 
centrated  on  investigating  the  spatial  and  temporal  distribution  of 
the  mesoscale  patterns.  Linking  the  Seasat  SAR  imagery  with  the 
"sea-truth"  measurements  is  a  task  which  requires  further  efforts. 

4.1  THE  JASIN  TEST  AREA 

During  the  previous  gravity  wave  analysis,  we  were  primarily 
concerned  with  the  area  in  and  around  the  Oceanographic  Intensive 
Array  (see  Figure  1).  During  the  analysis  of  mesoscale  features,  we 
studied  Seasat  SAR  imagery  collected  over  the  Northeast  Atlantic 
Ocean,  in  an  area  bounded  by  the  United  Kingdom  and  the  Faeroe 
Islands  on  the  east,  18°  W  longitude  on  the  west,  54°  N  latitude  on 
the  south  and  64°  N  latitude  on  the  north.  This  area  is  presented 
in  Figure  9,  which  was  reproduced  from  the  Institute  of  Oceanographic 
Sciences  Chart  Number  C6567  (for  a  description  of  this  chart,  see 
Roberts,  et  al . ,  1979).  In  conjunction  with  the'  JASIN  experiment, 
extensive  hydrographic  surveys  of  this  portion  of  the  Atlantic  were 
made,  resulting  in  a  bathymetric  chart  which  is  unusually  accurate 
for  a  deep  water  area. 

There  are  17  distinct  bottom  topographic  features  present  in  this 
area.  These  features  are  listed  in  Table  11,  along  with  the  shallow¬ 
est  contour  interval  charted  for  that  feature,  the  estimated  deor.h 
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FIGURE  9.  LOCATION  OF  MAJOR  BOTTOM  TOPOGRAPHIC  FEATURES  IN  NORTHEAST 
ATLANTIC. 
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TABLE  11 

BOTTOM  TOPOGRAPHIC  FEATURES  IN  THE  NORTHEASTERN  ATLANTIC 


Bottom  Feature 

Shallowest  Point 
(m) 

Depth  of 
Adjacent  DW 
(m) 

Slope 

(Degree 

Iceland-Faeroe  Rise 

300 

1000 

5 

Bill  Bailey's  Bank 

100 

1200 

13 

Lousy  Bank 

300 

1500 

28 

Hatton  Bank 

500 

2000 

18 

George  Bligh  Bank 

500 

1100 

1 1 

Faeroe  Bank 

100 

1000 

8 

Faeroe  Shelf 

300 

1000 

7 

Wyvil 1  e-Thomson 

R  i  dge 

400 

1200 

17 

Ymir  Ridge 

600 

1500 

31 

North  Feni  Ridge 

1200 

2000 

5 

Rosemary  Bank 

500 

2000 

23 

Rockall  Bank 

300 

2300 

20 

Anton  Dohrn  Seamount 

600 

2100 

28 

West  Shetland  Shelf 

300 

1000 

11 

Malin  Shelf 

200 

2000 

31 

Hebrides  Terrace 
Seamount 

1000 

2400 

27 

Hebrides  Shelf 

200 

1500 

23 
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of  the  adjacent  deep  water  area,  and  the  calculated  slope  for  the 
feature.  As  will  be  discussed  later,  the  interaction  between  these 
bottom  features  and  ocean  currents  results  in  internal  waves  and 
frontal  boundaries  in  this  area. 

In  order  to  study  the  relationship  between  surface  patterns  and 
Seasat  SAR  imagery  and  these  bottom  features,  the  following  procedure 
was  used:  First,  the  location  of  the  ground  coverage  of  the  Seasat 
SAR  was  obtained  by  consulting  the  satellite  data  record  (SDR)  for 
each  pass.  (The  SDR  contains  a  series  of  ground  latitude  and  longi¬ 
tude  coordinates  for  specific  times  during  the  Seasat  pass.  Since 
most  of  Seasat  SAR  images  contained  annotations  of  the  time  the  data 
were  collected,  the  location  of  a  particular  surface  pattern  could 
be  identified  by  noting  the  time  when  the  pattern  occurred  during 
that  orbit.)  Next,  the  ground  coverage  of  each  pass  was  marked  on 
an  overlay  placed  on  the  chart.  Finally,  the  locations  of  all  inter¬ 
nal  waves,  frontal  boundaries  and  weather  patterns  were  noted  and 
then  marked  on  the  overlay. 

Imagery  from  15  of  the  18  JASIN  Seasat  passes  was  studied.  The 
imagery  from  three  of  the  passes  did  not  have  accurate  time  annota¬ 
tions  and  contained  no  land;  therefore,  they  could  not  be  accurately 
located.  Figure  10  shows  the  total  SAR  coverage  of  the  Northeast 
Atlantic  by  the  15  passes. 

4.2  OCCURRENCE  OF  MESOSCALE  FEATURES 

Four  distinct  mesoscale  features  were  present  on  the  JASIN/Seasat 
SAR  imagery:  internal  waves,  frontal  boundaries,  areas  of  low  radar 
backscatter  believed  to  be  weather  patterns,  and  windrows.  In  this 
section  we  will  first  present  examples  of  these  surface  patterns  and 
then  discuss  their  spatial  and  temporal  distribution  on  the  JASIN/ 
Seasat  imagery.  The  locations  of  the  five  examples  of  SAR  imaces 
presented  in  this  chapter  are  shown  in  Figure  11. 
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FIGURE  11.  LOCATIONS  OF  JASIN  TEST  AREA  AND  AREAS  EXAMINED  DURING  DEEP 
WATER  SURVEY.  (After  Institute  of  Oceanic  Sciences  Chart 
No.  C6567,  Depth  in  Meters) 
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The  JPL-optical ly  processed  imagery  presented  in  Figure  12  was 
collected  during  Rev.  599  (7  August  1978)  as  Seasat  passed  over  the 
Wyv  ill  e-Thomson  Ridge.  The  corresponding  bathymetric  chart  for  this 
coverage  is  presented  in  Figure  13.  From  Figure  12,  it  can  be  ob¬ 
served  that  there  are  numerous  internal  wave  signatures.  Note  on 

the  SAR  image  that  the  internal  waves  are  both  range  (e.g.,  C/Dl  to 

C/D 2)  and  azimuth  (e.g.,  A9  to  09)  traveling.  The  internal  waves  in 
the  lower  left  corner  of  the  SAR  image  appear  to  occur  over  the  Ymir 
Ridge;  the  internal  waves  in  the  middle  of  the  image  appear  over  the 
Wyville-Thomson  Ridge;  and  those  at  the  top  of  the  image  occur  over 
the  Faeroe  Bank  Channel . 

Surface  measurements  collected  at  the  time  of  the  Seasat  over¬ 
pass  indicate  a  wind  of  6.9  m/s  from  the  north  and  wave  field  with  a 
significant  wave  height  of  1.1  meters,  a  dominant  wavelength  of  71 
meters  propagating  towards  20*(T). 

Figure  14  was  collected  by  Seasat  during  Rev.  791  (21  August 

1978)  as  it  passed  over  the  Anton  Oohrn  Seamount.  The  corresponding 

bathymetry  for  this  imagery  is  presented  in  Figure  15.  The  internal 
wave  signatures  collected  over  the  Anton  Dohrn  Seamount  are  quite 
subtle.  Three  groups  of  internal  waves  occur  over  this  seamount  (B6 
to  C6,  D/E6  to  D/E8,  and  A8),  which  is  quite  steep  sided,  as  can  be 
seen  by  the  narrow  spacing  of  the  contour  lines  in  Figure  15. 

The  sea  truth  collected  at  the  JASIN  test  area  at  the  time  of 
Rev.  791  indicates  the  winds  were  quite  strong,  with  a  speed  of  13 
m/s  out  of  the  southwest.  The  surface  gravity  wave  field  had  a  sig¬ 
nificant  wave  height  of  3.3  m,  a  dominant  wavelength  of  132  m  and 
was  propagating  towards  230 "(T). 

The  JPL-optical ly  processed  image  presented  in  Figure  16  was 
collected  as  Seasat  passed  over  the  Faeroe  Bank  Channel  and  the  edge 
of  the  Faeroe  Shelf  during  Rev.  762.  The  corresponding  bathymetric 
chart  for  this  area  is  presented  in  Figure  17.  The  curveo,  dark  line 
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FIGURE  13.  GROUND  COVERAGE  CORRESPONDING  TO  SEASAT  SAR  REV.  599 
(After  IOS  Chart  No.  C6567,  Contour  Intervals  in 
100-Meter  Increments.) 
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FIGURE  15.  GROUND  COVERAGE  OF  REV.  791  OVER  ANTON  DOHRN  SEAMOUNT. 

(After  IOS  Chart  No.  C6567,  Depth  Contours  in  100-Meter 
Increments.) 


RADAR  LOOK  DIRECTION 


FIGURE  15.  S LAS AT  SAR  IMAGE  OF  A  FRONTAL  BOUNDARY  COLLECTED 
OVER  FAEROE  BANK  CHANNEL  AND  FAEROE  RIDGE,  REV, 
762.  (Image  courtesy  of  JPL.) 
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beginning  at  A4  and  continuing  to  B9  is  believed  to  be  a  surface 
pattern  caused  by  a  frontal  boundary.  Its  location  appears  to  be 
over  the  edge  of  the  Faeroe  Shelf,  an  area  with  a  water  depth  between 
600  and  1000  meters. 

Patterns  on  the  ocean  surface  caused  by  meteorological  phenomena 
are  quite  common  on  Seasat  SAR  imagery  (see  e.g.,  Beal,  et  al.,  1931; 
Kasischke,  et  al.,  1982;  or  Fu  and  Holt,  1982).  Those  observed  on 
the  JASIN-Seasat  imagery  appear  to  be  the  result  of  either  high  or 
low  winds. 

High  winds  can  result  in  a  meteorological  phenomena  known  as  a 
windrow  or  row  vortice  (LeMone,  1973).  Windrows  are  secondary, 
crosswind  circulations  which  have  a  periodic  structure  in  the  cross- 
wind  direction,  i.e.,  parallel  to  the  wind  direction.  Their  hori¬ 
zontal  wavelength  is  approximately  two  to  four  times  the  thickness 
of  the  planetary  boundary  layer.  This  periodic  variation  of  the 

surface  wind  results  in  a  corresponding  alteration  of  the  ocean  sur¬ 
face  roughness,  which  in  turn  can  be  detected  by  a  SAR.  Figure  18 
was  collected  during  Rev.  642.  Sea  truth  measurements  indicate  a 
wind  speed  of  12.9  m/sec  from  150 ’(T)  was  present  at  this  time.  As 
shown  in  Figure  18,  the  direction  of  the  wind  is  parallel  to  the 

direction  of  the  periodic  structures  in  the  SAR  image. 

Areas  of  low  velocity  surface  winds  in  ocean  areas  result  in  dark 

areas  on  SAR  imagery.  This  is  due  to  the  fact  that  in  these  areas, 
insufficient  winds  exist  to  generate  the  short-period  capillary  waves 
responsible  for  the  microwave  backscatter  from  the  ocean's  surface. 
Figure  19  presents  a  SAR  image  collected  during  Rev.  556.  Althouqh 
a  surface  wind  of  3.5  m/sec  was  present  in  the  JASIN  area,  it  is  sus¬ 
pected  that  much  lower  winds  were  present  in  the  area  of  Figure  19, 
which  is  located  aoproximately  100  kilometers  southeast  of  the  JAST'i 
OIA. 
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The  analysis  of  the  spatial  and  temporal  distribution  of  internal 
waves  and  frontal  boundaries  on  the  JASIN-Seasat  SAR  imagery  was 
supported  partially  by  funding  provided  by  the  Naval  Research 
Laboratory  (see  Kasischke,  et  a!.,  1982).  The  main  goal  of  this 
effort  was  to  investigate  the  relationship  between  these  deep-water 
surface  patterns  and  the  seventeen  distinct  bottom  features  in  the 
area  (see  Table  11).  Most  of  the  internal  waves  and  frontal  boundary 
patterns  occurred  quite  near  these  features. 

Table  12  summarizes  the  occurrence  of  internal  wave  and  frontal 
boundary  patterns  on  Seasat  SAR  imagery  for  each  of  the  seventeen 
deep  water  bottom  features.  Table  13  presents  these  occurrences  for 
each  Seasat  revolution.  Given  in  Table  13  are  the  Seasat  revolution 
number,  the  total  number  of  deep  water  bottom  features  within  the 
ground  coverage  of  a  pass,  and  the  total  number  of  times  an  internal 
wave  or  frontal  boundary  pattern  occurred  over  a  bottom  feature. 
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Presented  in  Figure  20  are  the  locations  of  the  internal  wave 
features  which  appeared  on  the  Seasat  imagery.  The  length  of  the 
lines  indicate  the  two-dimensional  extent  of  each  internal  wave 
packet.  The  lines  with  the  arrows  indicate  which  direction  the 
crests  of  the  internal  waves  are  traveling.  If  no  direction  could 

be  determined,  then  an  arrow  was  placed  at  both  ends.  Plotted  on 
Figure  21  are  the  locations  and  extent  of  the  frontal  boundaries 

identified  on  the  JASIN  SAR  imagery. 

From  Table  13,  we  can  see  that  63  percent  of  the  internal  waves 

and  38  percent  of  the  frontal  boundaries  were  associated  with  a  bot¬ 

tom  feature.  An  interesting  result  of  the  mul ti -temporal ,  deep  water 
analysis  was  that  at  no  time  did  an  internal  wave  pattern  occur  over 
a  deep  water  area  (i.e.,  the  Rockall  Trough,  Iceland  Basin  or 
Norwegian  Basin).  These  results  strongly  suggest  a  correlation  be¬ 
tween  the  deep  water  bottom  features  and  the  internal  waves  and  also 
indicate  some  relationship  may  exist  between  the  occurrence  of 
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TABLE  12 

SUMMARY  OF  OCCURRENCES  OF  BOTTOM-RELATED  SURFACE  PATTERNS 
ON  SEASAT  SAR  IMAGERY  (Continued) 


Bottom  Feature 
Covered 
by  Seasat  SAR 


SAR 

Revolution 

Number 


SAR  Observed  Pattern 

Frontal 

Internal  Wave  Boundary 


Faeroe  Shelf 


Wyv ill e-Thomson  Ridge 


Ymir  Ridge 


North  Feni  Ridge 


Rosemary  Bank 


556 

599  X 

642  X 

757 

762  X 

556  X 

599  X 

642 

719  X 

757 

762  X 

556  X 

599  X 

642  X 

719  X 

757  X 

1049  X 

547  X 

556 

633 

757 

791  X 

834 

1006  X 

1049 

556  X 

599  X 

633  X 

757 

834  X 

1006  X 

1049  X 


X 

X 

X 

X 


X 


X 


X 

X 


X 

X 

X 


X 


X 
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TABLE  12 

SUMMARY  OF  OCCURRENCES  OF  BOTTOM-RELATED  SURFACE  PATTERNS 
ON  SEASAT  SAR  IMAGERY  (Concluded) 


Bottom  Feature 

SAR 

SAR  Observed 

★ 

Pattern 

Covered 

Revolution 

Frontal 

by  Seasat  SAR 

Number 

Internal  Wave 

Boundary 

Rockall  Bank 

547 

X 

556 

X 

757 

X 

791 

X 

958 

X 

X 

1044 

1087 

X 

Anton  Dohrn  Seamount 

547 

X 

X 

556 

X 

599 

X 

791 

X 

West  Shetland  Shelf 

642 

X 

762 

X 

Malin  Shelf 

547 

X 

633 

X 

791 

X 

834 

Hebrides  Terrace 

547 

X 

Seamount 

642 

X 

791 

Hebrides  Shelf 

599 

633 

642 

719 

762 

X 

834 

1006 

X 

1049 

TOTAL 

92 

58 

35 

An  x  indicates  that  surface  pattern  appeared  over  the  position 
of  the  bottom  feature  during  a  particular  Seasat  revolution. 


RADAR  DIVISION 


v  ERIM 


TABLE  13 

SUMMARY  OF  DETECTIONS  OF  DEEP  WATER  TOPOGRAPHIC  FEATURES 
BY  SEASAT  SAR  IMAGERY  USING  DATA  COLLECTED 
OURING  THE  JASIN  STUDY 


Seasat 

Revolution 

Total 

Number  of 

Deep  Water  Bottom 
Features  Covered 

Total 

Number 

of  Occurrences 
of  Internal 

Waves  Over 
Features 

Total 
Number  of 
Occurrences 
of  Frontal 
Boundaries 
Over  Features 

547 

7 

5 

3 

556 

8 

4 

4 

599 

6 

5 

1 

633 

5 

3 

2 

642 

6 

4 

3 

719 

6 

4 

2 

757 

8 

3 

1 

762 

6 

6 

4 

791 

8 

4 

4 

834 

5 

2 

0 

958 

5 

4 

2 

1006 

5 

2 

4 

1044 

6 

4 

2 

1049 

6 

4 

2 

1087 

6 

4 

1 

Total 

(Percent) 

89 

58 

(63) 

35 

(38) 

68 
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frontal  boundaries  and  deep  water  features.  They  also  suggest  that 
the  dissipation  of  deep  water  propagating  internal  wave  pockets  is 
sufficiently  high  to  preclude  their  detection  through  successive 
passes. 

Table  14  surmiarizes  the  occurrence  of  windrows  and  areas  of  low 
radar  return  on  the  JASIN-Seasat  imagery.  The  patterns  were  not  as 
widely  spread  nor  as  frequent  as  the  frontal  boundaries  and  internal 
waves.  Previous  analyses  of  Seasat  SAR  imagery  have  indicated  that 
areas  of  low  radar  return  collected  over  water  usually  indicate  areas 
of  low  wind  (see  Beal,  et  al.,  1981).  Without  any  sea  truth  from 
the  exact  areas  in  Table  14  where  dark  regions  occurred,  it  is  dif¬ 
ficult  to  categorically  state  that  these  are  low-wind  areas  also, 
but  it  is  our  belief  that  they  are. 

Of  the  four  times  windrows  were  present,  on  two  occasions  the 
windrows  were  oriented  parallel  to  the  wind  direction,  on  one  occa¬ 
sion  they  were  oriented  40°  to  the  reported  wind  direction,  and  on 
one  occasion  no  sea-truth  was  available. 

4.3  ANALYSIS  OF  INTERNAL  WAVE  AND  FRONTAL  BOUNDARY  PATTERNS 

Because  of  their  widespread  and  frequent  appearance,  the  internal 
wave  and  frontal  boundary  patterns  visible  on  the  JASIN/Seasat 
imagery  were  further  examined. 

4.3.1  INTERNAL  WAVE  SIGNATURES 

The  surface  environmental  conditions  present  at  the  times  of  the 
Seasat  overpasses  were  compared  to  the  occurrence  of  the  internal 
waves  to  determine  if  the  absence  of  patterns  could  be  explained  by 
the  theory  that  adverse  surface  conditions  (i.e.,  high  winds  and/or 
waves)  would  mask  the  internal  wave  patterns.  The  surface  conditions 
are  presented  in  Table  2. 
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TABLE  14 

SUMMARY  OF  THE  OCCURRENCE  OF  WINDROWS  AND 
AREAS  OF  LOW  RADAR  RETURN 


Seasat 

Revolution 

Windrows 

Areas  of  Low 
Radar  Return 

547 

No 

Yes 

556 

No 

Yes 

599 

No 

No 

633 

No 

Yes 

642 

Yes 

No 

714 

Yes 

No 

719 

No 

No 

757 

No 

No 

762 

No 

No 

791 

No 

No 

834 

No 

No 

958 

No 

Yes 

1006 

Yes 

No 

1044 

No 

Yes 

1049 

No 

Yes 

1087 

Yes 

No 

71 
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The  Seasat  Revolutions  were  divided  into  two  groups  for  this 
analysis.  Group  I  contained  all  those  revolutions  where  the  internal 
wave  detection  rate  was  greater  than  70  percent.  Group  II  contained 
all  revolutions  where  the  internal  wave  detection  rate  was  less  than 
70  percent. 

The  average  wavelength,  wave  height,  wind  speed,  and  direction 
of  wave  travel  with  respect  to  radar  look  direction  for  each  group 
are  presented  in  Table  15.  Also  presented  in  Table  15  are  the  range 
of  these  values  for  each  group.  Table  15  shows  that  no  strong  pat¬ 
terns  exists  in  the  data,  with  the  internal  waves  being  detected  over 
a  wide  variety  of  environmental  conditions,  including  high  surface 
winds . 

It  has  been  suggested  that  a  preferred  look  direction  for  viewing 
internal  waves  with  SARs  exists  with  internal  waves  whose  crest 
orientations  are  perpendicular  to  the  radar  1 ine-of-sight  (i.e., 
range  traveling)  being  detected  more  often  by  SARs  than  those  whose 
crests  are  parallel  (i.e.,  azimuth  traveling)  (Fu  and  Holt,  1982). 
Although  cases  exist  where  range  traveling  internal  waves  have  been 
imaged  by  SARs  (see  Vesecky  and  Stewart,  1982  or  Figure  11,  this 
report),  these  are  thought  to  be  the  exception  rather  than  the  rule. 

The  orientations  with  respect  to  the  radar  1 ine-of-sight  of  the 
crests  of  the  internal  waves  on  the  JASIN  SAR  imagery  were  deter¬ 
mined.  Four  orientation  classifications  were  used: 


1.  0o-30°:  The  crests  of  the  internal  waves  are  parallel  to 

the  radar  1 ine-of-sight  (i.e.,  azimuth  traveling). 

2.  30°-60°:  The  crests  of  the  internal  waves  are  half  wav  be¬ 

tween  range  and  azimuth  traveling. 

3.  60°-90°  The  crests  of  the  internal  waves  are  perpendicular 

to  the  radar  line-of-sight  (i.e.,  range  travel  ing). 

4.  0°-90°:  The  crests  of  the  internal  waves  are  oriented  in 

every  direction . 


TABLE  15 

COMPARISON  OF  INTERNAL  WAVE  DETECTION  GROUPS  TO 
ENVIRONMENTAL  PARAMETERS 


Environmental  Condition 


Group  I 

(Internal  Wave 
Detection  >70%) 


Group  II 
(Internal  Wave 
Detection  <70-0 


Wavelength 

Average 

169 

m 

175 

m 

Mini mum 

71 

m 

105 

m 

Maximum 

299 

m 

256 

m 

Wave  Height 

Average 

2.7 

m 

3.4 

m 

Minimum 

1.1 

m 

1.5 

m 

Maximum 

5.7 

m 

5.0 

m 

Wind  Speed 

Average 

8.9 

m/s 

10.3 

m/s 

Minimum 

4.5 

m/s 

3.6 

ml  s 

Maximum 

13.2 

m/s 

15.2 

m/s 

Wave  Direction  with 

Respect  to  Radar 

Look  Direction 

Average 

60 

Minimum 

25 

Maximum 

84 

47° 

14° 

85° 
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Table  16  summarizes  the  distribution  of  the  internal  wave  orien¬ 
tations  for  both  the  ascending  passes  and  the  descending  passes. 
Also  presented  in  Table  16  are  the  direction  of  orientation  for  the 
internal  waves.  A  direction  of  NE-SW,  for  example,  means  the  crest 
of  the  internal  wave  is  oriented  from  the  northeast  to  the  southwest. 
Note  from  the  data  in  Table  16,  that  44  internal  wave  packets  were 
observed  during  the  eight  ascending  orbits  and  44  internal  wave 
packets  were  observed  during  the  seven  descending  orbits. 

If  the  orientation  of  the  internal  waves  in  the  JAS IN  area  were 
equally  distributed  in  all  directions,  and  the  SAR  had  no  bias  in 
the  manner  it  detects  internal  waves,  then  we  would  expect  that  an 
equal  number  of  azimuth  (0o-30°)  waves  and  range  waves  (60°-90°) 
would  be  observed.  From  Table  16,  we  can  see  that  well  over  50  per¬ 
cent  of  the  detected  internal  waves  had  a  range  orientation.  How¬ 
ever,  we  can  also  see  that  the  majority  of  the  waves  were  oriented 
in  a  northeast  to  southwest  or  southeast  to  northwest  direction, 
indicating  the  waves  were  not  traveling  an  equal  amount  in  all  direc¬ 
tions.  It  therefore  appears  more  analysis  is  needed  to  interoret 
the  meaning  of  the  data  presented  in  Table  16. 

Throughout  the  analysis  of  the  JAS I N  data,  we  have  referred  to  a 
portion  of  the  anomalous  deep  water  surface  patterns  on  the  SAR 
images  as  internal  waves.  Although  we  have  yet  to  develop  direct 
evidence  that  these  surface  patterns  are  internal  waves,  the  circum¬ 
stantial  evidence  for  this  being  the  case  is  strong. 

The  generation  mechanisms  and  sources  responsible  for  internal 
waves  have  been  extensively  researched  during  recent  vears  ('.-lunch, 
1976;  Muller  and  Olbers,  1975;  Bell,  1975;  Baines,  1979).  Even  so, 
the  knowledge  of  the  physical  processes  for  internal  wave  generation 
is  still  very  limited  because  a  determination  of  the  sources  ana 
sinks  of  internal  waves  requires  a  continuous  record  of  the  .vave 
spectra  (Garrett  and  Munk,  1975  and  1979). 
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TABLE  16 

ORIENTATION  OF  INTERNAL  WAVES  OBSERVED  ON  SEASAT  SAR  IMAGERY 
COLLECTED  OVER  THE  JASIN  AREA 


★ 


Internal  Wave  Orientation 


Seasat  Revolution 

0 '-30 " 

J.  1  1  O  1  liu  1  H  U*  V, 

"  30'-60” 

i  i  i  i  o  u  i*  i  ■ 

0  °-90  ‘ 

Descending  Orbits 

E-W 

(Direction  of  Orientation) 

NE-SW 

N-S 

SE-NW 

All 

556 

1 

1 

3 

0 

599 

4 

3 

2 

1 

642 

5 

3 

2 

1 

757 

0 

1 

5 

0 

958 

1 

0 

3 

0 

1044 

2 

0 

2 

0 

1087 

0 

4 

0 

0 

Subtotal 

13 

12 

17 

2 

Ascending  Orbits 

E-W 

(Direction  of  Orientation) 

SE-NW 

N-S 

NE-SW 

All 

547 

0 

3 

4 

1 

633 

1 

1 

5 

1 

719 

2 

0 

4 

0 

762 

0 

1 

S 

0 

791 

0 

1 

4 

o 

834 

0 

0 

2 

0 

1006 

0 

0 

1 

0 

1049 

3 

0 

2 

0 

Subtotal 

6 

6 

30 

o 

Total 

19 

18 

47 

4 

N-S 

E-W  NE-SW 

SE-NW 

A 1  1 

Total  by 

Direction  of  Orientation 

4 

14  43 

9  9 

*A  0°  orientation  indicates  the  crests  of  the  internal  vaves 
are  aligned  parallel  to  the  radar  1  ine-of-siqht ;  a  °0  °  orientat ' 
indicates  they  are  aligned  perpendicular  to  the  raaar  1  ■» v 
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Several  sources  for  the  generation  of  internal  waves  have  been 
identified,  including  surface  atmospheric  effects,  internal  (oceanic) 
effects  and  topographic  effects  (Thorpe,  1975).  Surface  atmospheric 
effects  include  traveling  atmospheric  pressure  and  stress,  buoyancy 
flux,  surface  wave  interactions  and  Ekman  layer  instability.  Inter¬ 
nal  (oceanic)  generation  can  be  caused  by  decay  of  large  scale  cir¬ 
culations  and  breaking  baroclinic  instability  (of  which  there  is  very 
little  known).  The  generation  of  internal  waves  due  to  topographic 
features  is  similar  to  the  generation  of  Lee  waves  by  air  flow  over 
irregular  surfaces.  In  the  ocean,  internal  waves  can  be  generated 
when  a  current  (tidal  or  otherwise)  flows  over  a  deep  water  feature 
such  as  a  seamount,  guyot,  or  submarine  ridge  on  the  edge  of  a  con¬ 
tinental  shelf  (Wunch,  1975;  Bell,  1975;  Baines,  1979). 

A  recent  analysis  of  the  JASIN  imagery  gave  further  evidence  for 
the  occurrence  of  internal  wave  signatures  on  the  JASIN  Seasat  SAR 
imagery  [Liu  (1983)  and  Kasischke,  et  al .  (1983)].  The  following 
discussion  is  from  Liu  (1983). 

Figure  22  shows  the  Seasat  SAR  imagery  of  deep  water  features 
over  the  Wyvil le-Thomson  Ridge  collected  during  Revolution  762  on  19 
August  1978  at  approximately  06:45  GMT.  Figure  23  combines  line 
drawings  of  internal  waves  and  a  frontal  boundary  (discussed  below) 
observed  in  Figure  22  with  the  corresponding  bathymetric  chart  for 
this  area.  From  Figure  23,  it  can  be  seen  that  there  are  numerous 
internal  wave  signatures.  The  internal  waves  in  the  middle  of  the 
image  appear  to  be  generated  over  the  Wyvil  le-Thomson  Ridge.  Near 
the  lower  right  corner  of  the  image,  three  wave  packets  are  visible. 
It  is  to  be  noted  ‘hat  no  waves  were  generated  on  the  north  side  of 
the  r i dge . 

The  dominant  bottom  features  are  a  broad,  saddle-like  vallev  on 
the  East  side  of  the  ridge,  and  a  90°  sector  of  deeo  water  extendin'! 
out  'rom  the  ridge,  as  shown  in  Figure  24.  The  ridae  has  the 
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SEASAT  SAR  IMAGERY  OF  DEEP  WATER  FEATURES  OVER  THE  WYVILLE 
THOMSON  RIDGE.  REVOLUTION  762  (After  Liu,  1983) 


FIGURE  23.  LINE  DRAWING  OF  INTERNAL  WAVES  AND  A  FRONTAL  BOUNDARY 
OBSERVED  IN  FIGURE  22  SUPERIMPOSED  ON  THE  HYDROGRAPHIC 
CHART  OF  THE  NORTHEAST  ATLANTIC  OCEAN  (After  Liu,  1983) 


WE  I  LAND  ^ 
>ilELF 


W Y V I L  LE  -THOMSON 
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characteristic  scales  of  approximately  10  km  in  both  the  longitudinal 
and  transverse  direction.  The  minimax  of  the  bathymetric  saddle  is 
at  a  depth  of  about  450  m. 

Based  on  the  bathymetry,  it  seems  that  the  wave  packets  are  pro¬ 
duced  by  lee  wave  formation — a  mechanism  similar  to  the  undulation 
of  the  jet  stream  as  it  passes  over  a  mountain  ridge.  Figure  25 
schematically  illustrates  the  stages  that  are  hypothesized  during 
the  production  of  solitary  internal  wave  packets.  As  the  strong 
semidiurnal  tidal  current  flows  north  into  the  Faeroe  Bank  Channel, 
an  internal  lee  wave  forms  on  the  northern  edge  of  the  ridge.  As 
the  tidal  flow  goes  to  zero  six  hours  later,  the  lee  wave,  trying  to 
maintain  its  group  velocity  relative  to  the  current,  escapes  south 
over  the  ridge  barrier  back  into  the  deep  water,  emerging  as  a  Dacket 
of  solitary  waves.  This  phenomenon  has  previously  been  studied  by 
Holbrook,  et  al .  (1983)  in  the  Sulu  Sea.  The  initial  waveform  then 
can  be  viewed  as  a  localized  source  that  immediately  begins  to  prop¬ 
agate,  and  evolves  into  a  series  of  solitary  waves  while  undergoing 
radial  spreading,  encountering  variable  topography  and  slowly  losing 
energy  through  dissipation.  The  disappearance  of  internal  waves  on 
the  north  side  of  the  ridge  is  probably  because  there  is  no  strong 
tidal  current  flowing  south  across  the  ridge;  most  of  the  tidal  cur¬ 
rent  may  flow  along  the  Faeroe  Bank  Channel. 

Based  on  the  Seasat  SAR  image,  wave  packets  in  deep  water  over 
the  Wyvil le-Thomson  Ridge  have  the  following  properties: 

1.  Wave  packets  are  separated  by  distances  of  the  order  of  14 
km. 

2.  Wave  speed  is  estimated  to  be  0.31  m/s  based  on  the  assumo- 
tion  of  a  semidiurnal  tidal  origin. 

3.  Wave  packets  in  this  region  contain  about  five  waves  with 
packet  width  about  5  km. 
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4.  Wavelengths  are  between  500  m  and  1000  m. 

5.  The  crests  are  curved  in  a  horizontal  plane  with  their  convex 
sides  pointed  in  the  direction  of  propagation. 

Although  no  in  situ  measurements  are  available,  it  is  believed 
this  area  is  dominated  by  tidal  currents  with  velocities  of  less  than 
0.5  m/s  at  the  surface.  Figure  26  shows  the  profiles  of  potential 
temperature,  salinity  and  sigma  theta  generated  from  data  collected 
over  the  Icel and-Faeroe  Rise  during  the  Atlantic  Expedition  in  August 
1972  (Bainbridge,  1980).  Station  19  is  expected  to  be  the  measure¬ 
ment  closest  to  the  Wyv ill  e-Thomson  Ridge  area.  As  evident  from 
Figure  26,  a  mixed  layer  of  40  m  is  dominant  in  this  area  during  the 
summer.  Based  on  the  two-layer  finite  depth  model,  the  linear  wave 
speed  is  given  by 

[c  =  AS.  a 
L  0  0  k 

*3 

Near  the  Wyvil  1  e-Thomson  Ridge,  ao/d  -  0.6  x  10  ,  =  40  m, 

H2  =  260  m,  2Wk  -  500  m,  and  thus  the  wave  speed  is  approximately 
0.40  m/s,  depending  on  the  water  depths.  This  value  is  consistent 
with  the  estimated  wave  speed  of  0.31  m/s  from  the  SAR  image.  Based 
on  the  wavelength,  the  wave  amplitude  is  estimated  to  be  between  10 
and  15  m.  Note  that  the  effect  of  shear  is  important  in  this  case, 
since  the  wave  speed  is  slow  compared  with  the  background  current. 
Also,  the  bottom  topography  near  the  West  Shetland  Shelf  is  complex, 
which  influences  the  evolution  of  wave  packets,  as  shown  in  Figure 
23.  A  similar  solitary  wave  theory  developed  by  Liu,  et  al .  (1983) 
can  be  applied  here.  However,  due  to  the  limited  scope  of  this  in¬ 
vestigation,  the  detailed  analysis  and  calculations  will  have  to 
wait  until  sea-truth  data  are  available. 

Liu  (1983)  also  modeled  the  surface  strain  rate  produced  by  shal¬ 
low  water  (coastal)  internal  waves.  He  then  used  this  surface  strain 


i  i1/2 

coth  kH^  +  coth  kH^ 
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as  an  input  to  an  electromagnetic  model  and  simulated  a  SAR  image  of 
the  internal  wave,  which  compared  favorably  to  the  Seasat  image.  A 
similar  surface  strain  is  expected  to  occur  for  the  deep-water  inter¬ 
nal  waves,  resulting  in  the  patterns  observed  on  the  JASIN  imagery. 

4.3.2  FRONTAL  BOUNDARIES 

Even  though  they  were  not  as  prevalent  on  the  JASIN  Seasat  SAR 
imagery,  the  case  for  frontal  boundaries  also  being  related  to  deep 
water  bottom  features  is  strong. 

A  front  has  been  described  as  a  convergence  zone  between  Lwo 
ocean  water  masses  or  currents  of  different  velocities.  Hence,  a 
frontal  boundary  detected  on  SAR  imagery  is  usually  characterized  bv 
an  altered  surface  roughness  caused  by  the  formation  of  white  caps, 
the  bunching  of  the  surface  capillary  waves  due  to  water  velocity 
differences,  or  the  alteration  of  capillary  waves  by  surface  slicks 
or  foam  fronts. 

There  are  four  different  types  of  fronts  to  consider:  (1)  major 
ocean  currents;  (2)  mixing  of  water  from  estuaries  and  rivers; 
(3)  coastal  upwelling;  and  (4)  deep  water  upwelling.  Fronts  at  the 
boundary  of  major  ocean  currents  can  extend  over  thousands  of  kilom¬ 
eters.  The  fronts  are  associated  with  the  intrusion  of  either  warm, 
salty  water  of  tropical  origin  into  higher  latitudes  ( e . g . ,  the  Gulf 
Stream)  or  the  intrusion  of  cold,  less  saline  water  into  lower  lati¬ 
tudes  (e.g.,  the  Labrador  Current).  These  currents  are  generally  in 
quasi-geostrophic  balance. 

Two  smaller  scale  fronts  which  are  generally  restricted  to  con¬ 
tinental  or  near-continental  regions  are  those  found  at  river  or 
estuary  mouths  and  coastal  upwellings,  which  are  the  result  of  an 
offshore  surface  Ekman  transport  due  to  alongshore  wind  stresses. 
These  fronts  are  formed  in  continental  seas  and  estuaries,  and  are 
found  around  islands,  banks,  capes,  shoals,  and  generally  in  ’■he 
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boundary  regions  between  shallow,  wind  and  tidally  mixed  near-shore 
waters  and  stratified,  deeper  offshore  waters.  These  fronts  are  more 
influenced  by  local  acceleration,  bottom  stress  and  interfacial  fric¬ 
tion  than  by  Coriolis  force. 

The  JASIN  area  is  sufficiently  far  away  from  the  British  Isles 
to  rule  out  the  later  two  small  scale  fronts  as  those  observed.  Sur¬ 
face  measurements  of  currents  made  during  the  JASIN  experiment 
(Tarbell,  et  al.,  1979)  and  during  previous  years  (Bainbridge,  1980) 
indicate  that  no  persistent  large  scale  surface  currents  exist  in 
this  region.  It  is  believed  this  area  is  dominated  by  currents  with 
velocities  of  less  than  0.5  m/s  at  the  surface  (U.S.  Naval  Oceano¬ 
graphic  Office,  1977).  These  tidal  currents  can  also  generate  areas 
of  upwelling.  The  following  discussion  of  a  front  observed  in  the 
JASIN  imagery  again  follows  from  Liu  (1983). 

The  cold  front  that  forms  the  boundary  between  the  two  layer 
system  on  the  off-shore  side  and  the  vertically-mixed  water  in-shore 
has  been  observed  on  Seasat  SAR  images  of  the  Faeroe  Bank  Channel  . 
The  curved,  dark  line  beginning  at  48  and  continuing  to  55,  as  indi¬ 
cated  in  Figure  23,  is  believed  to  be  a  surface  pattern  caused  by  a 
frontal  boundary.  Its  location  appears  to  be  over  the  edge  of  the 
Faeroe  Shelf,  and  closely  matches  the  800  meter  bottom  contour. 

As  shown  in  Figure  26,  the  deep  water  upwelling  results  in  a 
frontal  boundary  over  the  Icel andic-Faeroe  Rise.  The  relationship 
of  the  front  to  local  upwelling  topographic  mixing  may  be  described 
by  the  following  process.  Figure  27  shows  the  cross-section  diagram 
of  bottom  topography  of  the  Wyv il le-Thomson  Ridge.  It  appears  that 
the  tidal  current  flowing  over  the  ridge  through  another  saddle-like 
valley  on  the  west-side  of  the  ridge  forces  the  deep  water  isotherms 
up  toward  the  surface.  To  illustrate  this  mechanism  for  frontal 
boundary  generation,  a  schematic  diagram  of  deep  water  upwelling  over 
the  Faeroe  Bank  Channel  is  shown  in  Figure  28.  Therefore,  the 
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This  report  further  analyzed  Seasat  SAR  imagery  collected  during 
the  J AS  I N  experiment.  Two  areas  were  investigated:  SAR  detection 
of  gravity  waves  and  SAR  detection  of  large-scale  features.  Fast 
Fourier  transforms  (FFTs)  from  all  18  JASIN/Seasat  orbits  were  gener¬ 
ated.  Estimates  of  dominant  wavelength  and  direction  of  propagation 
were  obtained  from  six  of  the  FFTs.  The  rest  did  not  produce  spectra 
which  could  be  interpreted.  Optical  Fourier  transforms  (OFTs)  gener¬ 
ated  interpretable  spectra  with  greater  success  than  Fast  Fourier 
Transformers  (FFTs).  This  is  most  likely  due  to  the  larger  aperture 
size  used  to  generate  the  OFTs. 

Comparison  of  the  wave  contrast  measurement  (peak-to-back ground 
ratio  or  PBR)  with  surface  truth  measurements  revealed  a  significant 
statistical  relationship  between  wave  contrast  and  dominant  wave¬ 
length  (x)  and  wave  height  (H ^3). 

SAR  data  processing  algorithms  which  compensate  for  the  motion 
of  surface  gravity  waves  were  evaluated  using  Seasat  SAR  imaqerv  from 
Revs.  762  and  1087.  Both  teleshift  (azimuth)  and  telerotation 
(range)  adjustments  were  tested.  It  was  found  that  Seasat  SAR 
imagery  of  gravity  waves  was  not  affected  using  these  adjustments. 
This  was  true  in  terms  of  both  wave  contrast  (i.e.,  wave  visibility! 
as  well  as  the  estimates  of  dominant  wavelenqth  and  direction  of 
propagation . 

An  algorithm  by  Thomas  (1982)  which  estimates  wave  heioht  *>o;T1 
digital  SAR  data  was  evaluated.  Althouqh  the  algorithm  showed  prom¬ 
ise,  it  raised  more  questions  than  it  answered.  Further  evaluations 
of  this  technique  should  include  varying  the  size  of  the  area  us-3: 
to  generate  the  wave  height  measure  and  aDDlyina  the  alcorith" 
wi der  range  of  wave,  wind  and  radar  look  directions. 
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Mesoscale  patterns  studied  on  SAR  imagery  collected  during  the 
OASIN  experiment  included  weather  patterns,  windrows,  frontal  boun¬ 
daries  and  internal  waves.  Large  areas  of  low  radar  return  present 
on  Seasat  imagery  are  suggested  to  be  areas  of  low  wind.  The  loca¬ 
tion  of  these  areas  were  compared  to  the  locations  of  major  weather 
fronts,  but  no  correlation  was  found.  Windrow  patterns  were  detected 
on  Seasat/JASIN  Imagery  and  were  found  to  be  aligned  parallel  to  the 
surface  truth  wind  direction. 

The  locations  of  all  frontal  boundary  and  internal  wave  patterns 
on  JASIN/Seasat  Images  were  plotted  on  a  bathymetric  chart  of  the 
Northeast  Atlantic  Ocean.  These  plots  revealed  that  these  patterns 
occurred  over  or  close  to  a  deep  water  bottom  feature,  such  as  a 
ridge,  seamount,  bank  or  shelf.  At  no  time  did  a  pattern  occur  over 
a  deep  water  basin.  An  Internal  wave  pattern  appeared  in  Seasat  SAR 
imagery  63  percent  of  the  time  Seasat  passed  over  a  deep  water  fea¬ 
ture,  while  frontal  boundary  patterns  appeared  38  percent  of  the 
time.  The  causes  of  internal  waves  and  frontal  boundaries  on  the 
Seasat  SAR  Imagery  have  been  identified  as  Lee  wave  generation  and 
deep  water  upwelling,  respectively.  Research  which  combines  oceanic, 
subsurface  measures  with  the  Seasat  SAR  imagery  is  definitely  war¬ 
ranted  to  further  study  these  patterns. 
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APPENDIX  A 
SEASAT/JASIN  PASSES 


Presented  In  this  appendix  Is  a  set  of  figures  to  assist  In  the 
visualization  of  the  Seasat  SAR  ground  track  with  respect  to  the 
oceanographic  Intensive  array  of  the  JASIN  experiment.  Each  JASIN 
Seasat  pass  Is  Illustrated  In  a  separate  figure.  Each  figure  contains 
the  location  of  the  Seasat  ground  track  superimposed  on  a  map  of  the 
JASIN  area  as  well  as  a  schematic  diagram  of  the  satellite  direction, 
radar  look  direction,  the  surface  wave  and  surface  wind  fields. 
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APPENDIX  B 

DIGITAL  FOURIER  TRANSFORMS  OF  JASIN 
SEASAT  SAR  DATA 

This  appendix  contains  copies  of  the  two-dimensional  FFTs  and  their 
resultant  one-dimensional  plots  for  these  passes  where  Interpretable 
spectra  were  obtained. 
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APPENDIX  C 

"THE  ESTIMATION  OF  HAVE  HEIGHT  FROM  DIGITALLY 
PROCESSED  SAR  IMAGERY"* 


*  M.H.B.  Thomas,  Int,  J,  Remote  Sensing.  3,  pp.  63-68,  1982. 
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The  estimation  of  wave  height  from  digitally  processed  SAR  imagery 
M.  H.  B  THOMASt 

Materials  Physics  Division.  AERE  Harwell.  Didcot.  OXI I  ORA,  England 


(Received  29  June  19HI .  revision  received  3  Sovemher  I9S1) 

Abstract.  A  simple  method  has  been  developed  for  estimating  wave  height  from 
synthetic  aperture  radar  (SAR)  imagery  obtained  by  the  Seasat  satellite.  The  method 
is  based  on  measuring  the  contrast  of  the  image  and  the  wavelength  of  the  dominant 
wave.  A  calculation  has  been  made  for  two  orbits  made  by  the  satellite  over  the 
North  Atlantic  in  1978.  using  digitally  processed  data  supplied  by  DFVLR  in  West 
Germany.  Comparison  with  sea  truth  measurements  shows  agreement  to  within 
about  20  per  cent. 

1.  Introduction 

As  part  of  a  programme  concerned  with  the  analysis  of  radar  imagery  of  the  sea. 
it  was  felt  worthwhile  to  investigate  the  possibility  of  obtaining  wave  height  from 
SAR  imagery.  Previously,  it  has  been  necessary  to  determine  this  using  a  radar  alti¬ 
meter.  For  satellite  applications,  this  has  an  obvious  disadvantage,  because  the  alti¬ 
meter  views  the  surface  directly  under  the  spacecraft  while  the  SAR  looks  out  to  the 
side.  Thus,  the  wave  height  is  being  measured  in  a  totally  different  area  from  that 
where  the  SAR  is  acquiring  its  imagery.  Clearly,  a  method  of  deriving  wave  height 
from  the  SAR  image  will  be  of  significant  practical  value. 

2.  Theory 

The  fundamental  interaction  between  radar  and  the  sea  surface  is  Bragg  scattering. 
Three  possible  mechanisms  have  been  proposed  (Alpers  et  al.  1981)  for  the  imaging 
of  waves  and  these  are : 

(a)  The  lifting  and  tilting  of  the  Bragg  scattering  waves  by  the  long-wavelength 
gravity  waves,  which  create  a  modulation  in  the  image  (tilt  modulation). 

( b )  The  asymmetry  of  capillary  waves  around  a  wave  crest  caused  by  non-linear 
wave  interactions  (straining). 

( c )  Doppler  shifts  associated  with  the  wave  motion  which  introduce  a  modulation 
into  the  SAR  imagery  (sometimes  known  as  velocity  bunching). 

The  expected  dependence  of  these  mechanisms  on  wave  direction  is  indicated  in 
figure  1  (Alpers  et  al.  1981).  Thus  for  waves  propagating  in  a  direction  which  is 
approximately  perpendicular  to  the  flight  direction,  mechanism  (c)  will  have  a  small 
effect.  W right  et  al.  (1980)  suggest  that  the  straining  mechanism  can  be  important  on 
the  basis  of  experiments  using  vertical  polarization.  However,  for  horizontal  polariza¬ 
tion  the  tilt  modulation  effect  is  considerably  greater  and  is  dominant  for  low  to 
moderate  sea  states  (Alpers  et  al.  1981).  For  the  present,  we  shall  assume  that  the 
dominant  imaging  mechanism  is  (a). 

If  we  denote  the  surface  displacement  from  the  mean  level  by  /(r.  r).  then  the 
’Present  address:  Logica  Ltd.  64  Newman  Street.  London.  W1A  4SE. 
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Figure  1.  Dependence  of  wave  imaging  mechanisms  on  wave  direction  (after  Alpers  el  al.  1981 ). 
significant  wave  height.  Ht.  is  defined  by 

H,=4n/</2>.  (1) 

For  a  narrow-band  swell-wave  spectrum  centred  on  wave  number  k0.  the  r.m.s.  slope 
of  the  sea  surface,  tan  0r,  is  given  by 

tan  0r _\k0\  it 

H,  4  2k0  ' 

By  rewriting  the  equation,  we  can  determine  Hs  if  we  know  tan  0,  and  x0.  From  a 
digital  image  it  is  straightforward  to  determine  k0  by  taking  the  Fourier  transform. 
The  determination  of  tan  8,  is  based  on  the  variation  of  the  backscattering  cross- 
section  of  the  sea  surface  with  incidence  angle. 

Figure  2  shows  the  theoretical  variation  of  backscatter  cross-section  for  horizontal 
polarization  (which  was  used  by  Seasat),  given  by  Valenzuela  (1978)  at  a  frequency  of 
4455  MHz  using  a  composite  surface  model ;  a  very  similar  curve  is  given  for  428  MHz. 
The  Seasat  frequency  of  1275  MHz  lies  between  these  values. 

The  tilt  modulation  mechanism  depends  on  the  change  in  the  local  angle  of  inci¬ 
dence  caused  by  the  swell  waves,  which  gives  rise  to  a  change  in  the  backscatter  cross- 
section.  This  intensity  modulation  then  allows  the  waves  to  appear  in  the  imagery. 
A  point  which  should  be  noted  here  is  that  the  variation  in  the  local  angle  of  incidence 
causes  a  variation  in  the  Bragg  scattering  wavelength.  The  backscatter  cross-section  is 
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Figure  2.  Variation  of  backscatter  cross-section  as  a  function  of  incidence  angle  for  horizontal 
polarization  at  4455  MHz  {Valenzuela  1978). 

proportional  to  the  power  spectrum  at  this  wavelength  and  this  has  been  included  by 
Valenzuela  by  assuming  a  k~*  dependence  (Phillips  spectrum). 

From  the  curve  of  figure  2,  it  is  possible  to  relate  a  given  variation  of  backscatter 
to  changes  in  the  local  angle  of  incidence  at  the  sea  surface.  The  immediate  problem 
then  is  to  determine  the  r.m.s.  variation  in  the  backscattered  power. 

This  can  be  determined  directly  from  the  image  because  the  pixel  intensity  in  the 
digital  data  is  proportional  to  the  square  root  of  the  backscattered  power.  In  fact,  in 
reading  data  from  tape  to  our  computer  disc,  the  pixel  intensities  have  been  squared 
so  that  they  are  now  directly  proportional  to  power.  The  standard  deviation  of  the 
backscattered  power  is  then  readily  obtainable.  However,  the  problem  is  complicated 
by  the  presence  of  speckle  in  the  SAR  image.  Speckle  is  characteristic  of  coherent 
imaging  systems  and  has  been  widely  studied  (Goodman  1975).  The  intensity  distribu¬ 
tion  of  a  speckle  pattern  is  a  negative  exponential  distribution  and  the  contrast  of  the 
pattern  can  be  defined  as  the  standard  deviation  of  the  intensity  divided  by  the  mean 
intensity.  For  the  negative  exponential  distribution,  this  is  equal  to  unity.  The  digital 
data  have  been  subjected  to  an  averaging  process  in  which  four  different  images  have 
been  incoherently  averaged,  specifically  to  reduce  the  speckle.  The  reduction  factor 
is  equal  to  the  square  root  of  the  number  of  separate  images,  so  that  the  contrast  of  the 
speckle  will  be  reduced  to  0-5. 

When  waves  are  present,  the  contrast  in  the  image  will  be  increased,  so  that  it 
should  be  possible  to  estimate  the  variation  in  backscatter  cross-section  from  the 
enhanced  contrast  of  the  image,  provided  that  the  speckle  contrast  can  be  accounted 
for.  The  local  average  intensity  of  the  speckle  is  proportional  to  the  backscatter  cross- 
section  at  that  point  in  the  image,  so  that  the  two  effects  combine  in  a  multiplicative 
manner  (Lee  1981). 

We  may  write 

where  /,  =  image  intensity,  ft  =  backscatter  cross-section.  /,  =  speckle  intensity,  and 
x  =  a  proportionality  constant.  Therefore 

</,2>  <n2xif> 

oy  w2  <a*>2</,>2 

if  ft  and  I,  are  independent. 
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We  may  use  the  relationship : 


where  x  =  a  measured  variable,  and  ax = the  standard  deviation  in  x.  It  now  follows 
that: 


(3) 


where  a,=  standard  deviation  of  the  image  intensity.  </,)=mean  image  intensity. 
0-,= standard  deviation  of  the  speckle  in  the  absence  of  any  waves,  </,>=mean 
speckle  intensity  in  the  absence  of  any  waves,  o,  =  standard  deviation  of  image  inten¬ 
sity  due  to  waves  in  the  absence  of  speckle,  and  </„>  =  mean  image  intensity  due  to 
waves  in  the  absence  of  speckle. 

Now  we  know  that  is  equal  to  0-5  for  a  four-look  image  and  we  can 

determine  <rf/</,>  directly  from  the  image ;  so  we  can  write 


+  1  =  1-25 


(4) 


from  which  we  can  calculate  <t„/<7w>.  From  here  it  is  straightforward  to  obtain  the 
value  of  the  r.m.s.  slope  using  Valenzuela’s  data  and  then  to  determine  the  significant 
wave  height 

We  form  the  quantity  V  =  10  log10[«/H,>  +  erJ/«/w>  —  a  J],  which  expresses  the 
backscatter  variation  in  dB  (Valenzuela's  curve  is  logarithmic).  This  allows  the 
determination  of  tan  6,  by  reference  to  the  curve  near  an  angle  of  22-6c,  which  is  the 
mean  angle  of  incidence  at  the  sea  surface  for  Seasat. 


3.  Results  and  discussion 

Two  Seasat  passes  have  been  investigated ;  these  were  pass  762  in  an  area  north¬ 
west  of  Scotland  and  pass  791  over  the  JASIN  area  in  the  Atlantic.  In  both  cf  these 
passes,  swell  waves  were  observed  travelling  with  an  azimuth  angle  <f>  (see  figure  1 ) 
approximately  equal  to  84-5°,  i.e.  almost  perpendicular  to  the  flight  direction  (see 
figure  3).  so  that  the  Doppler  mechanism  should  be  small.  The  contrast  of  the  image. 
at /</*>>  was  determined  and  then  aw/</w>  was  calculated  from  equation  (4).  This 
leads  directly  to  the  value  of  tan  6r.  The  dominant  wavelength  /.0  was  obtained  from  a 
Fourier  transform  and  the  significant  wave  height  then  calculated  from  equation  (2). 
Table  1  summarizes  these  results.  The  figures  for  tan  8r  are  averages  of  the  values 
derived  from  Valenzuela’s  curves  for  428  MHz  and  4455  MHz.  The  value  of  /.„  was 
obtained  by  calculating  the  centre  of  gravity  of  the  swell  part  of  the  power  spectrum. 
For  pass  762,  the  sea  truth  measurement  was  taken  from  a  shipbome  wave  recorder 
on  board  the  John  Murray  (Guymer  1980).  For  pass  791.  wave  rider  measure¬ 
ments  are  available  from  two  wave  riders  (Carlson  1980).  The  most  significant  errors 
in  the  calculation  are  associated  with  the  determination  of  A0  and  of  tan  6,  using 


Table  1.  Calculation  of  wave  height. 


Pass  No. 

»,/</,> 

°J<L> 

V.  dB 

tan  8, 

Aq.  m 

H,.  m 

Sea  truth,  m 

762 

0-798 

0-556 

5-45 

0-0492 

5-5 

791 

0-753 

0-503 

481 

00434 

2-95.  3-25 
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Figure  3.  Part  of  the  digitally  processed  image  for  Seasat  pass  762.  The  satellite  flight  path  is 

from  left  to  right. 

Valenzuela's  curve.  The  error  in  20  is  estimated  to  be  of  the  order  of  10  per  cent,  but 
we  have  no  figure  for  the  error  in  tan  0r.  It  is  unlikely,  however,  that  the  overall  error 
is  less  than  IS  per  cent. 

In  both  cases,  the  calculated  wave  height  is  rather  higher  than  the  measured  value : 
23  per  cent  for  pass  762  and  18  per  cent  for  pass  791.  The  similarity  of  these  two  errors 
may  indicate  some  common  effect  causing  the  overestimation  of  wave  height.  One 
possibility  is  that  either  the  straining  mechanism  or  the  velocity  bunching  is  making 
a  significant  contribution. 

An  expression  for  the  magnitude  of  the  straining  effect’s  contribution  to  the  modu¬ 
lation  transfer  function  is  given  by  Alpers  et  al.  (1981).  which  contains  an  unknown 
parameter  /r.  No  indication  of  the  magnitude  of  n  is  given,  but  even  assuming  ^  equal 
to  zero,  which  gives  the  maximum  possible  straining  effect,  this  effect  is  still  smaller  by 
a  factor  of  three  than  the  tilt  modulation  effect.  Also,  Alpers  et  al.  (1981)  remark  that 
the  tilt  modulation  is  dominant  for  low  to  moderate  states,  and  so  it  is  probably  safe 
to  neglect  it. 

The  velocity  bunching  mechanism  is  linear  for  a  range  of  azimuth  angles  about 
<j>  =  90°.  In  this  range,  Alpers  et  al.  (1981)  give  an  expression  for  the  contribution  to 
the  modulus  of  the  transfer  function  |  Rb  \ 

|K»j  =  * \k\co  cos<p g{0,  <j>), 

where  R  =  range  to  target.  V  =  platform  velocity,  k  =  wave  number  of  ocean  wave. 
cd = angular  frequency  of  ocean  wave,  6= angle  of  incidence,  and  g(0.  $) =0-9994  for 
Seasat.  For  Seasat.  the  range/velocity  ratio  is  128  s.  roughly  an  order  of  magnitude 
greater  than  for  an  aircraft.  Therefore,  as  <f>  deviates  from  90°.  i.e.  as  the  waves  acquire 
an  azimuthal  velocity  component,  the  bunching  mechanism  rapidly  becomes  im¬ 
portant.  However,  the  linear  regime  is  very  narrow  for  large  values  of  R/V  and  is 
defined  by  Alpers  et  al.  (1981)  as  the  range  for  which 

C  =  j  j  x  wave  amplitude  <  0-3. 

For  pass  761  C  is  equal  to  0-5,  and  so  we  are  well  outside  the  linear  regime.  For  pass 
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791.  C  equals  0-64.  However,  the  velocity  bunching  mechanism  is  not  entirely  under¬ 
stood.  For  example,  for  Seasai  pass  834  over  the  JASIN  area,  it  is  known  that  waves 
with  a  significant  wave  height  of  31  m  (Guymer  1980|  were  present  propagating  at 
an  azimuth  angle  4>  roughly  equal  to  60".  but  no  waves  are  visible  on  the  digitally 
processed  imagery.  On  the  other  hand.  Alpers  et  ai  (1981)  refer  to  passes  308  and 
1044  as  showing  azimuth  travelling  swell. 

Clearly,  the  position  is  not  fully  understood  but  it  seems  that  in  the  two  cases 
analysed  in  this  paper,  the  tilt  modulation  mechanism  is  able  to  account  in  large 
measure  for  the  image  modulation.  With  an  airborne  SAR.  the  velocity  bunching 
mechanism  would  be  considerably  less  than  for  a  satellite  system.  Indeed,  it  could  be 
eliminated  by  flying  parallel  to  the  wave  crests. 
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A8STRACT 

The  locations  of  all  frontal  bound¬ 
ary  and  internal  wave  patterns  on  21 
passes  of  Seasat  Imagery  collected  over 
the  eastern  North  Atlantic  Ocean  were 
plotted  on  a  bathymetric  chart  of  the 
region.  These  plots  revealed  that  the 
patterns  always  occurred  over  or  close 
to  a  deep  water  bottom  feature,  such  as 
a  ridge,  seamount,  bank  or  shelf.  At  no 
time  did  a  pattern  occur  over  a  deep 
water  basin.  An  Internal  wave  pattern 
appeared  In  Seasat  SAR  Imagery  67  percent 
of  the  time  Seasat  passed  over  a  deep 
water  feature,  while  frontal  boundary 
patterns  appeared  28  percent  of  the  time. 
The  causes  of  Internal  wave  and  frontal 
boundary  patterns  on  the  Seasat  SAR 
Imagery  have  been  Identified  as  Lee  wave 
generation  and  deep  water  ubwelllng, 
respectively.  Research  which  combines 
oceanic,  subsurface  measures  with  the 
Seasat  SAR  Imagery  Is  definitely  war¬ 
ranted  to  further  study  these  patterns. 


1.  INTRODUCTION 


The  primary  mission  of  the  Seasat  synthetic  aperture  radar 
waves.  Ice  and  other  oceanic  surface  patterns  that  are  a  result 
frontal  boimdarles,  oil  slicks  and  meteorological  occurrences. 


(SAR)  was  to  image  gravity 
of  Interval  waves,  currents. 

In  collecting  over  100  million 


‘Presented  at  the  Seventeenth  International  Symposium  on  Remote  Sensing  of  Environment,  Ann 
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square  kilometers  of  SAfi  imagery  from  space,  Seasat  provided  the  oceanographic  and  remote 
sensing  communities  with  an  abundance  of  data  to  examine. 

Numerous  studies  in  the  past  have  documented  the  presence  of  Internal  wavr  and  frontal 
boundary  patterns  on  SAR  images  in  shallow  {<100  m),  coastal  areas.  An  early  ^easat  study  of 
imagery  collected  during  Rev.  762  between  Scotland  and  Iceland  reported  numerous  large  wave- 
like  patterns  in  deep-water  (>200  m)  regions  {Shuchman  and  Kasischke,  1979).  Here,  we  report 
on  the  continuation  of  this  earlier  study  on  the  appearance  of  large-scale  surface  patterns 
on  Seasat  SAR  imagery  collected  over  the  northeast  Atlantic  Ocean.  These  large-scale  surface 
patterns  have  been  tentatively  identified  as  surface  manifestations  of  internal  waves  and 
regions  of  upwelling  which  result  in  frontal  boundaries. 

We  will  first  present  a  background  section  which  contains  discussions  on  the  Seasat  SAR, 
SAR  imaging  of  the  ocean  surface  and  past  investigations  on  Internal  wave  and  frontal  boundary 
patterns  on  SAR  imagery.  This  is  followed  by  a  section  which  presents  examples  of  deep-water 
internal  wave  and  frontal  boundary  patterns  observed  on  the  Seasat  SAR  imagery  and  discusses 
occurrence  of  these  features  on  the  21  passes  of  Seasat  data  examined  during  this  study. 
Finally,  two  first-order  hydrodynamic  models  are  presented  which  explain  why  the  observed 
patterns  were  detected  by  the  Seasat  SAR. 

2.  BACKGROUND 

The  data  to  be  presented  in  this  paper  were  collected  by  the  Seasat  satellite.  Among 
the  instrumentation  carried  by  Seasat,  which  was  launched  during  June  of  1978,  was  an  imaging 
synthetic  aperture  radar  (SAR).  This  satellite  collected  over  500  oasses  of  SAR  data  before 
suffering  a  catastroohic  power  loss  in  October  of  1978.  The  SAR  on  board  Seasat  was  an  L-band 
(23.5  cm  wavelength)  radar.  It  collected  25  x  25  m  resolution  imagery  with  a  ground  swath- 
width  of  100  km  and  a  length  of  up  to  4000  km,  and  viewed  the  surface  of  the  earth  with  an 
average  incidence  angle  of  20’.  For  a  detailed  description  of  the  Seasat  SAR  system  and  its 
mission,  see  Jordan  (1980);  Beal,  et  al .  (1981),  or  Fu  and  Holt  (1982). 

An  imaging  radar  such  as  the  Seasat  SAR  is  an  active  device  that  senses  the  environment 
with  short  wavelength  electromagnetic  waves.  As  an  active  sensor,  the  Seasat  SAR  provided 
its  own  illumination  in  the  microwave  region  of  the  electromagnetic  spectrum  and  thus  was  not 
affected  by  diurnal  changes  in  emitted  or  reflected  radiation  from  the  earth's  surface.  Addi¬ 
tionally,  the  23.5  cm  wavelength  utilized  by  the  Seasat  SAR  allowed  for  imaging  the  earth's 
surface  through  clouds  and  light  rain. 

The  principle  in  imaging  any  ocean  surface  with  a  radar  is  that  the  backscatter  of  micro- 
wave  energy  (echo)  received  by  the  radar  contains  information  on  the  roughness  characteristics 
(shapes,  dimensions  and  orientations)  of  the  reflecting  area.  Parameters  that  influence  the 
SAR  image  of  the  ocean  surface  Include  the  motion  of  the  scattering  surfaces,  coherent  spec¬ 
kle,  system  resolution  and  non-coherent  integration  as  well  as  the  surface  roughness.  In 
addition,  the  orientation  of  ocean  surface  patterns  with  respect  to  the  radar  "look"  direction 
can  be  influencing  parameters. 

Several  scattering  models  exist  that  attempt  to  explain  ocean  surface  image  formation 
with  synthetic  aperture  radars.  These  models  are  of  two  types:  static  models  that  depend  on 
instantaneous  surface  features,  and  dynamic  models  that  employ  surface  scatterer  velocities. 

Three  static  models  have  been  suggested  to  describe  the  radar  scattering  of  energy  from 
large  areas  on  the  ocean  surface.  These  three  scattering  models  Include:  (1)  the  specular 
point  mode!  which  is  most  aoprooriate  for  small  incidence  angles,  (2)  the  Sragg-Rice  scatter¬ 
ing  model,  described  below,  and  (3)  a  Rayleigh  scattering  model  which  is  often  used  in  terres¬ 
trial  earth  scattering  calculations.  There  Is  general  consensus  within  the  radio-oceanography 
scientific  community  that  a  8ragg-R1ce  scattering  theory  best  explains  the  $ARtobserved  back¬ 
scatter  values  obtained  from  the  ocean  surface  for  incidence  angles  between  20'  and  60’ 
(Shuchman,  et  al.,  1981).  The  Bragg-RIce  scattering  model  is  based  on  a  well  known  phenomena 
in  the  study  of  crystals,  grating,  and  periodic  structures.  If  one  considers  the  random  ocean 
surface  to  be  represented  by  a  combination  (l.e.,  spectrum)  of  periodic  surfaces,  then  the 
spectrum  region  which  satisfies  the  backscatter  phase  matching  condition  will  be  the  main  con¬ 
tributor  to  the  backscatter  cross  section.  Sometimes  In  the  literature,  this  phase  matching 


of  the  small  ocean  Bragg  waves  with  the  radar  electromagnetic  energy  is  termed  a  resonance 
phenomenon;  more  correctly  stated,  it  should  be  termed  a  constructive  interference  between 
the  electromagnetic  and  ocean  waves. 

Pioneering  theoretical  and  experimental  work  by  Wright  (1966)  at  the  Naval  Research 
Laboratory  ( NR L )  demonstrated  the  general  validity  of  a  Bragg  scattering  model  for  an  ocean 
surface  imaged  by  radar.  In  a  series  of  wave  tank  measurements  using  3  and  25  cm  wavelength 
continuous  wave  (CW)  Ooppler  radars,  Wright  demonstrated  that  8ragg  scattering,  that  is, 
transmitted  radar  energy  with  wave  number  K,  interacts  in  a  resonant  or  interference  fashion 
with  ocean  surface  waves  with  wave  number  such  that 

K,  -  2K  sin  ®,  (1) 

where  K*  u  2n/L  and  K  •  2 */x,  (L  and  x  are  the  wavelengths  of  the  surface  waves  and  the 
radar,  respectively),  and  »  is  the  incidence  angle.  Shuchman,  et  al .  (1981)  showed  that  a 
Bragg  scattering  equation  satisfactorily  explained  the  radar  backscatter  return  from  S AR  using 
data  collected  during  the  Marineland  experiment  (for  a  discussion  of  the  Marineland  experi¬ 
ment,  see  Shemdin,  1980).  It  should  be  noted  that  radar  data  of  large  ocean  areas  (1  x  1  km) 
were  averaged  in  that  analysis.  Thus  based  on  the  above  the  principal  radar  reflectivity 
mechanism  of  imaging  ocean  surfaces  is  via  the  capillary  and  small  gravity  waves  which  pro¬ 
duce  Bragg  scattering  (Raney  and  Shuchman,  1978). 

Synthetic  aperture  radars  are  also  sensitive  to  the  motion  of  scatterers  present  in  the 
imaged  scene  (Raney,  1971).  effects  of  scatterer  motion  on  SAR  imagery  may  include;  (1) 
image  displacement,  smearing  and  loss  of  focus  in  the  azimuth  direction,  and  (2)  loss  of  focus 
in  the  range  direction.  Some  of  these  effects  can  be  removed  during  processing  of  the  SAR 
signal  histories  by  making  appropriate  adjustments  to  the  processor  (Shuchman,  1981).  Effects 
vWiich  cannot  be  removed  during  processing  may  reduce  the  detectability  of  gravity  waves,  but 
are  not  expected  to  have  a  large  influence  on  other  ocean  surface  patterns. 

A  variety  of  processes  can  alter  the  surface  Bragg  waves,  resulting  in  a  distinct  pattern 
on  SAR  imagery.  These  include  oceanic  processes  (currents,  gravity  waves,  internal  waves, 
slicks,  local  water  depth  variations,  water  temperature  and  salinity),  climatic  processes 
(wind,  rain  and  air  temperature)  and  man-made  phenomena  (ships,  buoys  ana  oil  spills).  It  is 
the  hydrodynamic  interaction  between  several  of  oceanic  processes  and  a  distinct  bottom  fea¬ 
ture  which  allows  that  feature  to  be  detected  on  SAR  imagery.  Examples  of  these  bottom- 
induced,  deep-water  surface  patterns  will  be  presented  in  this  paper.  For  examples  of  other 
surface  patterns  on  Seasat  SAR  imagery,  see  Gower  (1981);  Beal,  et  al .  (1981);  or  Fu  and  Holt 
(1982). 

Internal  wave  patterns  have  been  detected  on  a  variety  of  remotely-sensed  images.  They 
have  been  detected  on  aerial  and  spacecraft  photography  (Osborne  and  Burch,  1980;  Gower  and 
Hughes,  1979),  Landsat  MSS  imagery  (Apel ,  et  al.,  1974;  1975;  1975a;  1976)  and  aircraft  and 
spacecraft  SAR  imagery  (Brown,  et  al . ,  1976;  Elachi  and  Aoel ,  1976;  Shuchman  and  Kasischke, 
1979;  Gower  and  Hughes,  1979;  Aoel,  1981;  Fu  and  Holt,  1982;  Ford,  et  al.,  1983;  Trask  and 
3riscoe,  1983;  Alpers  and  Salusti,  1983;  and  Hughes  and  Gower,  1983).  Most  of  the  above 
studies  do  not  present  any  direct  evidence  that  the  observed  patterns  are  due  to  internal 
waves,  but  the  patterns  are  widely  believed  to  be  internal  wave  surface  signatures.  An  ex¬ 
cellent  field  study  by  Hughes  and  Gower  (1983;  see  also  Gower  and  Hughes,  1979)  actually  mea¬ 
sured  a  ti daily-driven  internal  wave  packet  as  it  was  being  Imaged  by  an  aircraft  SAR  and 
Seasat.  This  study  verified  that  SARs  can  detect  surface  manifestations  of  internal  wave 
‘ielas. 


An  early  Seasat  study  by  Shuchman  and  Kasischke  (1979)  indicated  the  presence  of  internal 
"wave-like”  patterns  on  Seasat-SAR  images  in  deep-water  regions.  Since  then,  deep-water  in¬ 
ternal  wave  patterns  have  been  reported  on  Seasat  imagery  by  Fu  and  Holt  (1982)  and  SIR-A 
imagery  by  Ford,  et  al.  (1983). 

There  apoears  to  be  little  question  that  ocean  current  or  frontal  boundaries  can  fre¬ 
quently  be  observed  on  SAR  Images,  although  the  reason  for  their  aopearance  is  not  yet  under¬ 
stood.  Current  boundaries  were  first  noted  on  aircraft  SAR  imagery  by  Moskowitz  (1973). 

Since  then,  frontal  boundaries  have  been  extensively  studied  by  both  aircraft  SARs  (Larson, 
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et  a!.,  1976;  Shuchman,  et  al.,  1977;  1979;  Weissman,  et  al.,  1980;  Shuchman,  et  al.,  1981) 
ana  soaceDorne  SARs  (Hayes,  1981;  Shuchman,  et  al.,  1979;  Shemdin,  et  al . ,  1980a;  Mattie,  et 
al.,  1980;  Cheney,  1981;  Lichy,  et  al.,  1981;  Shuchman,  et  al.,  1981c).  Current  velocity 
gradients  (du/dx)  on  the  order  of  10-3  sec-'  appear  to  be  detectable  (Larson,  et  al., 

19761,  and  the  sensing  of  current  boundaries  apparently  occurs  for  a  wide  range  of  wind  speeds 
(3-10  m/ s)  and  air-sea  temperature  differences  (0.1  -  10'C). 

'he  appearance  of  frontal  boundaries  on  SAR  images  was  at  first  thought  to  be  due  to  a 
Ooppler-induced  image  displacement  (Larson,  et  al.,  1976),  but  this  notion  has  been  dis¬ 
counted  by  the  observation  of  current  boundaries  at  all  radar  look  direction  orientations  on 
Seasat  data  (Hayes,  1981).  Although  the  Ooppler  mechanism  may  be  a  contributing  factor  in 
some  cases,  it  is  now  felt  that  the  main  reason  for  the  appearance  of  current  boundaries  on 
CAR  imagery  is  either  a  stress-induced  variation  in  surface  roughness  in  the  area  of  the 
snear,  or  a  difference  in  the  surface  roughness  of  the  water  masses. 

3.  METHODS 

The  data  set  for  this  analysis  was  the  Seasat  SAR  Imagery  collected  over  the  eastern  half 
of  the  North  Atlantic  Ocean.  These  data  were  available  because  of  the  satellite  receiving 
station  which  was  operating  at  Oak  Hangar,  England  for  the  European  Space  Agency.  The  21 
passes  of  imagery  used  in  this  study  are  summarized  in  Table  1. 

'here  were  two  parts  to  the  present  analysis.  The  first  part  Involved  investigating  the 
relationship  between  the  internal  wave  and  frontal  boundary  patterns  and  the  deep-water  bottom 
features  in  this  region.  The  second  part  was  to  develop  a  first  order,  hydrodynamic  model  as 
to  why  the  patterns  appeared  on  the  Seasat  imagery. 

In  order  to  study  the  relationship  between  surface  patterns  on  Seasat  SAR  imagery  and 
the  deep-water  bottom  features,  the  following  procedure  was  used:  first,  the  location  of  the 
ground  coverage  of  the  Seasat  SAR  was  obtained  by  consulting  the  satellite  data  record  (SDR) 
-or  each  pass.  (A  SDR  contains  a  series  of  ground  latitude  and  longitude  coordinates  for 
specific  times  during  a  Seasat  pass.  Most  of  the  Seasat  SAR  images  contained  annotations  of 
the  time  the  data  were  collected.  Therefore,  the  location  of  a  part-'cular  surface  pattern 
could  be  identified  by  noting  the  time  when  the  pattern  occurred  during  that  orbit.)  The 
ground  coverage  of  each  pass  was  then  marked  on  art  overlay  placed  on  the  chart.  Next,  the 
'ocaticns  of  all  internal  waves  and  frontal  boundaries  detected  on  each  pass  were  noted  and 
tnen  these  locations  marked  on  the  overlay.  Finally,  a  surface  pattern  was  said  to  be  asso¬ 
ciated  with  a  bottom  feature  if  its  location  was  within  10  kilometers  of  tnat  feature. 

The  bottom  morphology  of  the  eastern  North  Atlantic  Ocean  is  quite  well  known  due  to  the 
extensive  hydrographic  surveys  conducted  In  this  region  (see  e.g.,  Roberts,  et  al.,  1979). 

"his  area  contains  a  large  number  of  submarine  ridges,  seamounts,  banks  and  shelfs.  Table  2 
lists  the  major  oottom  features  in  this  area,  along  with  the  shallowest  contour  interval 
charted  for  that  feature  and  the  depth  of  the  adjacent  deep  water. 

The  second  part  of  the  analysis  was  performed  by  analyzing  the  environmental  conditions 
which  were  most  likely  present  at  specific  sites  for  a  given  set  of  patterns  observed  in  a 
single  pass  of  Seasat  SAR  imagery,  and  then  formulating  a  first-order  hydrodynamic/ 
electromagnetic  model  which  explains  the  SAR-observed  surface  pattern. 

4.  OBSERVATIONS 

In  this  section,  we  will  first  present  examples  of  internal  wave  and  frontal  boundary 
patterns  fr—  the  study  area.  This  will  be  followed  by  a  summary  of  the  relationship  of  these 
patterns  to  the  deep-water  bottom  features.  Figure  1  is  a  diagram  of  the  positions  of  the 
Seasat  Images  presented  in  this  section. 

The  J?L-oot1cally  processed  SAR  imagery  in  Figure  2  was  collected  during  Rev.  599  (7 
August  1978)  as  Seasat  passed  over  the  Wyv  Ille-Thomson  Ridge.  The  corresponding  bathymetric 
chart  for  this  coverage  is  presented  in  Figure  3.  From  Figure  2,  it  can  be  seen  that  there 
are  numerous  internal  wave  signatures.  Note  on  the  SAR  image  that  the  internal  waves  are  both 
range  (e.g.,  C/01  to  C/02)  and  azimuth  (e.g.,  A9  to  09)  traveling.  The  Internal  waves  in  the 


lower  left  corner  of  the  SAR  image  appear  to  occur  over  the  Ymir  Ridge;  the  Internal  waves  In 
the  middle  of  the  image  appear  over  the  Wyville-Thompson  Ridge;  and  those  at  the  top  of  the 
image  occur  over  the  Faeroe  Sank  Channel. 

Surface  measurements  collected  at  the  time  of  the  Seasat  overpass  Indicate  a  wind  of  6.9 
m/s  from  the  north,  and  a  wave  field  with  a  significant  wave  height  of  1.1  meters,  and  a  dom¬ 
inant  wavelength  of  71  meters  propagating  towards  200 *( T) . 

The  SAR  Imagery  in  Figure  4  was  collected  during  Revolution  958  (1  September  1978)  as 
Seasat  passed  over  the  Icel and-Faeroe  Rise.  The  corresponding  bathymetric  chart  for  this  area 
is  presented  in  Figure  5.  The  Internal  wave  patterns  in  the  top  half  of  Figure  4  occur  in  an 
area  where  the  water  depths  are  between  600  and  1000  meters.  It  is  in  an  area  where  the  ocean 
floor  is  rising  from  the  deep  (>2500  m)  Iceland  Basin  to  the  relatively  shallow  (500  m)  Ice- 
land-Faeroe  Rise.  The  dark  areas  of  Figure  4  (J7  and  02)  are  believed  to  be  atmospheric 
events. 

Surface  measurements  obtained  near  the  test  area  indicates  the  gravity  wave  field  present 
at  the  time  of  Revolution  958  had  a  significant  wave  height  of  1.4  meters,  a  dominant  wave¬ 
length  of  147  meters  and  was  propagating  towards  176 "( T) .  Surface  winds  had  a  speed  of  7.2 
m/s,  but  no  direction  was  given. 

Figure  6  presents  Seasat  SAR  imagery  collected  during  Rev.  762  (19  August  1973).  The 
corresponding  bathymetric  chart  for  this  region  is  presented  in  Figure  7.  This  image  contains 
both  internal  wave  patterns  as  well  as  a  frontal  boundary  pattern.  There  are  two  distinct 
sets  of  internal  wave  packets  present.  One  set  (F5  to  J8)  can  be  seen  on  the  West  Shetland 
Shelf,  while  the  other  lies  along  the  southern  edge  of  the  Wyville-Thompson  Ridge  (C3  to  H4). 
The  frontal-boixidary  pattern  in  this  image  can  be  seen  to  lie  along  the  edge  of  the  Faeroe 
Shelf  (A6  to  89). 

Surface  measurements  obtained  near  the  time  of  the  Seasat  overpass  indicate  a  strong  wind 
(12.5  m/s)  from  the  south  and  a  swell  with  a  significant  wave  height  of  3.0  m/s  and  a  wave¬ 
length  of  210,  propagating  toward  60*(T)  was  present. 

Not  all  surface  signatures  believed  to  be  related  to  deep-water,  bottom  topographic  fea¬ 
tures  are  as  clear  and  distinct  as  the  examples  presented  previously.  Figure  8  was  collected 
by  Seasat  during  Rev.  791  (21  August  1978)  as  ft  passed  over  the  Anton  Dohrn  Seamount.  The 
corresoonding  bathymetry  for  this  imagery  is  presented  in  Figure  9.  The  internal  wave  signa¬ 
tures  collected  over  the  Anton  Dohrn  Seamount  are  quite  subtle.  Three  groups  of  internal 
waves  (36  to  C6,  0/E6  to  0/E8,  and  A8)  occur  over  this  seamount  ,  which  is  quite  steep  sided, 
as  can  be  seen  by  the  narrow  spacing  of  the  contour  lines  in  Figure  9. 

The  sea  truth  collected  near  the  test  area  Indicates  the  winds  were  quite  strong,  with  a 
speed  of  13  m/s  from  the  southwest.  The  surface  gravity  wave  field  had  a  significant  wave 
height  of  3.1  m,  a  dominant  wavelength  of  151  m  and  was  propagating  towards  45’(T), 

Table  3  summarizes  the  occurrences  of  the  internal  wave  and  frontal  boundary  patterns 
over  the  26  deep-water  bottom  features  in  the  study  area.  Table  4  summarizes  the  number  of 
times  each  Seasat  revolution  passed  over  a  deep-water  bottom  feature  and  the  total  number  of 
internal  wave  and  frontal  boundary  patterns  appeared  over  a  bottom  feature  in  that  pass. 

At  no  time  did  an  internal  wave  pattern  or  frontal  boundary  pattern  appear  on  the  Seasat 
imagery  over  a  deep-water  region  without  a  bottom  feature  in  the  immediate  vicinity.  From 
Table  3,  it  can  be  seen  that  67  percent  of  the  time  the  Seasat  SAR  Imaged  the  ocean  surface 
over  a  deep-water  bottom  feature,  an  internal  wave  pattern  was  present.  The  corresponding 
figure  for  frontal  boundary  oatterns  was  23  percent.  The  figure  for  Internal  waves  and 
frontal  boundary  patterns  combined  is  82  percent.  These  figures  suggest  a  strong  correlation 
between  the  bottom  features  and  the  SAR-observed  surface  patterns.  This  relationship  will  be 
further  analyzed  in  the  next  section. 
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5.  HYORODYNAMIC  MODELS 


The  Seasat  SAR  imagery  presented  in  Figure  6  was  used  as  a  base  to  formulate  first-order 
hydrodynamic  models  to  explain  why  the  internal  wave  and  frontal  boundary  patterns  occur. 
Figure  10  combines  a  line  drawing  of  the  internal  waves  and  a  frontal  boundary  observed  in 
Figures  with  the  corresponding  bathymetric  chart  for  this  area.  From  Figure  10,  it  can  be 
seen  that  there  are  numerous  internal  wave  signatures.  The  internal  waves  in  the  middle  of 
the  image  appear  to  be  generated  over  the  Wyville-Thomson  Ridge.  Near  the  lower  right  corner 
of  the  image,  three  wave  packets  are  visible.  It  is  to  be  noted  that  no  waves  were  generated 
on  the  north  side  of  the  ridge. 

The  dominant  bottom  features  in  this  region  are  a  broad,  saddle-like  valley  on  the  east 
side  of  the  ridge,  and  a  90*  sector  of  deep  water  extending  out  from  the  ridge,  as  shown  in 
Figure  11.  The  ridge  has  the  characteristic  scales  of  approximately  10  km  in  both  the  longi¬ 
tudinal  and  transverse  direction.  The  minimax  of  the  bathymetric  saddle  is  at  a  depth  of 
about  450  m. 


3ased  on  the  bathymetry,  it  seems  that  the  wave  packets  are  produced  by  lee  wave 
formation— a  mechanism  similar  to  the  undulation  of  the  jet  stream  as  it  passes  over  a  moun¬ 
tain  ridge.  Figure  12  schematically  illustrates  the  stages  that  are  hypothesized  during  the 
production  of  solitary  Internal  wave  packets.  As  the  strong  semidiurnal  tidal  current  flows 
north  into  the  Faeroe  8ank  Channel,  an  Internal  lee  wave  forms  on  the  northern  edge  of  the 
ridge.  As  the  tidal  flow  goes  to  zero  six  hours  later,  the  lee  wave,  trying  to  maintain  its 
group  velocity  relative  to  the  current,  escapes  south  over  the  ridge  barrier  back  Into  the 
deep  water,  emerging  as  a  packet  of  solitary  waves.  This  phenomenon  has  previously  been 
studied  by  Holbrook,  et  al .  (1983)  in  the  Sulu  Sea.  The  initial  waveform  then  can  be  viewed 
as  a  localized  source  that  Immediately  begins  to  propagate,  and  evolves  into  a  series  of  sol¬ 
itary  waves  while  undergoing  radial  spreading,  encountering  variable  topography  and  slowly 
losing  energy  through  dissipation.  The  disappearance  of  internal  waves  on  the  north  side  of 
the  ridge  Is  probably  because  there  is  no  strong  tidal  current  flowing  south  across  the  ridge; 
most  of  the  tidal  current  may  flow  along  the  Faeroe  Bank  Channel. 

Based  on  the  Seasat  SAR  Image,  the  wave  packets  in  deep  water  over  the  Wyville-Thomson 
Ridge  have  .he  following  properties: 

1.  Wave  packets  are  separated  by  distances  of  the  order  of  14  km. 

2.  Wave  speed  is  estimated  to  be  0.31  m/s  based  on  the  assumption  of  a  semidiurnal  tidal 
origin. 

3.  Wave  packets  in  this  region  contain  about  five  waves  with  packet  width  about  5  km. 

4.  Wavelengths  are  between  500  m  and  1000  m. 


5.  The  crests  are  curved  in  a  horizontal  plane  with  their  convex  sides  pointed  in  the 
direction  of  propagation. 

Although  no  in  situ  measurements  are  available,  it  is  believed  this  area  is  dominated  by 
tidal  currents  with  velocities  of  less  than  0.5  m/s  at  the  surface.  Figure  13  shows  the  pro¬ 
files  of  potential  temperature,  salinity  and  sigma  theta  generated  from  data  collected  over 
the  Iceland-Faeroe  Rise  during  the  Atlantic  Expedition  in  August  1972  (3ainbridge,  1980). 
Station  19  is  the  measurement  closest  to  the  Wyville-Thomson  Ridge  area.  As  evident  from 
Figure  13,  a  mixed  layer  of  40  m  is  dominant  in  this  area  during  the  summer.  Based  on  the 
two-layer  finite  depth  model,  the  linear  wave  speed  is  given  by 


Near  the  Wyville-Thomson  Ridge,  ao/o  *  0.5  x  1 0~3,  Hi  »  40  m,  H?  »  260  m,  2Wk  *  500  m, 
and  thus  the  wave  speed  is  approximately  0.40  m/s,  depending  on'the  water  depths.  This  value 
is  consistent  with  the  estimated  wave  speed  of  0.31  m/s  from  the  SAR  image.  Based  on  the 
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wavelength,  the  wave  amplitude  is  estimated  to  be  between  10  and  15  m.  Note  that  the  effect 
of  shear  is  Important  In  this  case,  since  the  wave  speed  is  slow  compared  with  the  background 
current.  Also,  the  bottom  topography  near  the  West  Shetland  Shelf  is  complex,  which  influ¬ 
ences  the  evolution  of  wave  packets,  as  shown  in  Figure  10.  A  similar  solitary  wave  theory 
developed  byliu,  et  al .  (1983)  can  be  applied  here.  However,  due  to  the  limited  scope  of 
this  investigation,  the  detailed  analysis  and  calculations  will  have  to  wait  until  sea-truth 
data  are  available. 

Liu  (1983)  also  modeled  the  surface  strain  rate  produced  by  shallow  water  (coastal)  in¬ 
ternal  waves.  He  then  used  this  surface  strain  as  an  input  to  an  electromagnetic  model  and 
simulated  a  SAR  image  of  the  internal  wave,  which  compared  favorably  to  the  Seasat  image.  A 
similar  surface  strain  is  expected  to  occur  for  the  deep-water  Internal  waves,  resulting  in 
the  patterns  observed  on  the  JASIN  Imagery. 

The  curved,  dark  line,  as  Indicated  by  the  dashed  line  in  Figure  10,  is  believed  to  be  a 
surface  pattern  caused  by  a  frontal  boundary.  Its  location  appears  to  be  over  the  edge  of 
the  Faeroe  Shelf,  and  closely  matches  the  800  meter  bottom  contour. 

As  shown  In  Figure  13,  the  deep  water  upwelling  results  in  a  frontal  boundary  over  the 
Icel andic-Faeroe  Rise.  The  relationship  of  the  front  to  local  upwelling  topographic  mixing 
may  be  described  by  the  following  process.  Figure  14  shows  the  cross-section  diagram  of  bot¬ 
tom  topography  of  the  Wyvll  le-Thomson  Ridge.  It  appears  that  the  tidal  current  flowing  over 
the  ridge  through  another  saddle-like  valley  on  the  west-side  of  the  ridge  forces  the  deep 
water  Isotherms  up  toward  the  surface.  To  illustrate  this  mechanism  for  frontal  boundary 
generation,  a  schematic  diagram  of  deep  water  upwelling  over  the  Faeroe  Bank  Channel  is  shown 
in  Figure  15.  Therefore,  the  frontal  boundary  seems  to  represent  simply  a  surface  intersec¬ 
tion  of  the  seasonal  thermoeline.  The  combination  of  frontal  slope  and  current  shear  could 
yield  information  on  density  contrast  across  the  front. 

6.  SUW  ARY 

Twenty-one  passes  of  Seasat  SAR  imagery  collected  over  the  eastern  North  Atlantic  Ocean 
were  examined  for  the  presence  of  Internal  wave  and  frontal  boundary  patterns.  All  passes 
contained  such  signatures,  and  the  geographic  position  of  all  signatures  were  found  to  be 
within  10  kilometers  of  a  deep-water  bottom  feature,  such  as  a  seamount,  submarine  ridge, 
bank,  or  shelf.  It  was  determined  that  on  67  percent  of  the  occasions  the  Seasat  SAR  passed 
over  the  geographic  location  of  a  deep-water  bottom  feature,  an  Internal  wave  pattern  was 
present  at  that  location.  The  corresponding  figure  for  frontal  boundary  patterns  was  38 
percent. 

A  first-order  hydrodynamic  model  was  developed  to  explain  the  appearance  of  the  Seasat 
SAR-observed  surface  patterns.  These  models  indicate  that  the  internal  wave  patterns  are 
formed  by  lee-wave  internal  waves,  which  are  generated  as  tidal ly-drlven  currents  flow  over 
the  deep-water  bottom  features.  The  frontal  boundary  patterns  were  shown  to  be  due  to  the 
upwelling  of  deeper  water  to  intersect  the  surface  layer.  This  upwelling  is  generated  by 
tidal  currents  flowing  against  a  deep-water  bottom  feature. 
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Figure  I.  Location  of  Major  Bottom  Features  In  Eastern  North  Atlantic 
and  Positions  of  Seasat  SAR  Images 
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lure  4.  Seasat  SAR  Imagery  Collected  Over  the  Iceland- 
Faeroe  Rise  During  Rev.  958 
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TABLE  1 

SUMMARY  OF  SEASAT  SAR  PASSES  COLLECTED  OVER  EASTERN  NORTH  ATLANTIC  OCEAN 


Seasat 

Time 

Revolution 

Date 

(GMT) 
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08:15 

1044 

8  September  1978 

00:18 

1049 

8  September  1978 

08:27 

1037 

10  September  1973 

00:30 

1149 

15  September  1978 
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TABLE  2 

BOTTOM  TOPOGRAPHIC  FEATURES  IN  THE  NORTHEASTERN  ATLANTIC 
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TABLE  3.  SUMMARY  OF  OCCURRENCES  OF  BOTTOM- RELATED  SURFACE  PATTERNS 
ON  SEASAT  SAR  IMAGERY  COLLECTED  OVER  THE  NORTHEAST  ATLANTIC 
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ABSTRACT 


The  SEASAT  satellite  system,  using  an  onboard  23  an  wavelength  synthetic 

aperture  radar  (SAR),  collected  ~  25  to  40  m  resolution  radar  images  of  the 

ocean  in  100  km  wide  swaths  ranging  In  length  from  ~  300  to  3000  km.  Here  we 

report  results  from  the  18  SEASAT  SAR  passes  during  the  Joint  Air-Sea 

Interaction  (JASIN)  experiment  conducted  off  the  west  coast  of  Scotland  in 

summer,  1978.  These  many  SAR  images;  when  coupled  with  the  intensive  ship, 

buoy  and  aircraft  measurements  of  the  JASIN  experiment;  provide  a  unique 

opportunity  to  assess  the  ability  of  satellite  SAR  to  measure  ocean  surface 

phenomena,  particularly  surface  wave  fields.  In  this  study  we  use  only 

optically  processed  SAR  images.  Although  gravity  waves  of  length  ~  80  to 

300  m  are  often  seen  in  SAR  imagery,  they  are  not  always  seen.  We  find  that 

SAR  resolution  and  waveheight  are  Important  criteria  for  determining  wave 

visibility  and  suggest  that  wind  velocity  is  also.  Comparisons  between  SAR 

and  buoy  estimates  of  dominant  wavelength  and  direction  (including  data  from 

14- 

several  other  experiments)  agree  to  within  about  ±1-3%  and  ±10°  respectively. 
We  use  a  focus  sharpness  algorithm  to  resolve  the  180°  directional  ambiguity 
of  SAR  Image  directional  estimates.  Correlation  of  buoy  measurements  of 
significant  waveheight  (Hjyj)  with  peak  signal  to  noise  ratio  in  Fourier 
transforms  of  SAR  Images  (r  *  0.7),  suggests  that  H^j  can  be  estimated  to  an 
accuracy  of  about  ±1  m.  However,  a  similar,  but  weaker,  correlation  exists 
between  signal  to  noise  and  wavelength,  i.e.  H^j  and  wavelength  are 
correlated  in  the  JASIN  data  set.  SAR  Image  power  spectra  are  in  rough 
agreement  with  buoy  measurements  of  omnidirectional  ocean  waveheight 
spectra  f(K)  and  correspond  less  closely  with  ocean  wave  slope  spectra  ¥'(K). 
However,  the  dominant  wavenumber  (Kpea|(}  in  SAR  spectra  typically  falls  below 
the  dominant  wavenumber  of  corresponding  buoy  spectra.  Further,  the  slope  of 
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SAR  spectra  for  K  <  Kp^  is  typically  steeper  than  corresponding  buoy 
measurements  of¥(K)  while  the  slope  of  SAR  spectra  for  K  >  Kpg^  is  typically 
less  steep  than  corresponding  buoy  spectra  ¥(K).  These  differences  between 
SAR  spectra  and  buoy  measurements  of  ^(K)  form  qualitative  support  for  the 
imaging  mechanisms  developed  by  Alpers,  Ross  and  Rufenach.  Analysis 
algorithms  based  on  this  theory,  but  including  emperical  modifications,  should 
substantially  improve  estimates  of  ^(K.Q)  using  SAR  images. 
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I.  INTRODUCTION 


Our  primary  objective  in  this  paper  is  to  assemble  a  number  of  comparison 
sets  of  ocean  gravity  wave  measurements  based  on  both  synthetic  aperture  radar 
(SAR)  and  surface  buoys.  These  comparison  sets  are  then  used  to  comment  upon 
the  ability  of  synthetic  aperture  radar  to  measure  ocean  gravity  waves  in  the 
period  range  from  about  7  to  19  s.  Comparisons  are  made  with  respect  to 
dominant  wavelength  and  direction  of  the  wave  field,  significant  wave  height 
and  the  directional  wavenumber  spectrum  of  waveheight  variance  *(K). 

SAR  imagery  collected  by  the  SEASAT  satellite  over  the  Joint  Air  Sea 
Interaction  (JASIN)  experiment  area,  west  of  Scotland  provides  an  excellent 
opportunity  to  make  the  desired  comparisons.  During  the  SEASAT-JASIN 
experiment  (August-September,  1978)  observations  of  many  air  and  sea 
parameters  were  made  by  ships,  buoys  and  aircraft  including  the  wind  and  waves 
which  are  of  most  interest  here  (Pollard,  1979).  Swaths  of  four  of  the  18 
SEASAT  SAR  passes  over  the  JASIN  area  are  shown  in  Fig.  1. 

The  23  cm  wavelength  SAR  carried  aboard  the  SEASAT  satellite  collected 
*  25  to  40  m  resolution  radar  Images  of  the  ocean  In  100  km  wide  swaths 
ranging  In  length  from  *  300  to  3000  km  along,  but  displaced  to  the  right  of, 
the  subsatellite  track.  The  SEASAT  SAR  Is  unique,  being  the  only  scientific 
SAR  ever  carried  aboard  a  satellite.  During  the  short  lifetime  of  SEASAT 
(June-September  1978)  some  10®  km^  (equivalent  to  about  20%  of  the  Earth's 
surface)  were  imaged.  A  description  of  the  instrument  Is  given  by  Jordan 
(1980)  and  preliminary  results  from  all  SEASAT  Instruments  were  given  in  the 
June  29,  1979  issue  of  Science.  Beal  jt^^aU  (1981)  and  Gower  (1981)  report 
further  SEASAT  SAR  results  and  Vesecky  and  Stewart  (1982)  assess  and  review 
the  ability  of  SEASAT  SAR  to  sense  ocean  surface  phenomena.  Here  we  report 
results  of  SEASAT  S AR  observations  during  the  JASIN  experiment. 
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SAR  images  collected  by  SEASAT  are  basically  high  resolution  maps  of  the 
radar  reflectance  (backscatter)  of  the  ocean  surface  modified  by  surface 
motion  effects.  Because  the  dielectric  properties  of  the  ocean  surface  are 
relatively  uniform,  variations  in  reflectance  are  due  primarily  to  variations 
in  surface  roughness.  In  particular  the  reflectance  is  due  to  resonant 
backscatter  of  the  \  *  23  cm  radar  waves  from  ocean  surface  waves  of  length 
a  »  \/2  sin  e  -  30  cm,  the  angle  of  incidence  Q  being  about  19  to  26°.  The 
resonant  mechanism  is  typical  of  radiation  scattered  from  a  lattice  and  is 
often  called  Bragg  scatter  (Bragg,  1933).  Because  SAR  uses  the  phase  of  the 
radar  echo  to  achieve  high  azimuthal  resolution  in  the  image  (i.e.  to 
correctly  locate  the  origin  on  a  given  echo)  motion  of  the  ocean  surface, 
which  introduces  a  phase  shift  in  the  scattered  signal,  causes  points  in  the 
Image  to  be  displaced  from  their  true  position.  Such  misplacement  is  aptly 
illustrated  by  moving  ships  being  displaced  from  their  wakes  (e.g.  see  Vesecky 
et  _al_. ,  1982  or  Vesecky  and  Stewart,  1982).  Spatial  modulation  of  surface 
roughness  and  surface  motion  (giving  rise  to  Doppler  shifts)  produce  the 
majority  of  ocean  surface  features  in  SAR  imagery.  Phillips  (1981)  describes 
a  number  of  ways  In  which  wind,  waves  and  current  can  modulate  decimeter  scale 
surface  roughness.  Often  the  difficulty  in  Interpreting  SAR  imagery  of  the 
ocean  Is  to  discover  which  of  the  several  candidates  is  really  the  underlying 
cause  of  the  surface  roughness  modulation  sensed  by  the  radar.  Ocean  gravity 
waves  are  easily  recognized  when  they  are  highly  coherent  and  present  an 
easily  recognized  large  scale  pattern  as  shown  In  Fig.  2. 

The  radar  echo  signal  as  received  on  the  SEASAT  satellite  and  relayed  to 
a  ground  station  can  be  processed  into  an  image  using  any  one  of  a  number  of 
different  algorithms.  Further  these  algorithms  can  be  implemented  using 
either  optical  (analogue)  or  digital  processing  apparatus.  The  relative 
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merits  of  optical  and  digital  processing  are  discussed  by  Vesecky  and  Stewart 
(1982).  In  general  optical  processing  can  be  done  more  quickly  and  at  less 
expense,  while  digital  processing  allows  a  better  quality  image  to  be 
produced.  The  SAR  images  used  in  this  study  have  all  been  processed  optically 
either  at  Jet  Propulsion  Laboratory  (JPL)  or  at  Environmental  Research 
Institute  of  Michigan  (ERIM). 


II.  THE  SEASAT-JASIN  EXPERIMENT 

The  SEASAT-JASIN  experiment  was  conducted  some  four  hundred  km  off  the 
west  coast  of  Scotland  during  the  summer  of  1978.  As  shown  in  Fig.  1  this 
area  of  open  ocean  is  generally  deep  with  numerous  banks  and  sea  mounts  rising 
to  depths  of  i  500  m.  The  oceanic  intensive  area  (OIA),  where  the  majority 
of  buoy  measurments  were  made,  has  a  depth  of  greater  than  1000  m  so  bottom 
topography  does  not  modify  the  deep-water  dispersion  relation  for  ocean 
gravity  waves.  Moderate  seas  prevailed  during  the  SEASAT-JASIN  experiment 
period  with  significant  wave  heights  (H^)  ranging  from  about  1  to  5  m. 

Winds  were  also  moderate,  ranging  from  about  3  to  15  ms"*.  Locations  of  wave 
rider  buoys  in  the  OIA  are  shown  In  Fig.  1  along  with  the  FIA  area,  where  most 
of  the  pitch-roll  buoy  measurements  were  made.  Further  details  are  given  in 
both  Fig.  1  and  Table  I. 

Both  pitch-roll  and  wave  rider  buoys  were  used  to  make  wave 
measurements.  Pitch-roll  buoys  were  deployed  from  both  the  R.V.  Atlantis  II 
(Woods  Hole  Oceanographic  Institution)  and  the  M.V.  Discovery  (Institution  of 
Oceanographic  Sciences)  on  8  and  2  occasions  respectively.  These  buoys 
provide  measurements  of  the  omnidirectional  waveheight  variance 
spectrum  *(w)  as  well  as  estimates  of  the  dominant  wave  direction  and  beam 
width.  The  Atlantis  II  buoy  is  described  by  Stewart  (1977).  Wave  rider  buoys 
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were  moored  by  the  Deutsches  Hydrographisches  Instltut  at  a  number  of 
locations  as  shown  in  Fig.  1.  These  buoys  measured  only  *(u)  and  provided  no 
directional  information.  The  SAR  data  collected  over  the  JASIN  area  and  used 
in  this  study  was  optically  processed  into  images  at  both  the  Jet  Propulsion 
Laboratory  (JPL)  and  the  Environmental  Research  Institute  of  Michigan 
(ERIM).  The  JPL  images  were  the  standard  SEASAT  SAR  data  product  having  a 
resolution  of  25  to  40  m.  The  ERIM  images  were  specially  processed  to  obtain 
focus  information  as  described  below  and  to  achieve  consistent  25  m 
resolution.  An  example  of  JPL  imagery  is  shown  in  Fig.  2  where  a  field  of 
about  170  m  waves  is  present.  Examples  of  ERIM  imagery  are  given  by  Kasischke 
et__al_.  (1982). 

During  the  experiment  period  (4  August  to  11  September,  1978)  18  SEASAT 
SAR  passes  were  made  over  various  portions  of  the  JASIN  area.  Table  I 
summarizes  the  data  collected  during  these  passes.  On  16  of  the  passes 
concurrent  or  nearly  concurrent  buoy  measurements  were  made.  On  three 
occasions  a  buoy  was  within  the  SAR  swath.  The  distances  between  the  SAR 
swaths  and  buoy  locations  in  the  other  cases  are  noted  in  Table  I.  Wind 
measurements  were  also  made  concurrent  with  most  of  the  SEASAT  SAR  passes  as 
noted  In  Table  I.  The  wind  speeds  are  averages  of  observations  by  3  to  5 
surface  instruments  within  the  OIA  area  of  Fig.  1  and  hence  near  the  wave  buoy 
locations. 

III.  VISIBILITY  OF  OCEAN  WAVES  IN  SAR  IMAGES 

Although  ocean  gravity  waves  of  length  -  80  to  300  m  are  often  seen  in 
SAR  Iruuges,  they  are  not  always  seen.  Wave  visibility  is  not  simply  a  case  of 
Inefficient  wave  height.  For  example,  one  notes  from  Table  I  that  waves  with 
*  3.6  m  In  a  6.6  ms”^  wind  were  not  seen  In  SAR  images  from  orbit 
1006.  Why  is  this  so? 
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We  suggest  that  three  criteria  are  important  in  determining  wave 
visibility: 

1.  SAR  resolution  (Including  the  degrading  effects  of  wave  orbital 
motion)  must  be  sufficiently  good  that  the  Nyquist  sampling  criterion,  i.e.  at 
least  two  resolution  cells  per  projected  ocean  wavelength,  be  satisified  in 
both  the  azimuth  direction  (parallel  to  the  SAR  flight  path)  and  the  range 
direction  (perpendicular  to  the  SAR  flight  path). 

2.  Hj/3  greater  than  about  1  m. 

3.  Wind  speed  greater  than  a  few  ms"1. 

In  this  study  we  have  used  only  SAR  images  processed  optically  at  the  Jet 
Propulsion  Laboratory  (JPL )  and  the  Environmental  Research  Institute  of 
Michigan  (ERIM). 

Alpers,  Ross  and  Rufenach  (1981)  point  out  that  orbital  velocities 
associated  with  surface  gravity  waves  produce  a  smearing  effect  on  SAR  imagery 
degrading  resolution,  especially  along  the  azimuth  (y)  direction.  They  also 
note  that  the  ocean  surface  coherence  time  (xg)  affects  resolution  (Raney, 
1980)  as  do  adjustments  (AV)  in  the  azimuthal  focus  relative  to  the  nominal 
focus  for  land.  A1  pers  _et__al_.  (1981)  argue  that  the  coherence  (xs)  term  can 
usually  be  neglected  and  we  do  so  here.  If  Ay  is  the  land  (unsmeared)  azimuth 
resolution  (~  25  to  40  m  for  SEASAT  SAR),  then  the  degraded  azimuth 
resolution  Ay*  for  N  looks  is  given  by 


2 

Ay*  -  NAyU  ♦  N'2[-5l(-Au»2g(e,*)cos(K  '*r  +  a')  +  (2VAV/R))]2  +  (-J-)2}^ 


(1) 


—  1 

where  T  is  the  one-look,  full  bandwidth  integration  time  (»  2.3  s  for  SEASAT 

SAR),  \  is  radar  wavelength,  r  the  location  of  the  target  in  the  (x,y) 

plane,  a*  *  tan“*(tan8  sin*)  and  A,  uj,  and  K  are  the  amplitude,  radian 

frequency  and  vector  wavenumber  of  the  dominant  ocean  wave.  V  is  the 

satellite  velocity  relative  to  the  Earth  and  aV  is  an  azimuthal  focus 

adjustment  parameter  which  we  take  to  be  zero  here  since  no  focus  adjustment 

was  made  in  imaging  the  data  used  here.  The  function 

g(8,«)  =  (sin2e  sin2*  +  cos2e^  ;  other  variables  are  defined  in  Fig.  3.  In 

★ 

our  calculation  of  ay  we  have  taken  1/u  as  the  typical  value  of  the  cosine 
★ 

term.  Knowing  ay  and  the  dominant  ocean  wavelength  along  the  y- 

di recti  on  a^  ■  A/cos*,  the  number  of  SAR  resolution  elements  or  samples  n^ 

along  the  azimuth  direction  can  be  estimated.  The  Nyquist  sampling  criterion 

(Bracewell,  1979)  demands  ny  >  2  for  the  waves  to  be  adequately  sampled.  Jain 

(1978)  and  Vesecky  _et__aj_.  (1981)  have  developed  wave  visibility  criteria  along 

these  lines.  While  an  analagous  criterion  nx  >  2  must  also  be  met  along  the 

range  direction,  resolution  degradation  along  the  x-di recti  on  is  much  less 

severe  and  we  have  neglected  it  here. 

Because  of  the  roles  played  by  waveheight ^and/or  slope  in  wave  imaging 

mechanisms  (Alpers  _et__al_. ,  1981)  one  expects  wave  visibility  to  decrease  with 

decreasing  waveheight  H^3  and  Indeed  we  find  waves  are  seldom  imaged 

for  ~  1  m.  This  waveheight  criterion  is  probably  set  by  the  influence  of 

slopes  on  the  wave  imaging  mechanism,  i.e.  on  the  spatial  modulation  of  the 

normalized  radar  cross  section  a  .  If  slopes  are  the  dominant  mechanism 

o 

modulating  o  and  if  slopes  are  small,  then  so  too  is  Ao  .  If  a o  is  less 

Oj  0  0 

than  the  noise  fluctuations  of  aQ  in  the  image,  the  waves  are  invisible. 
Interestingly,  theory  implies  that  there  may  be  both  a  minimum  below 
which  waves  are  Invisible  and  a  maximum  n  above  which  waves  traveling  near 
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the  azimuth  direction  are  not  imaged  because  azimuthal  resolution  is  severely 

degraded  by  large  wave  orbital  velocities.  The  two  criteria 

n  >2  and  H.  ,,  2  1  m  are  applied  to  the  SEASAT-JASIN  data  set  in  Fig.  4.  In 
y 

calculating  ny  N  was  set  at  four  and  Ay  at  6.25  m  corresponding  to  4  look 
optical  processing.  Where  possible  ocean  wave  data  are  drawn  from  buoy 
measurements  given  in  Table  1.  Cases  where  SAR  estimates  were  used  for  4  are 
noted  in  Fig.  4. 

Since  surface  winds  raise  the  ~  30  cm  ocean  waves  to  which  the  radar  is 
sensitive,  we  expect  wave  images  to  disappear  when  wind  speed  falls  below  a 
few  ms"*,  and  the  sea  becomes  calm.  A  wind  speed  of  3.5  ms"*  is  typical  of 
speeds  required  to  maintain  a  saturated  spectrum  of  gravity  waves 
for  A  £  30  cm  (Vesecky  and  Stewart,  1982).  A  SAR  image  of  a  calm  region 
containing  no  visible  waves  is  given  by  Beal  (1981).  During  the  SEASAT-JASIN 
experiment  waves  were  detectable  over  a  range  of  wind  speeds  of  from  3.5  to 
15.2  ms~*  while  for  wind  speeds  of  from  6.6  to  12.9  ms"*  waves  were  not 
detected.  Thus  wind  speed  was  apparently  not  a  factor  during  the  SEASAT-JASIN 
experiment  mainly  because  wind  speeds  were  always  above  3.5  ms"*. 

IV.  SAR  IMAGE  ANALYSIS  METHODS 

SAR  image  analysis  done  here  is  tied  to  a  very  simple  assumption,  namely, 
that  SAR  image  intensity  fluctuations  61  are  proportional  to  ocean  surface 
height  fluctuations  6h.  This  assumption  implies  that  the  Fourier  power 
spectrum  of  the  image  intensity  distribution  |Fj(K,e)|  3  |FtI(x,y)}|  is 

proportional  to  the  directional  waveheight  spectrum  *(K,y)  of  the  ocean 
surface  corresponding  to  the  image.  As  will  become  evident,  the  truth  is  not 
so  simple.  However,  this  assumption  is  not  far  wrong  and  provides  a 
convenient  working  hypothesis.  Using  the  assumption  dominant  wavelength  ana 


direction  can  be  estimated  simply  by  noting  the  values  of  K  and  u  where 
|F(K,a)|  peaks.  We  discuss  below  in  section  VI  how  the  constant  of 
proportionality  can  be  estimated  and  thus  i!(K,9)  estimated. 

Several  complicating  factors  which  qualify  the  simple  assumption  can  be 
summarized  in  the  following  equation  where  these  factors  enter  as  linear 
transfer  functions: 

|Fj  (K,0 )  |2  =>  jH|2  •  |  B  | 2  •  |  R  j 2  •  *(K,9)  .  (2) 

These  transfer  functions  H,  B  and  R  are  related  to  the  point  response  of  the 
SAR  system,  background  (non-wave-related)  fluctuations  in  the  image  and  the 
radar  wave-ocean  wave  interaction  (imaging)  mechanism  respectively.  The 
magnitude  signs  are  used  since  the  transfer  functions  are  in  general  complex. 

As  with  all  measurements ,  the  response  of  the  measuring  instrument  is  not 
perfect.  Thus  the  SAR  image  of  a  point  target  is  spread  over  some  region  of 
the  image  characterized  in  size  by  the  resolution  cell  dimensions  (Ax  and  Ay) 
and  in  shape  by  the  point  spread  function.  To  preserve  details  SAR  images  are 
digitized  at  sampling  intervals  small  compared  to  Ax  and  Ay.  When  an  image  is 
Fourier  transformed  and  the  energy  spectrum  calculated,  one  expects  the 
spectral  power  level  to  decrease  significantly  -(roll  off)  as  K  increases 
approaching  2ji/2Ax.  Here  we  take  ax  =  Ay  since  SAR  imaging  schemes  usually 
strive  toward  this  condition.  The  transfer  function  H(K,e)  can  be  used  to 
characterize  the  effects  of  the  point  spread  function  in  the  wavenumber  domain 
(for  reference  see  Bracewell,  1979  or  Champeney,  1973). 

Perusal  of  SAR  images  of  the  ocean  quickly  reveals  that  there  are 
numerous  features  on  all  size  scales  resolved.  Many  of  these  features  are 
clearly  not  related  to  ocean  gravity  waves  and  hence  constitute  a  background 
against  which  the  wave  images  are  viewed.  It  is  well  known  that  the  radar- 

backscatter-reflectance  of  a  uniform  rough  surface  is  dependent  on  angle  of 
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incidence.  This  fact  along  with  the  effects  of  antenna  pattern  and  automatic 
gain  control  (AGC)  introduce  a  large  scale  trend  in  image  brightness 
perpendicular  to  the  SAR  flight  path  and  hence  a  low  K  background  component  in 
the  image  transform.  Although  some  compensation  for  this  variation  can  be 
introduced  in  the  data  processing,  some  effects  remain,  especially  in 
optically  processed  images,  e.g.  see  Vesecky  and  Stewart  (1982).  We  lump  all 
effects  related  to  background,  i.e.  non-wave-related,  fluctuations  in  a  single 
transfer  function  8(K,0).  Clearly  this  is  a  crude  model  for  these  effects, 
but  it  proves  to  be  useful. 

The  physical  mechanism(s)  which  allow  ocean  gravity  waves  to  be  imaged 
are  discussed  by  a  number  of  authors  and  summarized  by  Alpers,  Ross  and 
Rufenach  (1981)  and  Vesecky  and  Stewart(1982) .  The  three  principal  mechanisms 
considered  are  as  follows: 

1)  tilt  mechanism  (Rt)  in  which  the  change  in  large  scale  surface 
slope  along  the  length  (A  ~  10' s  to  100's  of  m)  of  a  gravity  wave  changes  the 
angle  of  incidence  at  which  the  short  (A  ~  30  cm)  radar-resonant  waves  are 
viewed  by  the  radar  and  hence  spatially  modulates  the  radar  reflectance  as  a 
function  of  position  along  the  wave 

2)  hydrodynamic  mechanism  (R^)in  which  the  variation  in  acceleration 
with  position  along  a  long  gravity  wave,  i.e.  along  the  direction  of  wave 
travel,  induces  a  change  in  energy  density  of  the  radar-resonant  waves  and 
again  spatially  modulates  the  radar  echo  strength  as  a  function  of  position 
along  the  wave 

3)  velocity  bunching  mechanism  (R^)  in  which  the  variation  in  wave 

orbital  velocity  with  position  along  a  long  gravity  wave  causes  the  SAR 

imaging  process  (which  assumes  a  stationary  land  surface)  to  systematical ly 

misplace  image  brightness  along  the  azimuth  direction  thus  spatially 
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modulating  SAR  image  brightness  in  a  manner  corresponding  to  the  variation  in 
wave  orbital  velocity. 

Assuming  these  imaging  mechanisms  to  be  linear  we  can  write 

| R (K  , o )  | 2  =  | | 2 | Rh | 2 | Rb j 2 .  As  pointed  out  by  Alpers,  Ross  and  Rufenach 

(1981)  we  can  expect  the  tilt  and  hydrodynamic  mechanisms  to  be  modeled 

reasonably  well  by  linear  transfer  functions,  but  the  velocity  bunching 

mechanism  is  linear  only  over  a  small  range  of  directions  nearly  perpendicular 

to  the  SAR  flight  path.  Alpers,  Ross  and  Rufenach  have  estimated  the  K 

dependence  of  |Rt|^,  | R^  | ^  and  j R^ | ^  to  be  j R 1 1 2  ~  K^;  [ | 2  ~  to 

(u/V)K  for  n,  the  relation  time  constant,  increasing  from  0  to 
2  .2- 

n  >  w  and  ( R b |  ~(Kw'v  where  w  is  the  ocean,  wave  radian  frequency.  On  the 

2  CL 

basis  of  these  transfer  functions  alone  we  would  expect  |Fj|  ~  K  'i  where 
a  -  2  to  3.  This  K  dependence  of  Rb  applies  only  to  the  range  of  wave 
directions  near  O  =  90°  where  the  velocity  bunching  mechanism  is  thought  to 
be  linear. 

The  transfer  function  H  can  be  estimated  analytically  by  knowing  the 

radar  system  and  image  processing  parameters  as  done  by  Beal  (1981).  In  this 

study  we  have  taken  an  emperical  approach  using  the  variance  spectrum  of  a 

given  SAR  image  to  estimate  the  product  of  H  and  B.  Our  technique  is  to  first 

2 

consider  the  raw  image  spectrum  |Fj(K,9)|  finding  the  range(s)  of  9  where 

wave  energy  is  present.  If  we  further  assume  that  H  and  B  are  isotropic,  then 

|  F j  (K  ,9 )  j  averaged  over  directions  away  from  the  dominant  wave  direction(s) 

represents  background  fluctuations  B  modified  by  the  SAR  image  response  H  and 

2  2 

we  can  thus  estimate  |H|  *|B|  via 

I H | 2 •  |  B  j 2  -  < | Fj (K ,9  )  j  2>  (3) 

where  the  average  is  taken  over  9  away  from  dominant  wave  directions.  Having 

estimated  | H | ^  times  | B | ^  we  can  remove  their  effect  to  obtain  the  corrected 
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or  rectified  image  K  spectrum  |FR(K,tf)|^  =  |Fj(K,e)^/  | H [ 2  •  | B | 2 .  Because 

radar  system  and  processing  parameters  may  change  from  scene  to  scene  we 
2  2 

estimate  |H|  •  |B|  separately  for  each  image. 

The  omnidirectional  spectra  presented  in  this  paper  (Figs.  7  and  11)  were 
obtained  by  using  the  rectified  spectrum  |FR|^  averaged  over  directions  near 
(within  about  ±15°)  the  dominant  wave  direction(s).  Thus  these  spectra  are 
not  truly  omnidirectional.  However,  to  have  included  background  fluctuations 
from  directions  well  away  from  the  dominant  wave  direction  would  only  have 
added  noise.  We  .consider  |F^|^  as  shown  in  Fig.  7  to  be  our  best  present 
estimate  of  ¥(K)  for  comparison  with  buoy  estimates  of  this  quantity.  Note 

p 

that  we  have  not  made  any  adjustments  for  the  wave  imaging  mechanism  |Rp 
since  we  feel  it  is  not  sufficiently  well  known  at  this  time.  Rather  we  hope 
that  comparisons  of  | FR ) 2  with  buoy  measurements  of  ¥(K)  will  help  clarify 
which  if  any  of  the  proposed  SAR  imaging  mechanisms  correspond  to  reality. 

We  consider  results  from  SEASAT  orbit  547  (see  Fig.  1)  to  illustrate 
several  techniques  of  SAR  image  analysis  as  well  as  methods  for  displaying  the 
results.  In  Fig.  5  we  show  an  optical  Fourier  transform  (OFT)  produced  at 
ERIM.  In  this  method  the  coherent  optical  image,  formed  by  a  normal  SAR 
processing  optical  train,  is  run  through  a  converging  lens  to  obtain  the  two- 
dimensional  Fourier  transform  and  a  piece  of  film  in  the  focal  plane  records 
the  intensity  spectrum  which  corresponds  to  | F j | 2  discussed  above.  Practice 
has  shown  that  weakly  Imaged  waves  are  often  more  easily  detected  in  a 
photographic  display  such  as  Fig.  5. 

Digital  Fourier  transforms  (DFT's)  of  SAR  images  are  generated  by  first 

digitizing  a  SAR  image  and  then  performing  the  Fourier  transform  using  the 

fast  Fourier  transform  algorithm  on  a  digital  computer.  Generally  only  the 

2  2 

magnitude  squared  or  variance  spectrum  |Fj(K,tf)|  or  |FR(K,«)|  is 
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examined.  In  Fig.  6  we  see  a  contour  plot  of  | F . (K  ,K  )  j  performed 

1  X  Jr 

digitally  at  ERIM  using  SEASAT  data  from  orbit  547  (see  Fig.  1).  Fig.  7  shows 
another  method  of  displaying  OFT  results,  which  have  been  corrected  as 
discussed  above,  i.e.  ! FR ( ^  is  displayed.  Generally  OFT’s  have  been  used  for 
quick-look  analysis  and  to  detect  weakly  imaged  waves  while  DFT's  have  been 
used  for  more  quantitative  analysis. 

Practice  thus  far  indicates  that  the  threshold  for  detecting  waves  in  SAR 
images,  i.e.  estimating  dominant  wavelength  and  direction,  differs  between  the 
OFT  and  DFT  techniques.  However,  this  difference  is  not  fundamental  and  does 
not  always  hold.  OFT's,  optically  displayed  as  in  Fig.  5,  are  apparently  more 
sensitive  because  they  include  data  from  a  relatively  large  image  area. 

(OFT's  usually  use  a  relatively  small  image  area  to  avoid  excessive 
computation.)  Also  the  photographically  displayed  OFT  allows  the  eye  and 
brain  to  perform  the  pattern  recognition.  The  | H | 2  *  ( B [ 2  correction 
technique,  described  above,  aids  in  detecting  weakly  imaged  waves  in  DFT's, 
i.e.  waves  not  detected  in  displays  of  | F j | 2  are  evident  in  displays  of 
|Fr|2.  Waves  not  detectable  In  an  | F j | 2  plot  of  orbit  599  data  were  detected 
when  ( Fr | 2  was  computed  and  plotted. 

An  important  measurement  which  is  derived  from  OFT  or  DFT  data  is  the 
peak  signal  to  noise  (S/N )  or  peak  to  background  ratio  (PBR)  for  a  given 
image.  We  shall  use  (S/N)  to  denote  the  peak  signal  to  background  noise  level 
of  | FR  | 2  as  shown  in  Fig.  7.  The  PBR  ratio  refers  to  OFT  results  displayed 
photographically  and  denotes  the  ratio  of  film  density  at  the  peak  point 
(corresponding  to  the  dominant  wavelength  and  direction)  to  the  film  density 
at  a  specific  background  location  in  the  OFT  display,  e.g.  near  the  upper 
lefthand  corner  of  Fig.  5. 


V.  COMPARISON  OF  SAR  ESTIMATES  WITH  BUOY  MEASUREMENTS  FOR  THE 


WAVELENGTH,  DIRECTION  AND  BEAMWIDTH  OF  THE  DOMINANT  WAVES 
Comparisons  of  SAR  estimates  of  ocean  waveheight  spectra  with  buoy 
measurements  of  the  same  quantity  serve  two  purposes.  First,  they  allow  one 
to  judge  the  usefulness  of  SAR  observations  in  making  wave  measurements. 
Second,  the  comparisons  shed  light  on  the  physical  mechanisms  which  allow 
waves  to  be  imaged  by  a  SAR.  Better  understanding  of  the  imaging  mechanism 
should  lead  to  more  accurate  SAR  estimates  through  improvement  of  the  analysis 
algorithms,  i.e.  progress  beyond  the  simple  assumption  discussed  above.  Table 
1  summarizes  the  experiment  parameters  and  many  of  the  results  of  the  SEASAT- 
JASIN  experiment. 

It  is  clear  that  buoy  and  SAR  measurements  are  fundamentally  different  in 
that  the  SAR  measurements  are  largely  spatial  averages  at  a  point  in  time 
while  buoy  measurements  are  mainly  temporal  averages  at  a  point  in  space.  To 
compare  the  two  we  assume  that  the  ocean  surface  is  statistically  homogeneous 
and  stationary  over  scales  of  ~  20  km  in  space  and  ~  0.5  hr  in  time  so  that 
spatial  and  temporal  averages  are  statistically  equivalent,  i.e. 

'll  (u,  e)  du  »  i  (K,  9)  KdKThe  omnidirectional  slope  spectrum  4?  *  (K )  is  derived 
from  the  wave  height  spectrum  simply  by  multiplication  by  K2,  i.e. 
r(K)  «  K2'?(K)  (Kinsman,  1965). 

Estimates  of  dominant  wavelength  and  direction  are  obtained  simply  by 

noting  the  (K,9)  location  of  the  wave-related  maximum  in  | F j  | 2  or  |FR|2. 

Two  estimates  of  Kmav  and  9  „  were  made;  one  using  ERIM  OFT  data 

max  max 

corresponding  to  (Fj|2  and  one  using  Stanford  DFT  data  corresponding  to 

|FR|  .  In  making  comparisons  of  dominant  wavelength  and  direction  we  have 

simply  averaged  these  two  SAR  estimates.  As  one  would  expect  for  random 

errors,  in  the  two  estimates,  the  average  of  the  OFT  and  DFT  estimates 
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compares  more  favorably  with  buoy  measurements  than  either  SAR  estimate  taken 
alone.  This  leads  one  to  conclude  that  the  corrections  introduced  to  produce 
| FR | 2  do  not  make  emphatic  improvements  in  SAR  estimates  of  and  emax. 
However,  K  spectra,  such  as  those  in  Figs.  7  &  11,  were  substantially  improved 
by  these  corrections  especially  in  terms  of  distinguishing  between  the  ocean- 
wave-related  component  and  the  background. 

Fig.  8  shows  dominant  wavelength  comparisons  for  the  SEASAT-UASIN 
experiment  and  includes  2  similar  comparisons  from  the  DUCK  EX  experiment  and 
10  from  GOASEX  (Gonzalez,  et_al_. ,  1981).  Using  all  these  data  we  find  SAR 
estimates  agree  with  buoy  estimates  within  an  average  error  of  about  ±14%.  We 
also  note  that  the  SAR  estimates  of  wavelength  are  biased  high  by  about  10%. 
The  reason  for  this  bias  is  discussed  below.  A  simple  scaling  error  does  not 
appear  to  be  the  cause  since  the  same  bias  is  found  in  data  imaged  and 
analyzed  independently  by  three  different  research  groups  —  ERIM  and  Stanford 
(this  paper)  and  Gonzalez  et  al.(1981). 

In  Fig.  9  we  compare  SAR  estimates  of  wave  direction  with  buoy 
estimates.  Results  from  GOASEX  and  DUCK  EX  are  included  as  in  Fig.  8.  The 
SAR  measurements  show  no  significant  bias  and  the  mean  difference  between  the 
two  estimates  is  ±10°.  In  this  error  calculation  approximate  buoy 
measurements  (marked  *  in  Table  1  and  shown  by  diamonds  in  Fig.  9)  were 
excluded  and  hence  only  4  JASIN  cases  were  used. 

SAR  estimates  of  wave  direction  based  on  the  simple  assumption  discussed 

above  contain  a  180°  ambiguity,  l.e.  one  cannot  distinguish  between  waves 

traveling  along  £  and  those  traveling  along  (4  +  180°).  Shuchman  and 

Zalenka  (1978)  have  shown  that  this  ambiguity  can  be  resolved  using  SAR  data 

alone  by  special  focusing  of  a  SAR  image  during  the  imaging  process.  The 

difference  in  the  point  of  sharp  focus  for  waves  relative  to  the  sharp  focus 
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point  for  land  reveals  the  direction  of  wave  travel,  either  toward  or  away 
from  the  radar.  By  applying  this  procedure  to  the  JASIN  data  set,  we  are 
freed  from  relying  on  ad  hoc  assumptions  (such  as  waves  always  traveling 
toward  coasts  or  down  wind)  to  resolve  the  180°  ambiguity  in  dominant  wave 
di recti  on. 

As  shown  in  Fig.  6  SAR  images  yield  information  on  the  wave  directional 

distribution.  Beside  the  dominant  wave  direction  discussed  above  a  "beamwidth" 

for  the  distribution  can  also  be  calculated  by  noting  the  full  width  at  half 

maximum  in  a  SAR  image  directional  distribution  such  as  shown  by  Vesecky  et_ 

al .  (1981  and  1982).  The  SAR  directional  distribution  is  calculated  by 

2 

averageing  |F^(K,9)|  over  an  octave  range  of  K  centered  on  the  dominant 
wavenumber.  The  SAR  estimates  of  beamwidth  are  remarkably  close  to  the  pitch- 
roll  buoy  estimates  In  four  cases  —  average  difference  between  the  two  is 
~4-<’.  A  fifth  case  (orbit  834)  was  difficult  to  interpret  since  the  peak  SAR 
response  was  only  a  factor  of  two  above  the  noise  level.  We  note  here  that 
SAR  directional  distributions  often  show  considerable  structure.  Vesecky  et 
al .  (1980)  provide  evidence  that  at  least  some  of  this  structure  is  physically 
real.  Hence  SAR  measurements  of  wave  directional  distributions  having 
multiple  maxima  could  be  useful  in  locations  where  wave  sensor  arrays  are 
impractical . 

VI.  SAR  ESTIMATES  OF  SIGNIFICANT  WAVEHEIGHT 

As  discussed  above  and  Illustrated  In  Fig.  7,  the  signal  to  noise  ratio 

S/N  (at  the  peak  of  the  wave  K  spectrum)  can  be  calculated  from  digital 

transforms  of  SAR  data.  This  ratio  corresponds  to  the  peak  to  background 

ratio  (PBR)  obtained  from  optically  transformed  SAR  data  as  shown  In  Fig.  5. 

Nine  cases  during  the  SEASAT-JASIN  experiment  contain  both  SAR  (S/N)  and  buoy 
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estimates.  In  Fig.  10  we  plot  these  nine  cases  and  note  a  clear 
correlation  between  (S/N)  and  H^—correlation  coefficient  r  =  0.7.  A  linear 
least  mean  squares  fit  to  these  data  points  yields  (S/N)  =  0.02  +  1.64 
with  a  coefficient  of  determination  R2  =  0.53.  Based  on  this  data  set  the 
linear  fit  estimates  the  value  of  H^/3  with  a  mean  error  of  about  +1  m.  The 
±1  m  error  compares  well  with  the  accuracy  of  waveheights  determined  from  SAR 
data  using  observations  of  speckle  diversity  (Jain,  Medlin  and  Wu,  1982). 
Further,  the  correlation  of  (S/N)  and  Hj^  could  be  used  to  produce  SAR 
estimates  of  ¥(K)  which  y-ield  absolute  values  along  the  ordinates  as  well  as 
the  abscissae  of  Figs.  7  and  11,  i.e.  no  normalization  along  the  ordinates 
would  be  required.  However,  we  note  that  in  the  JASIN  data  set  H-,  n  and  A  are 

ivnd  h«*nc  c.  ^  S/ J 

themselves  correlated^  though  less  strongly  (R2  =  0.40)  (.£/*-  )  H  , 

VII.  COMPARISON  OF  SAR  IMAGE  SPECTRA  WITH  WAVEHEIGHT  AND 

WAVE  SLOPE  SPECTRA  MEASURED  BY  BUOYS 

In  Fig.  11  we  show  six  cases  in  which  SAR  images  from  SEASAT  and 

concurrent  buoy  spectra  were  collected  during  the  SEASAT-JASIN  experiment.  An 

additional  case  is  shown  in  Fig.  7.  The  SAR  spectrum  in  each  case  corresponds 

to  |Fr|2,  i.e.  the  correction  discussed  above  is  included.  In  several  cases 

the  buoy  was  located  within  the  SAR  image  while  in  other  cases  the  two  are 

displaced  as  noted  in  Table  1.  The  omnidirectional  waveheight  spectrum 

M00(K),  which  is  equivalent  to  Hf(K),  Is  shown  by  the  dotted  line  and 

corresponds  to  the  calibration  on  the  ordinate.  These  values  of  MQ0(K)  were 

measured  by  either  waverider  or  pitch-roll  buoys  as  noted  in  Table  1.  To 

represent  the  shape  of  the  wave  slope  spectra  *'(K)  we  have  simply 

multiplied  the  MQ0(K)  curve  by  K2  and  normalized  the  result  to  the  peak  level 

of  the  MQ0  curve.  The  SAR  spectra  are  normalized  along  the  ordinate  to  the 
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peak  level  of  MQ0.  No  normalization  was  done  along  the  abscissa.  As 
discussed  above  it  should  be  possible  to  make  a  direct  comparison  of  the  SAR 
and  buoy  estimates  of  MQ0  without  any  normalization  by  using  the  (S/N)  vs 
correlation. 

One  of  the  most  salient  features  of  the  comparisons  is  that  the  peak  of 
the  SAR  spectrum  is  often  displaced  (toward  lower  K)  from  the  peaks  of  MQ0  or 
K^Mqo  as  measured  by  buoys.  Ths  bias  of  SAR  spectra  toward  longer  wavelengths 
is  also  evident  in  Fig.  8.  Although  there  is  only  a  small  range  of  K  above 
the  peak  value  of  K  we  can  make  a  rough  estimate  of  the  slope  of  a  power  law 
fit  (K^ )  to  | FR | ^  for  K  >  Kpeak*  For  the  cases  shown  in  Figs.  7  and  11  the 
average  slope  in  this  region  is  <p>  — 2.3.  For  the  same  cases  the  average 

slope  of  ¥(K)  or  MQQ  for  K  >  Kpea(c  is  <3> - 3.1.  For  wave  slope  spectra 

the  slope  is  of  course  2  larger  than  for  wave  height  spectra,  i.e. 

<p>  — 1.7.  Thus  the  SAR  spectral  slope  for  K  >  Kpg^  is  approximately 
midway  between  the  average  slopes  for  waveheight  and  wave  slope  specta. 

Overall  the  SAR  spectra  correspond  slightly  better  with  the  waveheight  spectra 
than  with  the  waveslope  spectra. 

Both  of  the  spectral  features  described  above  can  be  explained,  at  least 
qualitatively,  by  noting  that  in  equation  (2)  there  is  a  transfer  function 
j R | ^  related  to  the  radar  wave  -  ocean  wave  interaction  (Imaging)  mechanism. 
Work  by  Alpers,  Ross  and  Rufenach  (1981),  discussed  In  section  IV  above. 
Indicated  that  we  can  expect  |R|^  to  be  proportional  to  to  K^.  Consider 
the  curves  for  'f(K),  or  In  Figs.  7  and  11.  If  one  were  to  apply  the 

transfer  function  J R ) ^  to  these  curves  (as  Indicated  by  eq.  2),  we  would 
expect  the  slope  of  the  resulting  SAR  spectrum  to  be  relatively  more  steep  on 
the  low  wavenumber  side  of  the  peak  (K  <  Kpea(()  and  relatively  less  steep  on 
the  high  wavenumber  side  of  the  peak  (K  >  Kpeak)  and  this  is  Indeed  what  we 


generally  find  in  the  SEASAT-JASIN  data  of  Figs.  7  and  11.  Further  the 

changes  in  slope  from  ¥(K)  to  |FR|  induced  by  (Rp  could  cause  the  peak 

wave  number  to  be  shifted  toward  relatively  lower  wavenumber  in  SAR  spectra  as 

is  observed.  However,  the^change  in  slope  from  ¥(K)  to  | FR (K ) |  is  not  as 

2  2  3 

dramatic  as  one  would  expect  from  |R|  ~  K  or  K  .  Thus  we  conclude  that  the 
theory  of  Alpers,  Ross  and  Rufenach  is  qualitatively,  though  not 
quantitatively ,  suppor+ed  bj  the  JASIN  data  set.  Since  the  transfer  functions 
derived  by  Alpers,  et__a1_.  are  linear  and  apply  rigorously  only  to 
monochromatic  waves,  quantitative  agreement  is  not  necessarily  expected  for  a 
spectrum  of  real  ocean  waves  in  which  non-linear  phenomena  may  play  an 
important  role.  Further  theoretical  development  along  the  lines  pursued  by 
Alpers,  Ross  and  Rufenach  together  with  emperical  modifications  where 
necessary  should  lead  to  better  algorithms  for  estimating  ^(K.e)  from  SAR 
Imagery. 


VIII.  SUMMARY  AND  CONCLUSIONS 

The  principal  results  and  conclusions  of  the  SEASAT-JASIN  experiment  are 
summarized  below.  The  results  are  based  on  optically  processed  SAR  imagery 
collected  by  the  SEASAT  satellite  and  surface  measurements  collected  by  buoys, 
ships  and  aircraft. 

1.  Wave  Visibility  in  SAR  Images:  Waves  observed  by  surface  buoys  are 
sometimes,  but  not  always,  Imaged  by  the  SEASAT  SAR.  Three  factors  appear  to 
be  of  primary  importance  for  waves  to  be  visible: 

a.  SAR  system  resolution  should  be  good  enough  to  resolve  at  least 
two  SAR  pixels  per  ocean  wavelength  along  either  the  along  track  (azimuthal) 
or  cross  track  (range)  directions.  SAR  resolution  is  adversely  affected  by 

ocean  wave  otion,  especially  resolution  along  the  azimuth  direction. 
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b.  Significant  waveheight  should  be  above  some  threshold.  The 
data  examined  here  are  insufficient  to  determine  this  threshold,  but  suggest 
it  is  less  than  or  approximately  equal  to  1  m. 

c.  Wind  speed  should  be  above  some  threshold.  In  the  data  examined 
here  wind  speeds  do  not  fall  below  3.5  ms-1  and  in  the  3.5  ms-1  case  waves 
were  visible  in  the  SAR  image.  Thus  the  threshold  must  be  below  3.5  ms”1, 
probably  a  few  ms-1. 

2.  SAR  Based  Estimates  of  Dominant  Wavelength  and  Direction:  Unadjusted 
SAR  estimates  of  dominant  wavelength  and  direction  for  the  SEASAT-JASIN 
experiment  are  accurate  to  about  +13%  and  ±10°  respectively.  In  computing  the 
latter  figure  approximate  buoy  measurements  (marked  -  in  Table  1)  were 
excluded  and  thus  only  4  JASIN  cases  were  used.  The  SAR  estimates  of  dominant 
wavelength  (for  the  JASIN  experiment  alone)  are  biased  high  by  about  10%.  If 
this  bias  were  to  be  removed  by  a  linear  least  mean  squares  fit  the  dominant 
wavelength  could  be  estimated  to  about  +  10%  using  SAR  data.  The  180° 
ambiguity  usually  present  in  SAR  estimates  of  wave  direction  can  be  resolved 
by  focussing  tests  during  the  imaging  process.  These  tests  remove  the  need 
for  ad  hoc  assumptions  such  as  waves  always  traveling  downwind  or  toward 
coasts. 


3.  Estimates  of  H|  q  Using  SAR  Images:  The  signal  to  noise  ration  (S/N) 
or  peak  to  background  ratio  (PBR)  at  the  dominant  wavenumber  of  a  SAR  Image  K 
spectrum  is  correlated  with  the  significant  wave  height  (R^  »  0.53).  On 
the  basis  of  the  data  reported  here  we  suggest  that  can  be  measured  to  an 
average  accuracy  of  about  tl  m  using  Fourier  spectra  of  SAR  images.  We  note, 
however,  that  in  the  JASIN  data  set  wavelength  is  also  correlated  to  (S/N), 
though  lass  ^  ■  0.40). 
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4. 


Estimates  of  Directional  Waveheight  Spectra  Using  SAR  linage 

Spectra:  SAR  image  spectra  are  in  rough  agreement  with  buoy  measurements  of 

omnidirectional  ocean  waveheight  spectra  *(K)  and  correspond  less  closely  with 

ocean  wave  slope  spectra  ¥'(K).  However,  the  dominant  wavenumber  in  a  SAR 

spectrum  typically  falls  below  the  dominant  wavenumber  of  a  buoy  spectrum. 

Further,  the  typical  slope  of  the  SAR  spectra  on  the  high  wavenumber  side  of 

the  peak  (K  >  Kp^)  is  ~  K  *  while  the  corresponding  slopes  for  buoy 
-3  7  -11 

spectra  are  ~  K  *  and  ~  K  ’  for  waveheight  and  wave  slope 
respectively.  These  differences  between  SAR  and  buoy  spectra  are  in 
qualitative  agreement  with  the  theory  of  Alpers,  Ross  and  Rufenach  (1981). 
Analysis  algorithms  based  on  this  theory,  but  including  emperical 
modifications  should  substantially  improve  estimates  of  *(K,9)  using  SAR 
images. 
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itlintk  11  mtrh  roll  buov :  b.  Illll  Have  rider  buoys  moored  in  0IA--see  fig.  1;  t.  Discovery  Pitch  roll  buoy 


triangles  use  only  buoy  data,  while  the  solid  triangles  indicate  that  a 
SAR  estimate  for  <&  was  used  in  the  absence  of  buoy  data. 

Fig.  5.  Optical  Fourier  transform  (OFT)  of  an  ocean  wave  field  from  orbit  547 
near  the  FIA  area  of  Fig.  1.  The  photographic  display  is  in  the  form  of  a 
polar  plot  with  wavenumber  K  measured  linearly  from  the  bright  central 
spot.  The  angular  coordinate  in  the  display  corresponds  to  wave 
direction.  Note  the  180°  ambiguity  in  direction  and  the  two  directional 
components.  The  bright  radial  lines  are  due  to  the  rectangular  aperture 
used  with  the  lens  system. 

Fig.  6.  Digital  Fourier  transform  (DFT)  of  an  ocean  wave  field  from  orbit  547 
near  the  FIA  area  of  Fig.  1.  The  contour  plot  is  in  polar  coordinates 
with  the  Kx  and  Ky  axes  corresponding  to  the  directions  shown  in  Fig.  3. 
This  DFT  was  computed  from  a  portion  of  the  same  image  used  to  form  the 
OFT  in  Fig.  5.  Note  the  two  directional  components  and  the  180° 
directional  ambiguity. 

Fig.  7.  Comparison  of  SAR  image  spectrum  (corrected  for  background 
fluctuations  and  SAR  point  response)  IF^K)!6  with  omnidirectional 
waveheight  and  waveslope  spectra  Y(k)  and  Y'(K)  as  measured  by  pitch-roll 
buoy.  Both  |FR(K)|^  and  Y'(K)  have  been  normalized  along  the  ordinate  to 
the  peak  of  Y(K). 


Fig.  8.  Comparison  of  dominant  ocean  wavelength  as  measured  by  SEASAT  SAR  ana 
surface  buoys.  The  SAR  measurements  are  biased  slightly  high  with  respect 
to  the  surface  measurements.  The  average  percent  difference  between  the 


* 
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Fig.  9.  Comparison  of  dominant  ocean  wave  direction  as  measured  by  SEASAT  SAR 
and  surface  pitch-roll  buoys.  The  SAR  measurements  show  no  significant 
bias  with  respect  to  the  surface  measurements.  The  mean  difference 
between  the  two  measurements  is  +U°.  If  two  JASIN  cases  involving 
approximate  buoy  measurements  (labeled  -  in  Table  1)  are  omitted,  the  mean 
difference  is  reduced  to  +ic°. 

Fig.  10.  Comparison  of  SEASAT  SAR  image  (S/N)  at  the  peak  of  the  K  spectrum 
with  buoy  measurements  of  H^j.  The  two  quantities  are  correlated  with  a 
determination  coefficient  R  *  0.54.  A  linear  least  mean  squares  fit  to 
the  data  is  shown  by  the  dotted  line.  Using  the  linear  fit  could  be 

measured  to  within  ±1  m  using  SAR  image  data.  Note,  however,  that 
wavelength  also  correlates  with  (S/N),  though  less  strongly  (R2  =  0.40). 

Fig.  11.  Comparisons  of  SAR  image  spectra  (corrected  for  background 
fluctuations  and  SAR  point  response)  |FR(K)|2  with  omnidirectional 
waveheight  and  waveslope  spectra  'if ( k )  (or  Moo)  and  ¥'(K)  (or  K  Mqo)  as 
measured  by  pitch-roll  buoys.  Both  |FR(K)|2  and  ¥'(K)  have  been 
normalized  along  the  ordinate  at  the  peak  of  ¥(K).  Note  that  the  peak  of 
| FR (K )  | 2  is  typically  displaced  toward  smaller  K  relative  to  the  peak 
of  ¥(K).  Further  notice  that  the  slope  of  |FR(K)|2  below  the  peak 
(smaller  K)  is  typically  more  steep  than  the  slope  of  ¥(K)  and  that  the 
slope  of  1 FR (K ) | 2  above  the  peak  (larger  K)  is  typically  less  steep  than 
the  slope  of  ¥(K).  These  six  comparisons  are  from  orbits  556  (a),  762  (b) 
&  (c),  791  (d),  1044  (e),  1049  (f).  Further  details  are  given  in  Table 
I.  For  orbit  762  (b)  refers  to  the  upper  entry  and  (c)  to  the  lower  entry 
in  Table  I. 
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J.  F.  Vosacky  and  H.  ft.  Assal 
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The  SEASAT  satelllt*  system.  up  Ins  «n  onboird  73 
ch  wavelength  synthetic  aperture  ndar  (SAP). 
coIlKtld  -  7S  to  *0  n  resolution  radar  load**  of 
land  and  oeaait.  w«  sunmerlit  multi  froai  IS 
SEASAT  SAP  pastil  during  thn  Joint  Air-Sea 
Intiractton  (JASIH)  esperlmtnt  conduct id  off  Che 
unit  coast  of  Scotland  In  suamer  1978.  Rrailty 
waves  of  lonotf!  -  80  to  3011  m  ere  ofttn  toon  In 
SAP  fnagot:  hoMxr.  Hitting  waves  an  not  always 
soon.  Hi  find  that  weveheloht.  SAP  molutlon  and 
Mind  velocity  an  Important  crltorla  for  deter- 
mining  navi  visibility.  Buoy  mosunaiints  of  H, 
corn  lata  uitb  poak  signal  to  nolsn  ratio  in 
four  1m  transform  of  SAP  Ins  ms.  Althouoh  thin 
an  significant  dlffinneis.  SAP  rttlnatts  of 
omnidirectional  Mtvoholght  spoctra  an  In  rough 
agrwmnt  ulth  buoy  spoctn.  Feature*  ms  won  no 
lot  mil  non  ppttons  an  almost  always  found 
ovar  seamounts  and  at  thi  »don  of  thn  continental 
shilf.  Ships  and  thotr  wakns  an  ortm,  hut  not 
always,  dotictabli  In  SEASAT  SAP  Imgns. 

tiywords:  Padlo  Ocoonography,  Remote  Sonslng. 

Gravity  navis.  Ships.  Hokes 


U  nmwouaion 

Tim  tynthitle  apirtun  radar  (SAP)  carrlid 
aboard  tho  SEASAT  fatal  11 ta  coll  act ad  radar  toagis 
of  tho  ocoon  In  IWI  k*  wida  swaths  ranglno  In 
lnngth  froai  about  300  to  3000  ka.  hurl  no  thn 
JASIH  liporlagnt  (Juna-Saptaabar.  1978)  surraci 
obsorvitlons  of  agny  air  and  sia  paraaatirs  win 
mdo  hy  ships,  buoys  and  aircraft,  a.g.  thi  wind 
and  wovi  conditions  which  an  of  nost  Intamt 
ban  (Pif.l).  During  August  and  Saptaahor.  1978, 
18  SEASAT  passos  producid  SAP  laogory  of  tho  JASIR 
ana.  Eaaaplas  of  such  passos  an  lllultratod  In 
Fig.  I.  Thiso  nltlpli  SEASAT  obsorvitlons 
envoi H  with  mtonstv*  surfaco  anasunamts  pro- 
vldl  a  unloua  opportunity  to  assay s  tha  nenan 
noasurawnf  eaoahilltias  of  satiltlti  Synthgtlc 
aportun  radar. 

SAP  Imagery  Is  capahli  (undor  sultabli  condi¬ 
tions)  of  sonslno  a  largo  vanity  of  own  surfaco 
phnnanna  including  long  gravity  want  (with  wan. 
imoth.  a  ,  gnattr  than  about  twlei  tho  SAP  riso- 
lvtlnn) .  winds,  furfaca  currant! .  Intimal  want . 
nenn  fronts,  undrrwatnr  typography,  slicks,  and 
ships  and  thoir  wskos  (Pgfs.  7.3).  Thg  SN  Inset 
•s  basically  a  high  msolutfon  two  of  tho  radar 


riflictlng  propartlas  of  the  ncaan  surface  modi- 
flad  by  surfaci  notion  ifficts.  Slnci  til  flielec- 
trlc  propirtlis  of  thi  ocian  surfaci  am  mlatlvi- 
ly  uniform,  variations  In  Surfaci  roughness  are 
prlnerf ly  responsible  for  changes  in  the  moar 
backscatter  strength.  In  particular  'resonant' 
ocean  waves  with  A  -  30  cm  Interacting  with  the  73 
cm  radar  waves  via  a  ’ragg  scatter  recnanlsm 
produce  the  dominant  radar  echo.  because  SAP 


Oi*  eeiMOinanmc  vim-  eowrowva 

inrtnsiva  amt*  wt  <  >i»i« 


Fig.  1.  Joint  Air-Sea  Interaction  (JAMR)  •.reel- 
nent  area  showing  typical  SFASAT  tee  inace 
Swaths  and  tactions  of  images  ns»a  for  wav* 
anolylls.  Have  rider  buoys  are  nenoten  h» 
CE,  CJ.  etc.  Pitch-roll  buoys  were  oe-rra 1 ' y 
deployed  in  the  F(A  area. 
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utilizes  the  ptiu  o f  the  radar  rcbo  to  achieve 
high  resolution  In  the  direction  parallel  to  the 
flight  path.  i.t.  to  correctly  locate  the  origin 
of  a  given  echo,  notion  of  tho  ocean  surface  can 
cauao  echo  energy  to  he  eltplaceO  In  the  laag*. 
Thia  misplacement  Is  Illustrate*  In  Fig.  a  where  a 
noving  ship  Is  displaced  from  its  wake.  Spatial 
■adulation  of  surface  roughness  and  notion  (giving 
rise  to  Doppler  shifts)  produces  the  najorlty  of 
ocean  surface  features  In  SM  Imagery.  Phillips 
(<)  describes  a  number  of  ways  In  which  wind,  wave 
and  current  variations  can  modulate  decimeter 
scale  surface  roughness. 

2.  GRAVITY  wave  ORSERVATIOftS 

2.1  have  Visibility  In  SAR  Images 

although  gravity  waves  of  length  «  80  to  300 
m  are  often  seen  In  SAP  imagery,  we  find  that 
waves  observed  by  surface  buoys  are  not  always 
visible.  In  ou r  study  we  have  used  only  imagery 
optically  processed  at  Jet  Propulsion  Laboratory 
(JPL)  and  Environmental  Research  Institute  of 
Michigan  (ERIh).  Study  of  JASIN  data  Indicates 
three  criteria  are  Important  in  determining  wave 
visibility:  wave  height.  SAP  resolution  and  wind 
velocity.  Alpert,  Ross  and  Rufenach  (5)  point  out 
that  the  orbital  velocities  associated  with  sur¬ 
face  waves  produce  a  smearing  effect  on  SAR 
imaqery  increasing  resolution  cell  site  especially 
along  the  atlmuth  direction,  (parallel  te  the  SAR 
flight  pathX.  An  expression  for  the  degraded 
resolution  cell  length  Ay*  along  the  atlmuth 
direction  Is  given  In  Ref.  $.  Knowing  by*  and  the 
ocean  wavelength  along  the  atlmuth  direction  Ay  • 
A/cosO.  the  number  of  samples  per  wavelength' n 
along  the  atlmuth  direction  can  be  calculated; 
Here  o  Is  the  angle  between  the  spacecraft  veloc¬ 
ity  vector  V  and  the  ocean  wave  vector  K.  The 
hyguttt  sampling  criterion  demands  n  »  2  In  order 
that  the  waves  be  adeguately  sampled.  Jain  (6) 
and  Vesecky  et  al.  (7)  have  developed  wave  visi¬ 
bility  crlterTa "JTong  these  lines,  tacause  of  the 
roles  played  by  waveholght  and/or  slope  In  wave 
Imaging  mechanisms  (Ref.  S)  one  aspects  wave  visi¬ 
bility  to  decrease  with  decreasing  significant 
waveholght  H.n  and  Indeed  we  find  waves  are 
seldom  Imaged  for  H...  <  1  m.  Applying  these  two 
criteria  n  >.2  and  h./j  l  1  "  to  the  JASIh  data 
set  we  found  only  one  case  out  of  IS  In  which 
these  criteria  did  not  work— nine  cates  where 
waves  were  visible  In  SAP  Images  and  sis  when 
waves  were  not  observed  in  SAP  Images.  All  ocean 
wave  data  were  drawn  from  buoy  measurements. 

Since  wind  raises  the  3S  cm  ocean  waves  which 
tho  radar  senses,  one  aspects  reduced  wave  visi¬ 
bility  at  low  wind  speeds.  During  the  SEASAT- 
JASIh  esperlment  waves  were  detectable  ever  a 
range  of  wind  speeds  of  from  3.S  to  IS. 2  m/t  while 
for  wind  speeds  of  from  S.S  te  12.9  m/s  waves  war* 
not  detected.  It  it  thus  apparent  that  wind  spend 
was  not  a  factor  during  the  SEASAT -JASIh  esperl¬ 
ment  mainly  because  mind  speeds  were  always  above 
J.S  m/s. 

z.2  jgwriwns^dMAP  and  Pupy  Estimates  of  wave 

2.2.1  wave  measurement  wet hod.  The  wave  measure¬ 
ments  from  sap  imagery  made  "ere  art  based  on  * 
simple  assumption,  te  wit.  that  fluctuations  In 
SAP  Image  Intensity  t(s.y)  ere  proportional  to 
ocean  wave  height  fluctuations.  The  (s.y)  plane 
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is  the  plane  of  tn*  avnranr  ocean  surface.  Under 
this  assumption  th»  directional  ocean  waveheight 
spectrum  »(K,0)  is  simply  proportional  to  the  mag¬ 
nitude  squared  of  the  two-dimensional  Fourier 
transform  ot  the  Intensity  fluctuations,  namely 
|F(I  («.y)  )|  .  Here  It  Is  2*  divided  by  the  ocean 
wavelength  and  0  Is  direction  in  the  («.y)  plane 
(0  is  usually  degrees  clockwise  from  true  north). 
As  will  become  clear  below,  this  assumption  Is  not 
precisely  correct.  However.  It  is  apparently  n0t 
far  wrong  (Ref.  8)  and  provides  a  convenient 
working  hypothesis.  Although  this  alooritnm  is 
straight-forward,  there  are  numerous  variations  to 
its  implementation.  Here  we  have  used  ootically 
Imaged  data  produced  at  the  Jet.  Propulsion 
Laboratory  (JPL)  and  at  the  Environmental  Research 
Institute  of  Wlehtgen  (ERIH).  finer  the  Imagery  is 
Fourier  transformed,  estimates  of  the  dominant 
wavelength  and  direction  can  be  made  on  the  basis 
of  our -assumption  by  simply  finding  the  peak  of 
|F(I)|'  and  noting  the  location  In  terms  of  wave¬ 
number  K  and  direction  3.  The  18D*  ambiguity  >n  3 
can  be  resolved  by  special  focussing  measurements 
during  the  optical  imaging  process.  Vesecky  et 
al.  (7)  note  that  the  peak  slanal  to  noise  power 
ratio  (S/K)  In  the  Fourier  transform  of  a  SAR 
Image  is  correlated  (r  -  0.8)  with  significant 
waveheloht  H,/3.  Peak  S/h  is  estimated  in  Fig.  2 
by  noting  tne  level  regions  in  IF|Z  at  very  low 
(K  1(1.015  m*‘)  and  very  high  (K  i  0.05  n*1)  wave- 
nuabert  and  averaging  the  two  levels  to  obtain  a 
background  noise  level,  with  this  method  Hj  ^  can 
probably  be  measured  to  t3/A  m  using  this  col-rela¬ 
tion. 


2.2.2  SAR  and  buoy  Estimates  of  Dominant  wave¬ 
length  and  Direction.  Eleven  estimates  by  sHI  and 
surface  buoys  during  the  JASIh  experiment  were 
compared.  The  average  distance  between  the  center 
of  the  SAR  Image  and  the  buoy  was  -  67  km  with  a 
range  of  0  to  1A6  km.  The  average  percent  differ¬ 
ence  In  wavelength  was  15X  with  a  consistent  bias 
nf  SAP  estimates  toward  longer  wavelengtns  by 
about  lOt.  Six  SAP  and  buoy  estimates  of  dominant 
wave  direction  were  compared  with  an  average  error 
of  13*.  However,  two  pitch-roll  buoy  estimates 
were  only  approximate  due  to  compass  problems.  If 
we  exclude  these  two  approximate  estimates,  the 
average  error  was  only  7*.  There  was  no  bias 
relative  to  the  range  direction.  I.e.  perpendicu¬ 
lar  to  the  subsatelllte  track. 

2.2.3  SAP  and  huov  Estimates  of  Omnidirectional 
MaveheldSt  Spectrum?  during  tne  JASIH  experiment 
estimates  of  the  omnidirectional  wave  height  spec¬ 
trum  H„  from  buoys  And  from  SAR  images  were  com¬ 
pared  Th  seven  cases.  Fig.  2  shows  a  typical 
example  of  theta  comparisons.  The  buoy  estimates 
wer*  computed  by  standard  methods  (Ref.  «)  while 
the  SAR  ettlmetes  are  slmoly  IF(I)|Z  ■  !F(k,s)|2 
averaged  over  a  30*  angular  sector  around  the 
dominant  wave  direction  and  nomalitnd  along  the 
ordinate  to  the  peak  of  H--.  from  buoy  measure¬ 
ments.  ho  normelltatlon  was  done  alona  tn*  ab¬ 
scissa.  Since  in  these  cases  there  is  only  , 
single  dominant  wave  direction,  data  from  direc¬ 
tions  away  from  the  dominant  wave  «irwc*-nn  con¬ 
tains  only  background  fluctuations:  its  ‘"elusion 
would  only  serve  te  raise  the  noise  level  in  the 
SAR  estimate  Of  HM. 
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?.2.4  IH f f effnce S  Petween  ?»  and  »ufflr  EsMmat»S 
of  t*nft  and  therr~ffi>T7natTqn  )?***?- "?QC*!*Al-.n 


inwqry.  r,«.  r  4  number'  orVeetumx 

TrtiTcE  are  typical  Of  comparisons  het ween  SAB  image 
Intensity  spectra  |F(K)|Z  •«*  buoy  wiiurwfltj  of 
ocean  wave  heioM  end  wave  slope  spectra.  x  *no* 
The  peak  of  the  SAP  jpnc try*"  is  -  3  times  t*e 
hackornund  indicated  hy  the  horlanntal  Ho#. 
Typically  this  would  correspond  to  Hl/3  -  3  m. 
Tht  peak  Of  |F(K)[Z  corresponds  to  «  wavelendth 
so**  10?  above  the  wavelength  of  tne  peak  «00  and 
20?  above  the  wiv#le«ath  of  the  neak  * **«*- 
ured  by  surface  hudys.  The  shaoe  of  the  >5"  spec* 
tru*»  corresponds  most  closely  with  the  wavehe»ght 

spectrum.  however,  at  wavenumbers.aDcve  the  peak. 

the  slope  of  the  SAP  spectrum  closer  to 

the  slope  Of  the  wave  slope  spectrum  thin 

to  that  of  the  wavenetght  spectrum  K"  *  .  This 
behavior  of  the  slope  corresponds  roughly  to  the 


Fig.  ?.  ^rwpart sen  of  SAP  image  spectra  (cor¬ 
rected  for  background  fluctuations  and  ?ap 
point  response)  with  omnidirectional  wave- 
hflnht  and  waveslope  Spec tra  r^00(K) 

respectively)  measured  by  moored  w* ve¬ 
nd#?  huoy.C2  m  Flo.  3.  *nth  the  $A»  sree- 
true  and  <zh_(J  h*ve  been  norm#  H  red  a»ono  the 
ordinate  at  ™h#  peak  of  *n<v  The  cap  and 
huOy  measurements  were  separated  hy  some  260 
km  alono  the  surfac*  and  some  A  hours  m 
time,  hut  under  stable  conditions  (Awe.  IP 
and  IP.  197P).  This  comparison  is  used  here 
because  it  illustrates  typical  features  so 
well,  bote  that  the  peak  of  the  SAP  spectrum 
is  displaced  toward  smaller  k  with  respect  to 
the  peaks  Of  both  *  and  farther 

notice  that  the  slope  nf  t*e  vP  spectrum 
heioy  the  oeak  (smaller  r)  M  mnre  steer  than 
the  slope  of  *00  and  that  the  slope  of  tee 
tAP  spectrum  above  t*e  rea*  Mercer  *M  is 

'ess  Sfeep  than  the  slnne  nf  v  . . 


theory  disc  iS***0  hy  At  nets  »t  el.  (r )  wh»-rr  they 
predict  that  the  spy  spectrum  slorn  «hn*«M  *o 
•  r.-  **?  to  V’  r  whwrr  the  *»av**h.» inht  spectrum 

slope  is  -  K’\  This  theory  also  presets  tKe 
Mas  Of  the  SAP  estimate  of  the  peak  of  M*0  toward^ 
longer  wavnlencths  which  we  observe. 

3.  lAPOF-SfAU  FfATlIMtS  hi  ip  IBP. 

T«F  SRASAT-JASIM  F*PCa|f**T 

3.1  3  •  r  s  i  s  t ence  of  Larce-*ca»e  features 

One  of  the  oreat  advantaces  of  SAP  iftiaoery  as 
a  remote  sens  •  no  tool  is  that  surface  r-orpro'  oey 
can  he  examined  over  larce  areas  (inn's  to  I'vtO  s 
rf  km)  in  virtually  a  snap  shot  'ashior.  <as 
imaoes  of  the  ocean  often  shew  'natures.  I  tc  *0 
km  In  extent  which  appear  to  he  internal  waves, 
oceanic  fronts,  slicks  and  spatial  variability 
(mistiness)  of  wind  Speer.  Various  examples  nf? 
the  it.?,  east  coast  are  niven  hv  "ea1  *T  aj .  n)- 
Visual  examination  of  SFASAT  SAB  imaoerv 
(1:1.000.00(1  scale)  collected  over  the  .1AM*  ex¬ 
periment  area  revealed  a  number  of  interest’nc 
characteristics. 

when  the  same  area  **as  imaoed  t**ce  within  a 
relatively  short  period,  no  hoM  larce  scale  fea¬ 
tures  were  found  to  persist  over  time  sca’es  of  * 

hours  (A  inace  pairs)  and  15  hours  '3  imaoe 

pa<*s).  Store  the  •~tane  oai^S  COu>  ’  rot  *>0  pre¬ 
cisely  aliened,  s*#1!  weak  features  could  nave 
persisted,  hut  been  overt  oo*e*i.  tr.  four  cases 
land  at  sore  point  could  he  used  to  raster  the 
f manes  to  within  a  few  kilometers.  To  the  other 
three  no  land  was  eresent  in  t*e  imsoery  eno 

superposition  was  less  accurate.  >e  4  orbits 
which  produced  these  pal*-*  of  i maces  dig  not  Have 
parallel  ground  tracks  anc  so  tne  boundaries  of 
the  imaoe  swaths  intersect  with  ancles  of  jhout 
no*.  Fig.  I  cortains  an  example  involve  orbits 
547  A  55ft.  Thus  the  viewinc  oeemet ry  d’ffyrs  Ctn- 
slderahly  between  the  two  i^aoes  in  a  ram.  Prom 
these  observations  we  conc’ude  that  eit**e»-  thy 
features  were  not  persistent  and  **nce  rot  ocean ic 
in  origin  or  that  the  conditions  for  sensiro  the 
phenomenon  did  not  persist. 

3 . 2  Internal  wave  Peat-jres 

Features  resemfliro  tne  internal  wav* 
patterns  of  iAN0$AT  visual  imagery  '®e'.  lb;  were 
almost  always  'ound  over  seameurts  and  at  the  »doe 
of  the  continental  shelf.  *he  patterns  «*•  erne«- 
eral  probably  deperainc  c«  t»*e  ttace  0'  t*e  tme. 
Fig.  3  Illustrates  such  patterns  occ.ir^a  tver  !"* 
steep  bathymetric  oradients  associated  w*th  t*e 
nyviu e-Thomson  3iflne  wr  ?0'  •».  «rv.  n^crr,s 
ranee  from  -  1*10  to  500  m  ‘r  this  rec'nr. 

3 . 3  *Hnd  Poll  vortices 

In  Imtanry  *r(y«  prM*  10?A.  rc*wj  of  f*reeij 
separated  hy  either  1.3  nr  ?.*  ;m  are  seen  rver  an 
area  of  at  least  1^0  x  ?t0  **.»,•  stre*.  s 

aliened  along  the  wind  direction  where  t»*e  m'^C 

speed  was  aenroximateW  12  */s.  t*»t 

these  chanoes  4"  fao  <mace  hr4Oh'r*sS  nv«r  j  iarc» 
area  were  due  *o  chann*s  m  wiM  velocity  >  "m- 
sonable  assumotim.  viz.  ®ef.  11'  co^c'une  ?h8- 
this  is  th#  first  clear  rhservaf^nn  z4  r cl* 
tlces  in  the  surface  Pc u« nary  layer  fvew  •»* 
ocean. 


Copt  wfl obte  te  PT1C  4oe»  ** 
pwmtt  fuBr  hgW»  »P»A«0,hm 
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3.4  Damping  Of  30  cm  kurfoce  Waves  by  $am 

Ox  orhit  $56  an  ami  roughly  SO  km  lx  sijo 
app eared  very  dark  (very  low  radar  ref tectivity). 
Several  tourers  of  data  IncluOlng  ship  loot. 
Infrarod  satellite  imagery  axd  weather  tyxoptxt 
indicated  that  thunderstorms  were  likely.  ue  con¬ 
clude  that  raix  in  (Hit  area  prohaOIy  smoothed  the 
tea  on  teala  titet  of  ~  30  cm  (to  which  the  radar 
signal  is  sansitlve)  thus  causing  the  dark  (nature 
«nen  in  the  imagery. 


4.  cm jut 

A  perusal  of  SEASAT  SAP  imagery  (e.g.  sen 
Pef.  2  or  Pig.  4)  reveals  many  imaoet  of  station¬ 
ary  ships  as  «elt  as  moving  ships  with  their 
attendant  wakes,  which  can  extend  -  20  km.  In 
only  a  few  cases  have  searches  heen  made  for  ships 
with  known  positions.  He  have  searched  for  ocean¬ 
ographic  ships  involved  <n  the  JASth  experiment  on 
orbits  $47  and  $90  using  imagery  which  was  opti¬ 
cally  processed  at  JPl.  The  ships  involved  and 
their  approximate  lengths  am;  P.v.  Attantts  It 
(64  n),  Gardline  Fndurer  ($}  a).  9.V.  meteor 

(62  m)  and  H.R.  m.s.  Tydgman  (A9  m).  The  ship 
1n.-a"nns  in  latitude  and  Inngi-.jdn  ef  rhn  times 


fig.  1.  SEASAT  SAP  Image  from  orbit  7*2  near 
lM4*;]f  ITT  on  Aug.  t«.  197*.  The  small  white 
dots  are  about  7.2  km  apart,  we  interpret 
the  wove  pockets,  most  prominent  at  center 
and  loner  left  center .  as  surface  roughness 
modulation  atsnciatnd  with  the  surface  cur¬ 
rent.  dlveruence  produced  hy  internal  waves. 
The  region  imgood  Is  centered  Over  the  steep 
bathymetric  gradients  nf  the  uyvtUe-Thgnsnn 
otdge  near  W  *»'  m  and  *'  ».  ikegths  are  * 

a»t  a**.*-  -*  *»-** 

'.*+  •  V  I  .*  %.#-»..«  Mb  .*•  •  H#  w«V*  *+C9’ 

.•  •,.»  *  .  mojbf  #1^  ’Aenjimil,  |l|a 

hathyawtric  nradimet  hnlnw.  The  arrow  below 
*!>•  (mane  indicates  true  north  and  the  tears 
was  «pti«*lly  prncmsend  at  th#  jet  frgpolsinn 
lahorernry. 


of  SFASAT  5Mt  imanrry  are  known  tn  an  accuracy  of 
a  few  ►"  nr  less,  witnout  special  processing  *e 
are  usually  ahle  to  locate  points  («  SEASAT 
imaoery  to  within  ahnuk  5  km.  Atlnwinn  an  addi 
Monai  ip  kn  for  unknown  systematic  errnr.  we 
eipect  fO  find  ship  images  within  15  kn  of  their 
estimated  position  in  the  imanrry. 

(In  orhtt  $47  only  the  *tient's  II  was  within 
the  SA»  swath  and  a  ship  imeet  was  found  12  km 

from  the  estimated  position  of  the  (mace,  hone  of 
tho  Other  JASIh  ships  mentioned  aoove  were  within 
the  SAP  swath  and  no  other  ship  manes  »»re  found 
within  1  IPO  km  along  the  sane  2$  km  wine  sub- 
swath.  hear  the  Atlantis  II  position  the  dominant 
smell  wes  of  wevelenoth  .1  -  17n  m  out  of  2$’  true 
with  H.  '  1.6  n  woile  the  wind  wes  out  of  IS*  et 
3.6  m/s. 

On  orbit  $90  all  four-  ships  were  within  the 
SAP  swath.  However,  only  the  largest  ship  (H.t. 

".$.  Tydemen)  can  even  be  tentatively  identified 
—  a  weak  target  some  7  km  from  the  estimated 

position.  hear  the  positions  of  a'l  tn»  shios, 
save  the  Gardline  Endurer.  the  dominant  swell  was 
A  -  222  m  out  of  2S0C  true  with  H,  .  -  1.2  m  unite 
the  wind  was  out  of  10*  at  (.1  «s-1.  *t  the 

position  of  the  Gardl  ine  Endurer  some  ISO  km  to 
the  hw  tho  wave  conditions  were  aoprpximately  the 
same  wnile  the  wind  speed  was  S.2  ms-1. 

This  evidence,  taken  with  an  unsuccessful 
search  for  thr  p.v.  nceenooraoher  during  the 
GOASFX  experiment,  leads  us  to  conclude  that 
•maces  Of  stationary  or  slowly  moving  ships  (no 
significant  mates)  are  fiff*rult  n r  impossible  (n 
Intepr  In  cr.ilAT  1  naneny  ,'nptlca'lv  procexsnn  at 
.*•1.)  fur  mortmratn  seas  ann  wing  velocities  Imater 
than  about  l.S  "s"  .  Thus  sucn  SAP  imaomry  cannot 
he  routinely  used  CO  determine  shio  Positions  as 
an  aid  to  performing  research  at  see.  re  e»ohe- 
stte  that  the  imagery  used  in  these  searches  »es 
the  steederd  'oroducc*nn  nin'  date  oroduct. 
Improvements  in  resolution  and  focussing  by  either 
soocialited  optical  or  digital  imeoing  would  pre¬ 
sumably  improve  detection  performance.  e'a.  A  ‘S 
an  example  of  hiob  aulllty  figitolly  proc»ss»o 
imagery.  “any  ships  are  evident.  *he  wakes  of 
foster  moving  ships  mekt  them  more  ev'oeot.  It  is 
intortstino  to  note  mat  a  oroup  of  four  ships 
just  to  the  lower  Mght  gf  center  are  *mageo  to 
tho  loft  of  their  wakes.  Veseckv  and  Stewart  ft ? 1 
point  out  that  this  displacement  (which  Is  slmc'y 
a  result  of  ship  movement  relative  to  the  average 
ocean  surface)  can  b»  used  to  measure  shio 
voloelty.  *  ship  hovlng  a  votoc'ty  component  v, 
perpendicular  to  the  satellite  grourfl  track  causes 
the  ship  to  bo  displaced  a  aistenee  o  *  »  x/v 
penile!  to  tbe  sote'nto  voloelty  vector  but 
lh  an  opposite  direction,  wner*  x  is  the  oerpen- 
d'eular  distance  from  the  satellite  surface  tracx 
to  tbe  shio  location.  In  Fio.  1  the  satellite 
surfoce  trac*  is  oara'lei  to  t"e  top  of  the  (mace 
amp  runs  left  to  right  some  24n  km  above  the  ton 
nf  the  image.  knowing  v  v  71P0  ns*‘.  0  ■  330  " 
and  x  •  277.009  "  for  tho  ship  nearest  the  bottom 
•n  the  group  .to  the  lower  right  of  center  we  fine 
v  •  9.7  ms*  .  From  the  ship  s  waie  we  twe  that 
it’is  moviho  very  nearly  alono  toe  x-direct’nn  to 
that  it»  speed  is  *14  tts— a  reasneaMe  value. 


Fig.  4.  $£4<A7  5A3  imaoe  from  orhit  '62  near  P6-P 

II?  cn  Auq.  19,  1974.  Vertically  upward  ir 
the  f*aure  is  alone  direction  24l4  t rue.  in* 
coast  i*ne  shown  evtenos  'rom  falais  (pricnt 
region  at  too  left'  to  ‘‘'unkerque  'priori 
reoion  alona  coast  at  'ower  ieft).  fur'ac- 
toughness  variations  associated  witn  hath y- 
metric  features  in  the  Struts  of  Ocver  cause 
linear  and  curved  features  »*  tne  sad  i~ac» 
eso ec’il’y  at  lower  left  and  upwards  a’on- 
the  coast.  4t  1  Owen  ricnr  center  a  group  n* 
Several  ships  and  their  wakes  ar»  jtfihl*. 
*iote  that  the  ships  are  rtiso'aced  #*,or»  th«ir 
wakes  during  the  <AR  imaging  process  because 
Of  the  ship  velocity  (-  10  n/sl  with  respec* 
to  the  average  surface.  This  was 

digitally  processed  at  the  "crweaMn  ''•fense 
Pstahl  <shr*wnt. 
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